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ABSTRACT 


The  research  work  reported  in  this  dissertation  deals  with 
the  inelastic  behavior  of  multistory,  planar  steel  frames.  The  study 
includes  a  review  of  different  types  of  formulation  for  stability  and 
strength  analyses.  A  method  to  compute  elastic  buckling  loads  for 
multistory  frames  is  discussed. 

A  general  approach  to  the  elastic  and  inelastic  nonlinear 
analyses  of  multistory  frames  is  presented  and  a  finite  element 
formulation  is  developed.  A  technique  to  determine  the  load- 
deformation  characteristics  of  frames,  using  the  finite  element 
equations,  is  discussed.  The  inelastic  analysis  includes  the  effect 
of  axial  loads  on  the  stiffness  of  the  frame,  gradual  penetration  of 
yielding  in  the  cross  section,  and  the  spread  of  inelastic  zones  along 
the  member  length.  The  effects  of  residual  stresses  and  strain  hard¬ 
ening  have  been  taken  into  account  in  the  analysis. 

The  features  of  the  elastic-plastic  response  of  frames  for 
first  and  second  order  analyses  are  presented.  Throughout  the  study 
numerical  examples  are  given  and  compared  with  available  results. 

In  the  final  section  the  behavior  of  a  number  of  frames  is 
examined  through  inelastic  and  elastic-plastic  analyses.  Comparisons 
are  made  between  the  two  approaches  and  the  effects  of  residual 
stresses  and  strain  hardening  on  the  behavior  of  frames  are  studied. 
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CHAPTER  I 


INTRODUCTION 


Almost  all  analysis  carried  out  in  engineering  practice  is 
based  on  linear  theory.  Linear  analysis,  rather  than  nonlinear,  is 
performed  because  of  its  relative  simplicity. 

In  recent  years,  there  has  been  a  considerable  amount  of  theore¬ 
tical  work  carried  out  on  the  effects  of  geometric  and  material  non- 
linearities  on  the  behavior  of  structures.  Nonlinear  analyses  by  the 
finite  element  method  were  first  introduced  for  geometric  nonlinearity 
by  Turner  et  at.  (1960)  ,  and  for  elastic-plastic  analysis  by  Gallagher 
et  at.  (1962) .  More  recent  work  has  dealt  with  large  strains  (Hofmeister 
et  at.  1971)  ,  creep  strains  (Greenbaum  and  Rubinstein  1968) ,  dynamic 
forces  (Stricklin  et  at.  1971) ,  and  variation  of  material  properties 
with  temperature  (Marcal  1972) . 

1 . 1  Types  of  Frame  Analysis 

Several  types  of  frame  analysis  may  be  illustrated  with  reference 
to  Fig.  1-1.  This  figure  shows  the  relationship  between  the  gravity  load, 
P,  and  the  lateral  deflection  at  the  top,  A,  for  a  simple  portal  frame, 
as  predicted  by  various  approaches. 

The  simplest  form  of  frame  analysis  is  one  which  predicts  the 
elastic  critical  load  of  the  frame,  P  .  In  this  method  the  material 
is  assumed  to  be  elastic,  the  horizontal  loads  applied  to  the  frame 
are  zero,  and  all  prebuckling  deformations  (e.g.  initial  out— of— 
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Load-Deformation  Relationships 
FIGURE  1-1  Types  of  Analysis 


straightness)  are  ignored.  When  the  deformations  due  to  lateral  loads 
are  included,  a  linear  load-deflection  curve,  A,  is  produced.  This 
solution  is  usually  referred  to  as  a  ’first  order  elastic  analysis’. 
When  nonlinear  effects,  such  as  nonlinear  strain-displacement  relation¬ 
ship  and  the  effect  of  deformations  on  the  equilibrium  equations,  are 
considered,  the  formulation  results  in  a  nonlinear  load-deflection 
response,  curve  B.  This  method  is  usually  called  a  ’second  order 
elastic  analysis’.  The  elastic  analyses  give  no  indication  of  the 
ultimate  capacity  of  the  frame  or  the  true  behavior  in  the  region  of 
ultimate  load. 

If  the  material  response  of  the  frame  is  no  longer  perfectly 

elastic,  another  critical  load  may  be  reached  before  the  elastic 

buckling  load.  This  critical  load,  P  ,  results  from  the  formation  of 

P 

a  plastic  mechanism,  ignoring  all  prebuckling  deformations.  This 
method  is  referred  to  as  a  first  order  rigid-plastic  analysis,  and  is 
illustrated  by  curve  C  of  Fig.  1-1.  A  second  order  rigid-plastic 
approach  considers  the  effect  of  the  mechanism  deflections  on  the 
equilibrium  equations,  and  the  loads  therefore  must  decrease  to  maintain 
equilibrium,  as  shown  by  curve  D. 

In  first  order  elastic-plastic  analysis  the  material  is  assumed 
elastic-perf ectly  plastic,  and  prebuckling  deformations  are  considered. 
The  load-deformation  response  for  such  a  solution  is  shown  as  curve  E. 
The  upper  limit  of  the  first  order  elastic-plastic  analysis  is  the  first 
order  rigid-plastic  limit,  P^.  When  the  effects  of  sway  deformations 
on  the  equilibrium  equations  are  considered,  the  formulation  results  in 
a  second  order  elastic-plastic  analysis,  shown  as  curve  F,  in  Fig.  1-1. 
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The  ’true’  behavior  of  the  frame  is  shown  as  curve  T.  The 
difference  between  the  results  of  a  second  order  elastic-plastic  analysis 
and  the  ’true’  behavior  is  due  to  gradual  penetration  of  the  yielded 
zones,  residual  stresses,  initial  imperfection,  and  strain  hardening. 

The  magnitude  of  the  difference  depends  on  the  geometry  and  stiffness 
of  the  frame,  material  properties,  and  loading  conditions. 

1 . 2  Inelastic  Frame  Analysis 

Extensive  research  has  been  conducted  in  recent  years  to 
investigate  the  behavior  of  multistory  frames  and  their  components. 

Much  of  this  work  has  been  aimed  at  the  development  of  practical  design 
procedures  which  could  account  for  material  nonlinearity  (Majumdar  and 
Adams  1971;  Davison  and  Adams  1974). 

Several  authors  (Hodge  1959;  Mcnamee  and  Lu  1972)  have  assumed 
that  yielding  takes  place  only  at  generalized  plastic  hinges  of  zero 
length  (the  concept  of  concentrated  plasticity).  The  researchers  have 
considered  multidimensional  yield  surfaces,  with  proposed  procedures 
for  modifying  the  elastic  stiffness  of  a  frame  member  to  account  for 
yielding.  General  computational  procedures  have  not  been  developed  and 
only  simple  cases  have  been  considered.  A  few  researchers  (Moses  1964; 
Kitipornchai  and  Trahair  1975)  considered  extended  regions  of  plastic 
deformations  rather  than  discrete  hinges,  but  no  conclusions  were 
reached  regarding  their  effect  on  the  behavior  of  frames. 

Although  steel  exhibits  some  strain  hardening,  most  investiga¬ 
tors  have  assumed  that  the  material  is  ideally  elastic-plastic.  The 
influence  of  strain  hardening  on  the  ultimate  load  was  studied  by 
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Hrennikof f  (1965) ,  Sawko  (1965) ,  and  Horne  and  Medland  (1966) .  In 
most  cases  an  idealized  moment-curvature  relationship,  as  shown  in 
Fig.  1-2,  was  assumed.  However,  a  rigorous  solution  that  takes  strain 
hardening  into  account  is  at  best  a  cumbersome  and  impractical  approach. 

A  few  authors  (Alvarez  and  Birnstiel  1969)  considered  the 
effect  of  axial  load  on  the  stiffness  of  columns  and  girders.  Others 
(Korn  and  Galambos  1968)  developed  stability  functions  for  columns  on 
the  assumption  that  the  axial  forces  in  the  girders  were  small  enough 
to  neglect  their  effect  on  girder  stiffness.  Parikh  (1966)  modified 
the  column  moment-curvature  relationship  to  compensate  for  the  decrease 
in  bending  stiffness  due  to  the  yielded  condition  of  the  cross  section, 
for  axial  loads  greater  than  0.7  P  ,  where  P^  represents  the  yield  load 
of  the  column. 

The  influence  of  residual  stresses  was  considered  by  Parikh 
(1966)  and  Alvarez  and  Birnstiel  (1969) .  The  latter  formulated  the 
member  stiffness  matrix  so  as  to  account  for  the  gradual  penetration  of 
yielding,  the  presence  of  residual  stresses,  the  spread  of  inelastic 
zones  along  the  member  length,  and  strain  reversal  in  previously  yielded 
fibers.  However,  their  method  was  applicable  only  to  relatively  small 
structures . 


Extensive  experimental  work  has  been  performed  to  investigate 
the  ’true’  behavior  of  frames.  Results  of  combined  gravity  and  lateral 
loading  tests  on  large  scale  multistory  frames  have  been  reported  by 
Yura  (1965),  Yarimci  (1966),  and  Majumdar  et  at.  (1970).  Schilling 
et  at.  (1956) ,  and  Arnold  et  al.  (1968)  performed  some  tests  on  single 


story  steel  frames. 
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FIGURE  1-2  Idealized  Moment-Curvature  Relationship 
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1.3  Objective 

The  purpose  of  this  investigation  is  to  develop  a  nonlinear 
method  of  frame  analysis,  based  on  large  deformation  theory,  that  is 
applicable  to  both  elastic  and  inelastic  solutions  of  plane  frame 
problems  under  any  conditions  of  loading  and  geometry. 

In  the  inelastic  formulation  the  effect  of  axial  loads  on  the 
stiffness  of  the  structure  is  considered.  Gradual  penetration  of 
yielding,  the  spread  of  inelastic  zones  along  the  member  length,  the 
presence  of  residual  stresses,  and  strain  hardening  of  the  material  are 
also  accounted  for. 

Numerical  solutions  are  obtained  by  developing  a  set  of  finite 
element  equations  applicable  to  stability  problems.  However,  general 
purpose  computer  programs  have  been  developed  to  solve  a  variety  of 
problems.  The  most  advanced  of  these  can  handle  very  complex  nonlinear 
analyses,  but  does  not  require  an  excessive  amount  of  computer  time. 

1 . 4  Outline  of  Contents 

Chapter  2  reviews  stability  theory  and  discusses  different 
types  of  formulation  for  stability  and  strength  analyses. 

Chapter  3  is  divided  into  two  parts.  An  elastic  buckling 
formulation  is  presented  in  the  first  part,  with  a  solution  technique 
to  determine  the  elastic  critical  load  of  frames.  Numerical  results 
for  a  variety  of  problems  are  compared  with  those  obtained  through 
classical  techniques.  The  detailed  derivation  of  the  equations  for 
this  approach  is  carried  out  in  Appendix  A. 
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The  second  part  of  Chapter  3  presents  the  nonlinear  formula¬ 
tion,  based  on  large  deformation  theory.  The  equations  are  then 
specialized  to  elastic  response,  and  a  finite  element  model  is  developed 
for  this  case.  A  solution  method  to  solve  for  the  response  of  the 
frame  using  the  finite  element  equations  is  discussed  and  a  variety 
of  problems  are  solved  to  determine  the  accuracy  and  efficiency  of  the 
proposed  technique.  The  detailed  derivation  of  the  equations  for  this 
method  is  presented  in  Appendix  B,  with  a  large  displacement  trans¬ 
formation  matrix  derived  in  Appendix  C. 

Features  of  the  elastic-plastic  response  of  frames  for  first 
and  second  order  analyses  are  presented  in  Chapter  4.  Two  techniques 
are  discussed,  namely,  the  incremental  method  and  the  iterative  technique 
The  two  methods  and  their  results  are  compared  with  some  published  data. 

The  extension  of  the  approach  presented  in  part  2  of  Chapter  3, 
to  the  case  of  inelastic  response,  is  presented  in  Chapter  5.  The 
finite  element  equations  are  formulated  with  respect  to  local  reference 
axes  through  the  original  centroid.  Some  experimental  results  are 
compared  with  those  obtained  using  the  inelastic  formulation.  The 
detailed  derivation  of  the  section  properties  and  stress  resultants 
for  an  inelastic  element  is  given  in  Appendix  D,  and  details  of  the 
inelastic  formulation  are  shown  in  Appendix  E. 

Chapter  6  illustrates  a  behavioral  study  for  a  number  of 
frames  subjected  to  combined  vertical  and  lateral  loads.  Comparison 
is  made  between  the  inelastic  and  the  elastic-plastic  methods  of  analysis 
The  effects  of  residual  stresses  and  strain  hardening  on  the  behavior  of 
frames  are  also  studied. 

A  summary  of  the  investigation  and  the  conclusions  reached  are 


presented  in  Chapter  7. 


I 


' 


CHAPTER  II 


REVIEW  OF  STABILITY  THEORY 


2 . 1  Introduction 

Knowledge  of  structural  stability  theory  is  of  paramount 
importance  to  the  practicing  structural  engineer.  In  many  instances, 
buckling  is  a  primary  consideration  in  the  design  of  various  structural 
configurations.  Because  of  this,  and  because  the  formulation  for 
stability  problems  is  more  complex  than  for  strength  problems,  few 
other  fields  within  structural  engineering  have  such  a  varied  history 
as  the  theory  of  stability  of  metal  structures  (Bleich  1952) . 

This  chapter  presents  a  review  of  the  basic  concepts  of 
stability.  The  various  methods  of  analysis  used  in  solving  stability 
problems  are  discussed. 

2 . 2  Introduction  to  Linearly  Elastic  Formulation 

In  elastic  first  order  strength  analysis  two  main  assumptions 
are  considered.  The  material  of  the  structure  is  assumed  to  be  linearly 
elastic,  and  the  equilibrium  equations  are  formulated  on  the  undeformed 
shape  of  the  structure.  In  elastic  first  order  stability  analysis  the 
material  is  assumed  to  behave  linearly  elastic,  but  the  formulation  of 
the  equilibrium  equations  must  include  the  effect  of  displacements. 

2.2.1  Bifurcation  and  Buckling 

A  structure  is  considered  to  be  in  a  state  of  stable  equilibrium 
when  slight  changes  in  loading  do  not  produce  disproportionate 
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distortions  of  the  system.  A  load-displacement  history  such  as  the 
one  shown  in  Fig.  2-1  is  called  an  equilibrium  path.  Each  point  on  the 
path  represents  an  equilibrium  configuration  of  the  structure.  In 
linear  elastic  analysis  all  equilibrium  paths  are  straight  lines  that 
pass  through  the  origin  (Brush  and  Almroth  1975). 

Consider  a  straight,  simply  supported  column  of  length  L, 
subject  to  a  concentric  applied  compressive  load,  P,  as  illustrated  in 
Fig.  2-2a.  Under  the  assumptions  of  "small"  deflections,  v,  and 
constant  moment  of  inertia,  I  ,  the  equilibrium  equation  of  the  elastic 

X 

curve  is  given  by  Timoshenko  and  Gere  (1961)  as 

El  +  P  v  =  0  (2.2.1) 

x  dzz 

Introducing  the  notation 

k2  =  (2.2.2) 

X 

Equation  2.2.1  can  be  written  in  the  form 

v"  +  k2  v  =  0  (2.2.3) 

For  constant  k,  this  equation  is  a  linear  homogeneous  differ¬ 
ential  equation  the  general  solution  of  which  is  the  homogeneous 
solution  and  can  be  written  as 

v  =  A  sin  kz  +  B  cos  kz  (2.2.4) 

For  simply  supported  ends  the  boundary  condition  equations  are 


V  =  0 

at 

z  =  0 

(2.2.5a) 

v  =  0 

at 

z  =  L 

(2.2.5b) 

. 
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FIGURE  2-1  Load-Displacement  Relation  for  Beam 
Clamped  at  Both  Ends 
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(a)  Euler  Column 


(b)  Load-Deflection  Relationship 


FIGURE  2-2  Equilibrium  Paths  for  Initially  Straight  Column 
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Substituting  condition  2.2.5a  into  Eq.  2.2.4  leads  to  B  =  0. 

Consequently 

v  =  A  sin  kz  (2.2.6) 

From  condition  2.2.5b  one  obtains 

A  sin  kL  =  0 

Equation  2.2.7  can  be  satisfied  in  one  of  two  ways, 

A  =  0 

or 

sin  kL  =  0  (2.2.8b) 

If  A  =  0,  k  and  consequently  P  can  have  any  value.  But  if  A  =  0,  the 
displacement,  v,  of  Eq.  2.2.6  is  identically  zero.  This  result  is 
known  as  the  trivial  solution,  in  which  a  column  is  in  equilibrium  under 
any  value  of  axial  load,  P,  as  long  as  the  member  remains  perfectly 
straight . 

There  are  an  infinite  number  of  values  of  k  for  which 
sin  kL  =  0.  These  are 


(2.2.7) 

either 

(2.2.8a) 


k  L  =  n  TT 
n 


(2.2.9) 


Substituting  this  expression  into  Eq.  2.2.2  gives 
nontrivial  solutions  of  Eq.  2.2.1  can  be  obtained 


n2  TT2  El 


n 


the  loads  for  which 
as 


(2.2.10) 


Problems  of  this  type  are  called  eigenvalue  problems  (Ziegler  1968) . 

The  values  k  or  the  corresponding  loads  P  are  called  eigenvalues  of 
n  n 

the  problem,  and  the  corresponding  displacements 

v  =  A  sin  n  ^  Z  (2.2.11) 

n  n  L 

are  called  eigenfunctions. 

The  only  significant  value  of  P  arises  from  the  lowest  eigen¬ 
value  obtained  by  setting  n  equal  to  1  and  is  known  as  the  "Euler 
buckling  load": 

tt2  El 

Pe  =  (2.2.12) 

Equilibrium  paths  for  this  column  may  be  obtained  by  plotting 
P  versus  midspan  deflection  (v  at  z  =  ~)  .  Such  a  plot  for  n  =  1  is 
shown  in  Fig.  2-2b.  In  this  figure  the  primary  equilibrium  path 
(i.e.  the  trivial  solution)  is  intersected  by  a  secondary  equilibrium 
path.  The  point  at  which  equilibrium  paths  intersect  is  called  a  bi¬ 
furcation  point.  At  such  a  point  the  equilibrium  equations  have  multipl 
solutions,  one  corresponding  to  each  branch.  This  leads  to  the  classi¬ 
cal  concept  of  bifurcation  as  a  condition  in  which  two  (or  more) 
adjacent  equilibrium  configurations  exist  at  the  same  load. 

In  theory,  the  change  in  deformation  that  takes  place  when 
the  structure  passes  from  a  point  on  the  primary  equilibrium  path  to  a 
point  on  a  secondary  equilibrium  path  is  generally  known  as  buckling. 

For  a  linearly  elastic  stability  formulation,  buckling  occurs  at  the 
bifurcation  point.  In  general  the  load  at  which  this  happens  is  known 


as  the  critical  load. 
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Bifurcation  of  rectilinear  framed  structures  occurs  only  if 
the  members  remain  straight  during  loading.  Thus,  if  a  column  with  an 
initial  imperfection  vq(z)  is  considered,  as  illustrated  in  Fig.  2-3a, 
the  moment  at  any  section  along  the  column  is 


M  =  P  (v  +  v  ) 

o 


(2.2.13) 


and  hence  the  differential  equation  of  equilibrium  takes  the  form 


El  +  P  (v  +  v  )  =  0 

x  dz^  o 


(2.2.14) 


Dividing  by  El  and  using  the  notation  of  Eq.  2.2.2  leads  to 

X 


v"  +  k2v  =  -  k2v 

o 


(2.2.15) 


For  constant  k,  this  equation  is  a  linear  inhomogeneous 
differential  equation,  and  the  solution  consists  of  two  parts.  The 
first  is  the  homogeneous  solution,  identical  to  Eq.  2.2.4,  which  can  be 
obtained  by  setting  the  right  hand  side  in  Eq.  2.2.15  equal  to  zero. 

The  second  term  is  any  particular  solution  which  satisfies  Eq.  2.2.15. 

If  the  initial  shape  of  the  axis  of  the  bar  is  assumed  to  be 

v  (z)  =  a  sin  — ■  (2.2.16) 

o  L 

in  which  a  is  the  initial  midspan  deflection,  then  Eq.  2.2.15  becomes 

vM  +  k2v  =  -  k2  a  sin  (2.2.17) 

Li 

and  the  general  solution  of  this  equation  is  (Timoshenko  and  Gere  1961) 


■ 


Load,P 
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tp 

(a)  Imperfect  Column 


L 


FIGURE  2-3  Load-Deflection  Relation  for  Imperfect  Column 
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=  A  sin  kz  +  B  cos  kz  + 


7TZ 


(7T2/k2L2)  -  1 


a  sin  —  (2.2.18) 


To  satisfy  the  boundary  conditions  (v  =  0  at  z  =  0  and  z  =  L)  for  any 
value  of  k,  A  and  B  both  must  be  equal  to  zero.  If  the  ratio  between 
the  applied  load  and  the  critical  load  is  defined  as 


3 


P  =  P 

P  tt2EI/L2 

cr 


k2L2 

TT2 


(2.2.19) 


then  Eq.  2.2.18  becomes 

v  =  — — - — r-  a  sin  ~~  (2.2.20) 

1  —  p  L 

and  the  final  ordinates  of  the  deflected  curve  can  be  obtained  by  adding 
Eq.  2.2.16  to  Eq.  2.2.20,  which  yields 

V  =  sin  Y  (2.2.21) 

The  load-deflection  curve  for  this  imperfect  column,  as  P 
varies,  is  nonlinear  as  shown  in  Fig.  2-3b.  The  deflection  increases 
without  limit  as  P  approaches  P  . 

The  solution  to  the  homogeneous  equation  (Eq.  2.2.1)  is 
unique  and  identically  zero  until  the  load  reaches  the  critical  load 
which  corresponds  to  the  lowest  eigenvalue  of  the  differential  equation. 
At  this  point  the  magnitude  of  the  deflection  is  undefined,  because 
Eq.  2.2.7  is  satisfied  for  arbitrary  values  of  A.  The  solution  to  the 
inhomogeneous  equation  (Eq.  2.2.15)  is  unique  and  nontrivial,  until 
the  load  reaches  the  critical  load  of  the  homogeneous  problem,  at  which 
point  it  increases  without  limit.  Therefore,  for  this  type  of  problem. 


. 


18 


the  critical  load  of  the  homogeneous  solution  is  adequate  to  determine 
the  stability  limit  of  the  structure. 

It  should  be  noted  that  the  above  conclusion  is  valid  for  any 
arbitrary  initial  imperfection,  since  it  can  always  be  expanded  in  a 
Fourier  series  (Timoshenko  and  Gere  1961) .  The  conclusion  also  remains 
valid  for  a  column  with  transverse  loads,  since  these  simply  lead  to  a 
different  form  of  the  inhomogeneous  term  in  Eq.  2.2.15. 

2.2.2  Linearly  Elastic  Formulation  for  Frames 

The  study  of  frame  stability  is  associated  with  the  determina¬ 
tion  of  the  buckling  condition  of  a  system  of  interconnected  members 
as  well  as  with  the  determination  of  the  maximum  capacity  of  such  a 
system.  Frame  critical  loads  have  been  the  subject  of  numerous  investi¬ 
gations  (Lu  1962) ,  and  many  classical  techniques  are  available  for 
evaluating  such  critical  loads  (Bleich  1952).  The  following  discussion 
attempts  to  place  these  methods  in  a  general  context. 

The  most  fundamental  method  of  determining  critical  loads  for 
frames  is  a  direct  analytical  solution.  It  is  analogous  to  the  solu¬ 
tion  of  the  buckling  problem  for  a  pin-ended  column  (see  Sect.  2.2.1). 

Equation  2.2.1  is  valid  only  for  pin-ended  columns.  If  a 
column  is  fixed  at  both  ends  and  loaded  as  shown  in  Fig.  2-4,  the 
differential  equation  can  be  expressed  as 

EItt  =  -  M(z)  (2.2.22) 

dz 

in  which  M(z)  is  the  moment  at  any  point  on  the  column,  due  to  the 
end  effects  and  the  transverse  loads.  The  boundary  conditions  influence 
M(z)  in  an  indeterminate  manner.  However,  as  shown  by  Timoshenko  and 


. 


I- 


A  v 

P 


FIGURE  2-4  Fixed  End  Column 


20 


Gere  (1961) ,  a  single  fourth  order  equation,  applicable  to  any  prismatic 
column  regardless  of  the  boundary  conditions,  can  be  employed. 

Differentiating  Eq.  2.2.22  twice  with  respect  to  z  yields 


El 


d4v 

dz4 


+  P 


d2v 

dz2 


q(z) 


(2.2.23) 


Dividing  by  El,  and  using  Eq.  2.2.2,  leads  to 


d4v 

dz4 


+  k2 


dfv 

dz2 


(2.2.24) 


The  general  solution  of  Eq.  2.2.24  is 

v  =  C..  +  C„z  +  C0  sin  kz  +  C.  cos  kz  +  v  (z)  (2.2.25) 

1  Z  J  4  p 

in  which  the  integration  constants  C^,  C^,  and  are  determined 

from  the  boundary  conditions  of  the  case  under  investigation. 

In  most  structures  the  ends  of  the  columns  are  neither  hinged 
nor  fixed.  In  frames  the  columns  may  be  rigidly  connected  to  the  other 
members,  which  permits  a  limited  amount  of  rotation  to  occur  at  the 
ends  of  the  columns.  Supports  of  this  type  are  referred  to  as  elastic 
restraints  (Chajes  1974). 

Both  flexibility  solutions  and  stiffness  solutions  of  frame 
buckling  problems  can  be  obtained  directly  from  the  solution  of  the 
differential  equation  expressed  by  Eq.  2.2.25.  All  the  classical 
elastic  approaches  which  use  stability  functions,  become  special  cases 
of  this  method.  A  general  treatment  of  the  problem  may  be  carried  out 
as  follows. 

Equation  2.2.25  may  be  written  in  matrix  form  as 
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v(z)  =  <  4) ( z)  >  {c}  +  v  (z) 

P 


in  which 


(2.2.26a) 


<  <P  >  =  <  1,  z,  sin  kz ,  cos  kz  > 


(2.2.26b) 


and 

{C}  =  <  cv  c2,  c3,  c4  >T 


(2.2.26c) 


Using  Eqs.  2.2.26  to  evaluate  the  transverse  displacements  and  rotations 
at  the  ends  of  the  member,  and  adopting  the  notation  of  Fig.  2-5a, 
yields 


/  \ 

f  \ 

5(0) 

<  0(0)  > 

vp(0) 

6(0) 

A 

o 

n — * 

V 

v’(0) 

< 

>  = 

{c}  +  ^ 

P  r 

5  (L) 

A 

/— s 

►J 

V 

Vp(L) 

0(L) 

A 

to 

-©■ 

V 

v’(L) 

V. 

P 

which  may  be  written,  symbolically,  as 


(2.2.27a) 


{A}  =  [G1]  (C>  +  {Ap}  (2.2.27b) 

Similarly,  the  end  forces  may  be  evaluated  from  Eqs.  2.2.26.  Using  the 
fundamental  relations  for  moment  and  shear  (Timoshenko  and  Gere  1961)  , 
the  following  expressions  can  be  written 


and 


M  =  -  El  v" 


(2.2.28a) 


V  =  -  El  vMf  -  P  v'  (2.2.28b) 

Adopting  the  notation  of  Fig.  2-5a,  Eqs.  2.2.28  yield 


f 

v(0) 

~ei  < 0"'(O)>  -  P  < 0’  (0)  >" 

El  v’"  (0)  -  P  v’  (0) 

P  P 

M(0) 

-El  <  0"(O)  > 

-  El  v"  (0) 

< 

►  = 

{c} 

p  r 

V  (L) 

-El  <0,M  (L)>-P<0*  (L)> 

-  El  v"  ’(L)  -  P  v*  (L) 

P  P 

M(L) 

N  J 

El  <  0"(L)  > 

El  v"(L) 

(2.2.29a) 
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FIGURE  2-5  Behavior  of  Frames 
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Using  the  notation  of  Eq.  2.2.2  leads  to 


t  > 

r 

v(0) 

<4>”’  (0)  >  -  k2  <  <£’  (0)  > 

vMf  (0)  -  k2  v’  (0) 

P  P 

M(0) 

A 

/■ — s 

o 

V 

i 

-  v"  (0) 

< 

*  =  EX 

{c}  +  EI  < 

P  r 

V  (L) 

-  <  4)'"  (L)  >  -  k2  <  4)1  (L)  > 

-  v™  (L)  -  k2vr  (L) 

P  P 

M(L) 

+  <  4>m(l)  > 

v"  (L) 

s.  y 

— 

l  P  J 

(2.2.29b) 


which  may  be  written,  symbolically,  as 

{M}  =  [Gj  {C}  +  (M  }  (2.2.29c) 

2  P 

Equations  2.2.27  and  2.2.29  become  the  basis  for  the  following 
developments . 


2.2.2a  Flexibility  Approach 

Solving  Eq.  2.2.29c,  for  {c},  yields 

(C>  =  [ G2 ]_1  {M -  Mp}  (2.2.30) 

Substituting  this  result  into  Eq.  2.2.27b  yields 

{A}  =  [G1] [G2]-1{M-  Mp}  +  {Ap}  (2.2.31a) 

or 

{A}  =  [F]  (M-M  }  +  (A  }  (2.2.31b) 

P  P 

The  matrix  [F]  in  Eq.  2.2.31b  is  the  (4x4)  member  flexibility  matrix 
for  the  displacement  coordinates  and  their  associated  forces  (Murray 
et  al.  1977) ,  as  indicated  in  Fig.  2-5a. 

A  classical  flexibility  analysis  can  now  be  carried  out  by 
arranging  all  the  element  vectors  appearing  in  Eq.  2.2.31b  sequentially 
into  global  vectors,  and  writing  the  result  as 
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{r}  =  [F]  {R*}  +  {r  }  (2.2.32) 

P 

For  the  frame  illustrated  in  Fig.  2-5b  and  2-5c  the  vector  {r}  would 
appear  as 


<r>  =  <  6^,  0^,  6^,  0^ ,  6^,  0^,  6^,  6^*  6^,  0^,  02  > 

(2.2.33) 

The  vector  {r^}  is  a  similar  array  of  particular  solution  displacements, 
whereas  the  vector  {R*}  is  the  associated  array  of  end  forces  and 
moments,  in  which  each  component  represents  the  difference  between 
the  total  value  and  that  associated  with  the  particular  solution.  The 
matrix  fFj  contains  the  element  flexibility  matrices  of  Eq.  2.2.31b 
on  its  principal  diagonal. 

Compatibility  equations  for  the  frame  illustrated  in  Fig.  2-5b, 

are 
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1 

1 


5 


1 

2 


+ 


0 


6 


1 

2 


0 


e 


2 

2 
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Equations  2.2.34  and  2.2.35  may  be  written  as 


(2.2.34) 

(2.2.35a) 

(2.2.35b) 

(2.2.35c) 


[A]  {r} 


{0} 


(2.2.36) 
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in  which  all  elements  of  the  compatibility  matrix  [A]  are  either 
0,  1  or  -1.  Substituting  Eq.  2.2.32  into  Eq.  2.2.36  yields 

[A]  fFJ  {R*}  =  -  [A]{r  }  (2.2.37) 

This  compatibility  equation  cannot  yet  be  solved  for  the  indeterminate 
forces  {R*},  since  there  are  only  9  equations  but  there  are  12  unknowns. 

To  complete  the  solution,  the  left-hand  side  of  Eq.  2.2.36 
may  be  considered  to  define  a  set  of  relative  displacements  { q } ,  such 
that 


{q}  =  [A]  {r}  (2.2.38) 

A  set  of  generalized  forces  {Q*}  may  be  associated  with  these  relative 
displacements  and,  for  work  equivalence;  it  is  necessary  that 
(Przemieniecki  1968) 

{R*}  =  [A]T  {Q*}  (2.2.39) 

Substituting  Eq.  2.2.39  into  Eq.  2.2.37  yields 

[A]  m  [A]T  {Q*}  =  -  [A]  {rp}  (2.2.40a) 

which  may  be  written  symbolically  as 

[F]  {Q*}  =  -  {qp}  (2.2.40b) 

Equation  2.2.40b  is  a  set  of  compatibility  equations  which 
permit  the  solution  for  the  generalized  redundant  forces  {Q*},  which 
restore  continuity  to  the  frame  by  counteracting  the  incompatibility 


■  J|  L  i 
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displacements  { } .  The  matrix  [F]  is  the  flexibility  matrix  associated 
with  these  relative  displacements. 

For  the  condition  where  distributed  loads  are  applied  to  the 
frame,  the  particular  solution  v^(z)  of  Eq.  2.2.25  is  non-zero. 

Hence,  {qp}  of  Eq.  2.2.40b  is  non-zero.  For  any  value  of  distributed 
load  and  axial  forces,  Eq.  2.2.40b  can  be  solved  for  {Q*};  {R*}  can 
be  determined  from  Eq.  2.2.39;  {c}  from  Eq.  2.2.30,  and  hence  the 
solution  of  Eq.  2.2.26a  is  known  for  each  member.  However,  as 

|  |  F  |  |  ->  0  (2.2.41) 

in  which  ||  F  ||  indicates  the  determinant  of  [F] ,  the  redundants 
{Q*},  and  hence  all  displacements,  increase  without  limit.  Similarly, 
if  all  Vp(z)  are  identically  equal  to  zero,  only  the  trivial  solution 
is  possible  until 

| |  F  | |  =  0  (2.2.42) 

at  which  time  the  eigenvectors  of  the  homogeneous  form  of  Eq.  2.2.40b 
can  be  determined,  and  forces  and  the  associated  displacements  of 
arbitrary  magnitude  may  exist.  Thus,  a  bifurcation  solution  of  the 
homogeneous  problem  and  unbounded  solutions  of  the  inhomogeneous  problem 
arise  in  the  same  manner  as  for  a  simple  column,  and  permit  the  same 
interpretation,  as  illustrated  in  Fig.  2-5d. 

The  above  technique  is  a  general  approach  and  has  been  used 
for  simple  structures.  Normally  an  analyst  imposes  the  boundary 
conditions  of  Eq.  2.2.34  on  the  differential  equations  prior  to 
writing  the  intermember  compatibility  conditions  of  Eq.  2.2.35 
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(Timoshenko  and  Gere  1961) .  This  reduces  the  size  of  the  flexibility 
matrix  in  the  illustrative  problem  from  9  x  9  to  3  x  3.  In  non¬ 
computerized  analysis  this  could  reduce  the  problem  from  an  intract¬ 
able  solution  to  a  tractable  one.  However,  determining  critical  loads 
by  the  flexibility  method  is  still  formidable.  It  is  necessary  to 
assume  or  calculate  the  value  of  the  axial  force  in  each  member,  in 
order  to  determine  the  values  of  k  in  the  general  solution,  Eq.  2.2.25. 

It  may  then  be  assumed  that  these  axial  forces,  and  hence  the  k- 
values  of  the  members,  are  proportional  to  the  applied  loads.  If 
the  load  factor  from  the  reference  load  is  denoted  by  A,  the  stability 
condition  (Eq.  2.2.42)  becomes  a  transcendental  function  of  A.  It  is 
necessary  to  find  the  smallest  root,  A  ,  of  this  'characteristic' 
equation . 

The  most  practical  method  of  doing  this  is  by  trial  and  error, 
and  this  is  called  the  'determinant  search'  technique.  Values  of  A 
are  assumed,  and  ||  F  ||  is  evaluated  for  each  assumed  value.  By  plotting, 
or  interpolating  between  pairs  of  values,  the  value  of  A  can  be 
determined  (Bleich  1952) . 

2.2.2b  Stiffness  Approach 

A  classical  stiffness  formulation,  exactly  parallel  to  the 
flexibility  formulation,  may  be  derived  as  follows.  For  a  stiffness 
formulation  Eq.  2.2.27b,  instead  of  Eq.  2.2.29c,  is  solved  for  (c), 
to  yield 

{C}  =  [GjT1  {A  -  Ap} 


(2.2.43) 
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Substitution  into  Eq.  2.2.29c,  yields 

{M}  =  [Gj  [G.]”1  {A  -  A  }  +  {M  }  (2.2.44a) 

2  1  p  p 

which  may  be  written  symbolically  as 

{M}  =  [k]  (A  -  A  }  +  (M  }  (2.2.44b) 

P  P 

in  which  the  matrix  [k]  is  the  (4  x  4)  member  stiffness  matrix, 
expressing  the  member  end  forces  of  Fig.  2-5a  in  terms  of  the  corres¬ 
ponding  displacements. 

Arranging  the  element  vectors  appearing  in  Eq.  2.2.44b 
sequentially  into  global  vectors,  Eq.  2.2.44b  for  all  members  may 
be  combined  into  the  single  matrix  equation 

{R}  =  fKj  {r*}  +  {Rp}  (2.2.45) 

For  the  frame  illustrated  in  Figs.  2-5b  and  2-5c  the  vector  {r}  would  be 

<R>  =  <  vj,  M^,  vj;,  m£,  vj,  M^,  v2z,  m^,  v^,  mJ,  v^,  > 

(2.2.46) 

The  vector  (R  }  is  a  similar  array  of  particular  solution  end  forces, 

P 

and  the  vector  {r*}  is  the  corresponding  array  of  end  displacements, 
where  each  term  represents  the  difference  between  the  total  value  and 
the  particular  solution  value  of  the  variable.  The  matrix  fkj 
contains  the  element  stiffness  matrices  of  Eq.  2.2.44b  on  its  princi¬ 
pal  diagonal. 

Equilibrium  equations  may  now  be  written  which,  for  the  frame 


of  Fig.  2-5b,  are 
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M1  +  M2  +  V2  *  L1  =  0  (2.2.47a) 

+  V22  •  L2  =  0  (2.2.47b) 

•  L  =  0  (2.2.47c) 

V2  +  V1  =  °>  V1  +  V2  =  °»  V1  +  V2  =  °  (2.2.47d) 

=  0,  =  0,  =  0  (2.2.47e) 

Equations  2.2.47  may  be  expressed  symbolically  as 

[A]  (R)  =  {0}  (2.2.48) 


in  which  the  elements  of  the  equilibrium  matrix  [A]  are  either  0,  1, 
or  -1. 

Substituting  Eq.  2.2.45  into  Eq.  2.2.48  yields 

[A]  fKj  {r*}  =  -[A]  {Rp}  (2.2.49) 

Since  Eq.  2.2.49  contains  only  9  equations,  and  there  are  12  unknowns 
in  {r*},  it  cannot  be  solved  directly.  However,  Eq.  2.2.48  can  be 
considered  to  define  a  set  of  unbalanced  forces  {Q},  such  that 

{Q}  =  [A]  {R}  (2.2.50) 

for  which  the  associated  work  equivalent  set  of  generalized  relative 
displacements  {q*}  are  related  to  the  displacements  {r*}  by  the 
equation 

{r*}  =  [A]T  {q*}  (2.2.51) 

Substituting  Eq.  2.2.51  into  Eq.  2.2.49  yields 
[A]  fKj  [A]T  {q*}  =  -  [A]  {R  } 


(2.2.52a) 
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which  may  be  written  symbolically  as 

[K]  {q*}  =  -  {Q  }  (2.2.52b) 

Equations  2.2.52  are  a  set  of  equilibrium  equations  that 
permit  the  evaluation  of  the  generalized  displacements  {q*}  which 
establish  equilibrium  between  members  of  the  frame.  The  matrix  [K] 
is  the  stiffness  matrix  associated  with  the  generalized  displacements. 

For  any  distributed  loading  and  axial  forces,  the  particular  solution 
forces  and  the  vector  {Q  }  may  be  formed.  Similarly,  the  matrix  [K] 
may  be  constructed.  Once  {q*}  has  been  determined  from  Eq.  2.2.52b, 

{r*}  may  be  determined  from  Eq.  2.2.51  and  all  member  end  forces  from 
Eq.  2.2.45. 

The  stability  condition  is  determined  in  the  same  manner 

\ 

as  was  done  for  the  flexibility  formulation  of  Sect.  2.2.2a.  A  unique 
set  of  finite  displacements  arises  from  the  solution  of  the  equilibrium 
equations  (Eq.  2.2.52b).  However,  as  axial  forces  increase  and 

||  K  ||  -*■  0  (2.2.53) 

displacements  increase  without  limit.  For  the  homogeneous  case,  the 
solution  to  Eq.  2.2.52b  is  the  trivial  solution  unless 

||  K  ||  =0  (2.2.54) 

in  which  case  nontrivial  eigenvectors  of  arbitrary  magnitude  exist 
as  solutions.  The  interpretation  of  these  solutions  is,  again,  as 
illustrated  in  Fig.  2-5d.  The  determination  of  the  critical  load 
can  be  accomplished  by  a  determinant  search  technique  as  discussed  in 


Sect.  2.2.2a. 
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The  above  stiffness  technique  represents  a  general  approach. 

It  has  been  widely  used  in  the  more  conventional  form  of  the  slope- 
deflection  technique  or  the  moment  distribution  procedure.  The  element 
slope-deflection  matrix  arises  directly  from  Eq.  2.2.44b,  as  shown  in 
the  following: 

Equation  2.2.44b  can  be  written  as 

{M}  =  [k]  {A}  +  (M  }  -  [k]  {A  }  (2.2.55) 

P  P 

For  zero  member  displacements  the  total  forces,  called  the  fixed  end 
forces  and  denoted  by  {MF} ,  are  obtained  from  Eq.  2.2.55  as 

{MF}  =  (M  }  -  [k]  (A  }  (2.2.56) 

P  P 

Consequently  Eq.  2.2.55  may  be  written  as 

{M}  =  [k]  {A}  +  (MF)  (2.2.57) 


In  the  slope-deflection  method,  only  the  end  moments  are  selected  from 
the  vector  (m) ,  to  yield 
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(2.2.58) 


in  which  <k> .  is  the  ith  row  of  the  matrix  [k] .  Recognizing  that  the 

l  — 

coefficients  of  6^  and  in  these  rows  are  numerically  the  same  but 
opposite  in  sign  permits  Eq.  2.2.58  to  be  written  as 
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(2.2.59) 


in  which  C  and  S  are  the  standard  stability  functions  (Galambos  1968) . 

Assembly  of  Eq.  2.2.59  into  the  equilibrium  equations  Eq.  2.2.47 
by  a  standard  slope  deflection  technique  gives  Eq.  2.2.52b,  and  the 
solution  procedure  then  follows  as  described  above.  Thus,  the  magni¬ 
tude  of  the  axial  forces  which  satisfy  Eq.  2.2.54  are  sought  by  a 
determinant  search  technique. 

Equation  2.2.59  is  also  the  starting  point  for  a  moment  distri¬ 
bution  search  for  critical  loads  (Hoff  1941) .  In  this  technique  the 
magnitudes  of  the  axial  forces  are  assumed  and  the  stiffnesses,  S, 
and  carry-over  factors,  C,  are  evaluated.  Distribution  factors  are 
computed  with  these.  At  the  critical  load  the  moment  distribution 
process  diverges,  leading  to  end  forces  and  displacements  which  become 
arbitrarily  large  (Winter  et  at.  1948).  This  is,  in  effect,  a  test 
that  Eq.  2.2.54  has  been  satisfied. 


2.3  Elastic  Nonlinear  Formulation 


An  elastic  nonlinear  formulation  may  be  carried  out  to  deter¬ 
mine  the  behavior  of  the  structure  after  the  deflections  are  large 
enough  to  cause  significant  changes  in  the  geometry.  This  phenomenon 
is  referred  to  as  geometric  nonlinearity.  To  account  for  the  effects 
of  changes  in  the  geometry  as  the  applied  loading  is  increased,  an 
exact  expression  for  the  curvature  must  be  used. 
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In  Sect.  2.2.1  the  deflection  of  a  column  was  found  to  be 

indeterminate  at  the  critical  load.  This  is  strictly  correct  only  as 

long  as  the  deflection  remains  "small"  and  may  be  better  understood 

by  recognizing  that  the  governing  equation,  Eq.  2.2.1,  is  based  on  the 

d2v 

approximate  expression  for  the  curvature  of  the  buckled  column  . 

If  the  exact  expression  is  used,  the  solution  will  not  exhibit 
indefiniteness  in  the  value  of  the  deflection.  The  shape  of  the 
elastic  curve  of  a  simple  column,  when  developed  on  the  basis  of  the 
"exact"  differential  equation,  is  called  the  'elastica'  . 

2.3.1  Large-Deformation  Formulation  For  Columns  (The  Elastica) 

To  obtain  a  solution  for  the  ’elastica’,  consider  the  simply 
supported  column  shown  in  Fig.  2-6a.  Except  for  the  now  invalid 
assumption  of  "small"  deflections,  all  idealizations  are  the  same  as 
those  of  the  Euler  column.  If  the  y-z  coordinate  system  is  taken  as 
shown  in  Fig.  2-6a,  and  the  column  is  in  equilibrium  in  a  bent  con¬ 
figuration,  then  the  moment  (Pv)  at  any  section  is  equal  to  the  resisting 
moment  (-EI$) .  Thus, 

-  EI<j>  =  Pv  (2.3.1) 


in  which  <j)  is  the  angle  change  per  unit  length.  If  the  slope  is  given 
by  9,  and  the  distance  along  the  curve  from  the  origin  to  the  point 
is  s,  then 


d_9 

ds 


(2.3.2) 


Substituting  this  expression  into  Eq.  2.3.1  leads  to 


d9 

El  -j—  +  Pv  = 
ds 


0 


(2.3.3a) 
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(a)  Large  Deflections  of  a  Column 


FIGURE  2-6  Column  Behavior  for  Large  Deflections 
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or 


+  k2v  =  0 

ds 


(2.3.3b) 


dv 

Differentiating  Eq.  2.2.3b  with  respect  to  s  and  replacing  —  by 

Q  S 

sinQ,  yields 


d29 

ds2 


+  k2  sin9  =  0 


(2.3.4) 


which  is  a  nonlinear  differential  equation  in  9. 

Multiplying  Eq.  2.3.4  by  (2  d9)  and  integrating  gives 


(-j^-) 2  -2  k2  cos9  =  C 
ds 


(2.3.5) 


The  integration  constant,  C,  is  evaluated  from  the  boundary  conditions 
of  Fig.  2-6a,  namely. 


J  A 

9  =  a  and  —  =0  at  z  =  0 

ds 


(2.3.6) 


Hence 


C  =  -2k2  cosa 


(2.3.7) 


and  Eq.  2.3.5  becomes 


€>• 


=  2k2  (cos9  -  cosa) 


(2.3.8) 


or 


d9 

ds 


±  k  /  2  /  cos9  -  cosa 


(2.3.9) 
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Solving  for  ds  gives 


ds  =  - 


d9 


k  nr  r  cos9  -  cosa 


in  which  the  positive  sign  has  been  dropped  because  9  always 


as  s  increases. 


The  total  length  of  the  column  is  obtained  from 


-CL 


L  = 


ds 


d9 


+a 


k  nr  r  cos9  -  cosa 


Equation  2.3.11  may  be  written  as  (Timoshenko  and  Gere  1961) 


L  = 


2K 


/  P/EI 


in  which 


K 


f+a 

d0 


,  /  .  o  a  .20 

;-a  4  /  smz  ~  sm  — 


is  a  complete  elliptic  integral  of  the  first  kind,  which  is 
tabulated  in  mathematical  handbooks  (Jahnke  and  Emde  1945) . 
Solving  Eq .  2.3.12  for  P  yields 

P  =  4K2 

P  IT2 

cr 


(2.3.10) 

decreases 


(2.3.11) 


(2.3.12) 


(2.3.13) 


(2.3.14) 
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in  which  P  =  tt2EI/L2 
cr 

If  the  deflection  of  the  member  is  very  small  the  value  of 
TT 

K  approaches  -j.  Equation  2.3.14  then  becomes 


P 


P  = 
cr  L2 


(2.3.15) 


The  nonlinear  theory  thus  leads  to  the  same  critical  load  as  the  linear 
theory. 

The  load-deformation  characteristic  of  the  column  can  be 
obtained  as  follows. 

Noting  that  dv  =  sinQ  ds,  and  making  use  of  Eq.  2.3.10,  gives 


,  sinQ  d9 

dv  =  -  - 

k  /  2  /  cos0  -  cosa 


(2.3.16) 


The  midheight  deflection,  5,  is  obtained  as 


5 


o 


dv 

■  6 


a 

_ sin9  d9 _ 

k  Z~2~  J~  cos9  -  cosa 
o 


(2.3.17) 


from  which  it  can  be  shown  that  (Timoshenko  and  Gere  1961) 


5_  _  2  sin  a/2 

L  tt  /  P/P 

cr 


(2.3.18) 


Using  Eqs.  2.3.14  and  2.3.18  it  is  possible  to  compute  for 
various  values  of  a  the  corresponding  values  of  P/P  and  <$/L.  The 
load-deformation  characteristic  is  shown  in  Fig.  2-6b.  Comparison 


« 


between  this  figure  and  Fig.  2-2b  shows  that  large  deformation  theory 
leads  to  the  same  critical  load  as  linear  theory,  but  predicts  a  very 
slight  increase  in  load  with  increasing  deflection  above  the  critical 
load.  The  increase  in  load  becomes  pronounced  only  after  consider¬ 
able  deformation  has  taken  place.  This  indicates  that  the  prediction 
of  bending  at  constant  load,  as  implied  by  the  linear  theory,  is  a 
good  approximation  for  a  considerable  range  of  deformations. 

2.3.2  Elastic  Nonlinear  Analysis  of  Frames 

In  a  structural  framework  made  from  a  linearly  elastic 
material,  nonlinearity  arises  when  the  effect  of  the  deformations  on 
the  equilibrium  equations  is  considered. 

In  Sect.  2.2.2  an  approximate  expression  for  the  curvature 
was  used,  and  the  axial  force  in  each  column  was  assumed  to  remain 
constant  during  buckling.  These  assumptions  result  in  a  linear 
formulation  of  the  problem.  The  classical  solutions,  as  presented 
in  Sect.  2.2.2,  may  be  modified  to  account  for  geometric  nonlinearity. 
In  this  case  an  exact  expression  for  the  curvature  is  used,  resulting 
in  an  ’exact'  differential  equation  (See  Sect.  2.3.1).  The  axial 
force  also  may  be  assumed  to  be  changing  in  each  column  during 
buckling.  A  modified  classical  formulation,  however,  would  not  be 
easy  to  solve.  A  few  authors  have  considered  the  nonlinear  strain- 
displacement  relationship  (Turner  et  at.  1960;  Martin  1965),  and 
many  others  have  incorporated  the  effect  of  joint  displacement  on  the 
overall  behavior  of  the  structure  (Nair  1975) . 

An  approximate  method  of  analysis  for  elastic  nonlinear 
problems  has  been  described  by  Adams  (1972).  The  effect  of  joint 
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displacements  on  the  overall  behavior  of  a  multistory  frame  is 
accounted  for  by  modifying  the  results  of  a  standard  (linear  formula¬ 
tion)  analysis.  The  method  is  usually  referred  to  as  the  P-A  method, 
and  the  secondary  effects  as  the  P-A  effects. 

The  mechanics  of  the  P-A  method  can  be  described  as  follows, 
using  Fig.  2-7  (Johnston  1976). 

1  -  The  loads  are  applied  to  the  frame,  and  the  lateral  displacements 

A_^  due  to  primary  moments  are  computed. 

2  -  The  additional  story  shears  that  are  caused  by  the  displaced 

vertical  loads  are  calculated  as 


3 


V! 

l 


A.) 


(2.3.19) 


i 

in  which  is  the  additional  shear  in  story  i  due  to  the  sway 

forces,  ZP .  is  the  sum  of  the  column  axial  loads  in  story  i,  h. 

l.  l 

is  the  height  of  story  i,  and  A^+^  and  A^  are  the  lateral  dis¬ 
placements  of  the  frame  at  levels  i+1  and  i,  respectively. 

~  The  fictitious  sway  forces,  H^,  are  computed  at  each  floor  level  as 


H!  =  V!  -  -  V!  (2.3.20) 

1  1-1  1 


4  ~  The  sway  forces,  H^,  are  added  to  the  applied  loads  and  a  new 

linear  analysis  is  performed. 

5  ~  When  the  deflected  shape  remains  unchanged  from  one  iteration  cycle 

to  the  next,  the  solution  has  converged,  and  the  resulting  forces 
now  include  the  secondary  moments. 

In  practical  applications  the  convergence  is  fast,  and  the 
first  iteration  produces  acceptable  results  (Springfield  and  Adams 
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FIGURE  2-7  Sway  Forces  Due  to  Vertical  Loads 
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1972) .  Slow  convergence  is  a  sign  that  the  structure  is  excessively 
flexible,  and  lack  of  convergence  indicates  that  the  structure  is 
unstable  (Wood  et  at.  1976) . 


2 . 4  Inelastic  First  Order  Stability  Theory 

Early  tests  of  steel  columns  of  practical  proportions  showed 
that  failure  occurred  at  loads  less  than  the  Euler  load.  The  difference 
can  be  attributed  to  the  assumptions  made  in  deriving  the  Euler  load. 

In  each  of  the  preceding  sections  it  has  been  assumed  that  the 
material  behaves  linearly  elastic,  but  in  order  for  this  to  be  valid, 
the  total  stresses  in  the  member  must  remain  below  the  proportional 
limit  of  the  material.  However,  in  most  columns  the  elastic  limit  is 
exceeded  before  the  load  reaches  the  Euler  load.  The  results  of  the 
elastic  analysis  therefore  are  not  valid  and  the  buckling  load  must 
be  determined  by  taking  inelastic  behavior  into  account. 


2.4.1  Tangent  Modulus  Theory 

Engesser  (1889)  suggested  that  if  column  failure  occurred 
at  a  stress  above  the  proportional  limit  of  the  material,  the  column 
strength  could  be  obtained  by  simply  replacing  E  by  Et»  the  tangent 
modulus,  in  the  Euler  formula  (Eq.  2.2.12).  The  slope  of  the  tangent 
to  the  stress-strain  curve  of  the  material  at  any  point  A  is  E^ ,  as 
shown  in  Fig.  2-8.  That  is 


E 


t 


da 

de 


(2.4.1) 


To  justify  the  substitution  of  E^  instead  of  E,  it  is 
assumed  that  the  column  remains  straight  and  that  the  stress-strain 


Load,  P 
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FIGURE  2-9  Load-Deflection  Curve  of  Shanle/s  Model 


characteristics  are  the  same  throughout  the  length  of  the  member. 
Thus,  the  average  stress  at  the  tangent  modulus  load  for  a  short 
column,  such  as  the  one  shown  in  Fig.  2-2a,  is 


o 

cr 


(L/r)2 


(2.4.2) 


2.4.2  Reduced  Modulus  Theory 

When  a  column  begins  to  bend  at  the  critical  load,  there  is 
a  possibility  that  stresses  on  the  convex  side  decrease.  Referring 
to  Fig.  2-8,  any  incremental  decrease  in  the  compressive  strain  from 
point  A  involves  a  decrease  in  stress  equal  to  (Ae)E,  whereas  an 
incremental  increase  in  the  strain  involves  an  increase  in  stress 
equal  to  (AeOE^..  This  line  of  reasoning  is  the  basis  for  the  reduced 
modulus  theory,  or  as  it  is  sometimes  called,  the  double  modulus 
theory. 

Based  on  the  above  reasoning,  Engesser  (1895)  realized  the 
theoretical  contradictions  inherent  in  the  tangent  modulus  theory,  and 
suggested  the  use  of  a  reduced  modulus,  E  ,  for  E  in  the  Euler  formula. 
The  magnitude  of  E^  lies  between  those  of  E  and  E  .  The  expression 
for  E^  and  the  position  of  the  axis  of  bending  at  critical  load  are 
determined  so  as  to  satisfy  the  traditional  buckling  concept  that 
the  load  remains  constant  during  buckling. 

Von  Karman  (1910)  revised  the  theory,  and  determined  the 
expressions  for  the  reduced  modulus  for  a  rectangular  cross  section 
and  for  the  idealized  I-section.  His  subsequent  experimental  work 
showed  that  the  actual  buckling  loads  were  closer  to  the  tangent 
modulus  values  than  to  the  reduced  modulus  values.  In  contradiction 


to  the  understanding  of  column  behavior  at  that  time,  this  problem  was 
not  resolved  for  another  35  years. 

2.4.3  Shanley’ s  Contribution 

Shanley  (1947)  resolved  the  apparent  contradiction  between 
theory  and  tests  through  a  revised  tangent  modulus  model.  He 
demonstrated  that  upon  reaching  the  tangent  modulus  load,  P  ,  there 
is  nothing  to  prevent  the  column  from  bending  and  at  the  same  time 
support  an  increasing  axial  load. 

Figure  2-9  illustrates  the  load-deflection  curve  according 
to  Shanley’ s  model.  He  concluded  that  the  column  load  may  exceed 
the  tangent  modulus  load,  but  cannot  be  greater  than  the  reduced 
modulus  load,  P^_.  Shanley  also  noted  that  the  tangent  modulus  theory 
predicts  more  accurately  the  maximum  load  that  an  inelastic  column 
may  support. 

2.4.4  Inelastic  Behavior  of  Beam-Columns 

Inelastic  action  must  be  considered  in  the  determination  of 
the  ultimate  strength  of  a  beam-column.  The  behavior  of  beam-columns 
is  different  from  the  behavior  of  either  columns  or  beams.  The  term 
defines  a  member  subjected  to  axial  load  as  well  as  bending  moment. 

A  set  of  tests  was  performed  by  Ketter  et  at.  (1955)  on  W8  x  31  beam- 
columns  loaded  up  to  failure.  The  results  showed  that  the  fully 
plastic  moment  capacity  of  a  beam-column  cannot  be  reached.  The  study 
concluded  that  the  maximum  moment  capacity  of  members  is  a  function 
of  the  axial  force  as  well  as  the  moment  applied  to  the  member. 

Galambos  and  Ketter  (1961)  presented  theoretical  ultimate 
strength  interaction  curves  for  beam-columns  with  various  end 


conditions,  correlating  the  analysis  with  earlier  tests  by  Mason 
et  at.  (1958).  One  set  of  such  curves  is  shown  in  Fig.  2-10, 
representing  strong  axis  bending  of  a  W8  x  31,  subjected  to  equal 
end  moments.  The  curves  are  a  reasonable  approximation  for  most 
typical  wide-flange  sections.  The  validity  of  such  interaction 
curves  has  been  verified  by  extensive  experiments  (Van  Kuren  and 
Galambos  1964),  and  the  correlation  between  tests  and  theory  is  good. 

The  strength  of  members  subjected  to  bending  moments  and 

axial  loads  can  be  approximated  by  interaction  formulas  in  terms  of 

the  ratios  P/P  and  M  /M  as 

u  o  u 


D  CM 

— - 1 - m.  ■  .° -  <  IQ 

P  M  (1  -  P/P  )  -  i,U 

u  u  e 


(2.4.3) 


in  which  P  is  the  axial  load  at  failure,  P  is  the  ultimate  load  for 

u 

a  concentrically  loaded  pin-ended  column,  Mq  is  the  maximum  applied 
moment  (not  including  second  order  effects),  and  M^  is  the  ultimate 
moment  capacity  in  the  absence  of  axial  load.  The  coefficient 
is  the  equivalent  moment  factor  that  accounts  for  unequal  end 
moments,  and  Pg  is  the  elastic  critical  load  for  buckling  in  the 
plane  of  the  applied  moments.  The  term  1/(1  -  P/Pg)  in  Eq.  2.4.3  is 
the  amplification  factor,  covering  the  effect  of  the  axial  load 
acting  on  the  deflected  column. 

Equation  2.4.3  is  based  on  failure  of  the  beam-column  due  to 
instability.  It  approximates  the  experimental  interaction  curves, 
and  is  the  basis  for  the  Canadian  Standards  Association  Specification 
(CSA  1974). 
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FIGURE  2-10  Maximum  Carrying  Capacity  Interaction  Curves 
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2.4.5  Simple  Plastic  Analysis 

The  objective  of  the  plastic  methods  of  analysis  is  to 
predict  the  loads  at  which  a  structure  may  fail  by  the  development  of 
excessive  deflections.  It  is  assumed  that  whenever  the  moment  at  any 
section  reaches  the  plastic  moment  capacity,  Mp,  a  plastic  hinge 
forms,  which  can  undergo  extensive  rotation  while  the  moment  remains 
constant.  The  arrangement  of  plastic  hinges  together  with  real 
hinges  which  are  functioning  at  collapse  is  known  as  the  collapse 
mechanism  or,  simply,  mechanism  (Neal  1977).  The  plastic  limit  load 
for  a  frame  is  the  maximum  load  it  can  carry,  assuming  that  no 
increases  will  take  place  after  a  mechanism  has  formed. 

Two  approaches  were  commonly  used  in  plastic  analysis  in 
early  studies  (Pippard  and  Baker  1968)  namely,  the  equilibrium  method 
and  the  mechanism  method. 

The  equilibrium  method  of  analysis  consists  of  constructing  a 
moment  diagram  for  the  structure  after  it  has  been  rendered  statically 
determinate  by  replacing  a  sufficient  number  of  the  redundant  moments 
with  the  fully  plastic  moment  Mp  (Phillips  1956) .  Since  it  is  not 
always  possible  to  choose  the  correct  locations  of  the  plastic  hinges, 
the  method  essentially  is  based  on  trial  and  error.  If  Mp  is  not 
exceeded  at  any  point  in  the  moment  diagram,  but  is  attained  at  a 
sufficient  number  of  points  to  form  a  mechanism,  the  proper  plastic 
moment  diagram  has  been  found.  If  a  moment  greater  than  Mp  appears  at 
any  point,  it  becomes  necessary  to  assume  another  moment  distribution. 
When  the  correct  moment  diagram  has  been  found,  the  plastic  limit 
load  may  be  evaluated  by  statics. 


The  mechanism  method  is  often  less  cumbersome  than  the  equili¬ 
brium  method.  It  consists  of  equating  the  external  work  done  by  the 
loads  acting  on  the  deflected  structure  to  the  internal  work  performed 
by  plastic  hinge  rotations  during  a  virtual  displacement  of  an 
assumed  mechanism.  As  in  the  case  of  the  equilibrium  method,  a  trial 
and  error  procedure  is  needed  where  there  are  two  or  more  possible 
mechanisms.  The  mechanism  which  yields  the  smallest  collapse  load 
governs  the  solution  of  the  problem. 

The  above  methods  predict  the  failure  load,  but  give  no 
indication  of  the  load-deformation  characteristics  of  the  structure. 

The  methods  also  are  difficult  to  use  for  large  frames,  since  the 
number  of  possible  mechanisms  may  be  very  large. 

The  complete  load-deformation  response  of  complex  frames  can 
be  determined  by  a  computerized  step-by-step  analysis  (Wang  1963) . 

In  this  approach  the  difficulty  of  locating  and  verifying  plastic 
mechanisms  for  large  structures  is  avoided.  During  the  first  step 
a  linear  elastic  computation  is  performed,  based  on  the  original 
geometry  and  stiffness  of  the  structure.  The  search  for  the  location 
of  the  first  plastic  hinge  is  conducted,  and  the  loads  on  the 
structure  are  increased  until  the  section  with  the  highest  moment 
has  reached  the  plastic  moment  of  the  member  at  that  point.  The  first 
plastic  hinge  now  has  formed,  and  the  stiffness  of  the  structure  must 
be  modified  by  changing  the  particular  member  stiffness.  The  search 
for  the  next  hinge  to  be  formed  is  then  conducted  by  the  same  procedure, 
which  is  continued  until  a  plastic  mechanism  has  been  formed. 

Two  different  step-by-step  approaches  can  be  utilized.  One 
is  an  elastic-plastic  analysis,  where  elastic  behavior  is  assumed 
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until  the  hinge  has  formed,  and  the  moment-curvature  relationship 
shown  in  Fig.  2-lla  is  used.  According  to  this  assumption,  the 
segments  of  the  frame  between  the  hinges  will  remain  elastic.  The 
resulting  load-deformation  characteristic  is  termed  the  first  order 
elastic-plastic  curve,  and  is  illustrated  by  curve  E  in  Fig.  1-1. 

The  second  approach  based  on  the  step-by-step  method  is 
a  rigid-plastic  analysis,  which  gives  a  first  estimate  of  the 
ultimate  load  of  the  frame  (Beedle  1958).  According  to  this  theory 
the  frame  remains  undeformed  until  a  mechanism  has  formed.  This  is 
consistent  with  the  moment-curvature  relationship  assumed  in 
Fig.  2-llb.  The  behavior  is  termed  first  order  rigid-plastic  and 
is  illustrated  by  curve  C  in  Fig.  1-1. 

2 . 5  Inelastic  Nonlinear  Analysis 

The  determination  of  the  ultimate  strength  of  a  frame  is  a 
nonlinear  bending  problem  (Alvarez  and  Birnstiel  1969).  The  non¬ 
linearity  results  from  changes  in  stiffnesses  of  the  members  due  to 
axial  force  and  plastif ication,  and  from  changes  in  the  geometry  of 
the  structure  due  to  displacements.  The  true  ultimate  strength  has 
been  reached  when  the  combination  of  yielding,  axial  force,  and  joint 
displacements  reduce  the  stiffness  of  the  structure  to  such  an  extent 
that  the  frame  cannot  carry  any  more  load. 

Experiments  have  shown  that  the  ultimate  strength  is  less 
than  that  predicted  by  simple  plastic  theory  based  on  the  plastic 
hinge  concept.  The  actual  and  the  predicted  values  differ  because 
nonlinear  effects  are  neglected  in  the  simple  plastic  theory. 


Moment,  M  Moment, 


J  .• 


(a)  Elastic-Plastic  Material 


Mp 


- ► 

Curvature,  0 


(b)  Rigid-Plastic  Material 


FIGURE  2-1 1  Idealized  Moment-Curvature  Relationship 


2.5.1  The  Step-by-Step  Method  of  Analysis 


The  step-by-step  method  described  by  Korn  and  Galambos  (1968) , 
and  outlined  in  Sect.  2.4.5,  can  be  modified  to  account  for  nonlinear 
effects.  In  this  case  the  analysis  is  formulated  on  the  deformed 
structure  thus  incorporating  both  strength  and  stability  effects  into 
the  analysis.  This  can  be  done  by  performing  a  P-A  analysis, 
discussed  in  Sect.  2.3.2,  on  a  frame  with  plastic  hinges. 

If  the  moment-curvature  relationship  is  of  the  type  shown  in 
Fig.  2-lla,  indicating  an  elastic-plastic  material,  a  second  order 
elastic-plastic  analysis  can  be  performed.  Such  a  model  can 
adequately  represent  complex  inelastic  behavior  within  the  limita¬ 
tions  imposed  by  confining  inelastic  effects  to  individual  cross 
sections.  The  results  of  the  analysis  is  the  complete  load-deformation 
response  of  a  frame,  and  is  shown  as  curve  F  in  Fig.  1-1.  A  second 
order  rigid-plastic  analysis  can  be  performed  on  the  frame  for  which 
the  moment-curvature  relationship  for  the  material  is  as  illustrated 
in  Fig.  2-llb.  In  this  approach  the  frame  remains  undeformed  until 
a  mechanism  is  developed,  at  which  point  the  frame  exhibits  lateral 
deformation  with  a  decrease  in  load.  This  behavior  is  shown  in 
Fig.  1-1  as  curve  D.  The  plastic  collapse  load  obtained  from  the 
rigid-plastic  analysis  will  always  give  an  upper  bound  value  for  the 
true  maximum  load. 

2.5.2  Development  of  Second  Order  Plastic  Analysis 

A  technique  for  the  inelastic  stability  analysis  of  portal 
frames  with  hinged  bases,  neglecting  the  effect  of  axial  forces,  has 
been  presented  by  Ang  (1960).  Ojalvo  and  Lu  (1961)  developed  a 


procedure  for  the  analysis  of  portal  frames,  symmetrically  loaded 
into  the  inelastic  region.  Moses  (1964)  showed  a  procedure  for  the 
inelastic  analysis  of  portal  frames  with  hinged  bases,  considering 
the  spread  of  the  inelastic  zones  and  the  effects  of  axial  forces. 

In  this  procedure  the  frame  is  given  a  lateral  displacement  and  the 
compatible  forces  and  moments  are  determined  by  iterations.  Adams 
(1964)  presented  a  similar  method,  but  neglected  the  axial  force  in 
the  beam,  as  did  Chu  and  Pabarcius  (1964).  The  latter  also  assumed 
elastic  behavior  of  the  beam. 

The  effect  of  lateral  joint  displacements  on  the  inelastic 
behavior  of  a  hinged  base  portal  frame  was  included  in  a  paper  by 
Yura  and  Galambos  (1965).  The  influence  of  beam  axial  force  was 
neglected.  This  is  the  method  called  the  P-A  approach,  as  discussed 
in  Sect .  2.3.2. 

Horne  and  Majid  (1966)  incorporated  both  design  and  analysis 
features  into  a  second  order  program  capable  of  considering  reductions 
in  plastic  moment  capacity  due  to  axial  forces,  but  restrictions  are 
placed  on  the  column  plastic  hinges.  The  program  is  limited  to 
considering  only  bending  deformations. 

Davison  and  Adams  (1971)  considered  the  effect  of  a  finite 
column  width  on  the  lateral  stiffness  of  the  structure*  Considering 
the  width  of  vertical  members  reduces  the  clear  span  of  girders,  thus 
increasing  the  bending  stiffness  and  decreasing  the  fixed  end  moments. 

McNamee  and  Lu  (1972)  examined  the  frame  buckling  problem, 
and  presented  an  engineering  solution  for  inelastic  buckling  loads  of 
multistory  frames.  In  this  approach  only  flexural  deformations  are 
considered,  and  small  deformation  theory  is  used.  The  plastic  hinges 
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are  assumed  to  take  place  at  discrete  points,  and  an  idealized  elastic- 
plastic  behavior  is  implied. 

Davison  and  Adams  (1974)  developed  a  method  to  analyze  the 
stability  of  braced  and  unbraced  frames.  The  formulation  is  based 
on  the  slope-deflection  equations,  modified  to  incorporate  the  effect 
of  plastic  hinging  and  finite  column  widths.  The  material  is  assumed 
elastic-perfectly  plastic,  and  P-A  effects  are  covered.  The  plastic 
moment  capacity  of  the  column  is  reduced  for  axial  force,  and  the 
effect  of  axial  force  on  the  stiffness  and  carry-over  factors  for  a 
column  is  considered  by  using  elastic  stability  functions. 

Frame  collapse  analysis  has  been  presented  by  Tranberg  et  at. 
(1976).  The  concepts  of  geometric  and  tangent  stiffness  are  used, 
together  with  an  initial  strain  procedure.  An  attempt  was  made  to 
incorporate  most  of  the  factors  that  influence  frame  collapse  into 
a  single  analysis.  Extended  plastic  zones,  axial  deformations,  and 
strain  hardening  are  considered,  but  residual  stresses  are  ignored. 

Cohn  and  Rafay  (1977)  investigated  second  order  analysis 
using  the  P-A  method.  The  analysis  made  use  of  the  plastic  hinge 
concept,  as  well  as  considering  the  axial  forces  in  the  members. 
Residual  stresses  and  strain  hardening  were  neglected,  and  only  stress 
distribution  at  collapse  is  obtained,  with  no  information  about  the 
complete  load-deformation  characteristics. 


CHAPTER  III 


ELASTIC  FORMULATIONS  AND  SOLUTIONS 

In  Chapter  2  a  number  of  different  kinds  of  classical  formu¬ 
lations  were  reviewed  for  the  solution  of  member  and  frame  stability 
problems.  Relatively  recently,  the  finite  element  method  has  proven 
itself  as  one  of  the  most  versatile  techniques  available  for  the 
numerical  solution  of  complex  structural  problems  (Gallagher  and 
Padlog  1963) .  In  this  chapter  the  elastic  solutions  of  frame  stability 
problems,  using  the  finite  element  method,  are  compared  with  those 
obtained  by  the  classical  solutions,  reviewed  in  Chapter  2. 

The  finite  element  method  as  a  solution  technique  is  reviewed 
in  Sect.  3.1.  Section  3.2  is  devoted  to  the  basic  equations  and 
assumptions.  The  formulation  of  stability  equations  of  equilibrium 
is  derived  in  Sect.  3.3  together  with  a  solution  technique  to  solve 
for  the  elastic  critical  load  of  frames.  Numerical  results  for  a 
variety  of  problems  are  presented  in  Sect.  3.4.  Formulation  of  non¬ 
linear  equilibrium  equations  as  well  as  finite  element  models  are 
developed  in  Sect.  3.5,  which  also  contains  a  solution  technique  to 
solve  for  the  load-deformation  response  of  a  framed  structure  using 
the  finite  element  equations.  A  variety  of  nonlinear  problems  are 
checked,  for  accuracy  and  efficiency  of  the  solution  technique,  in 
Sect.  3.6. 
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3 . 1  Review  of  the  Finite  Element  Method 

The  development  of  the  finite  element  method  as  an  analysis 
tool  was  initiated  with  the  advent  of  the  electronic  digital  computer. 
In  the  solution  of  a  continuum  problem  it  is  necessary  to  establish 
and  solve  a  system  of  governing  differential  equations.  Using  the 
finite  element  method  and  a  digital  computer  it  is  possible  to  discre¬ 
tize  and  solve  the  governing  equations  for  complex  system  in  a  very 
effective  way  (Zienkiewicz  1971). 

In  the  finite  element  method  structures  can  be  visualized 
as  an  assemblage  of  structural  elements  interconnected  at  a  discrete 
number  of  nodal  points.  In  an  elastic  continuum  the  true  number  of 
interconnection  points  is  infinite  and  therein  lies  the  difficulty 
of  its  numerical  solution.  The  concept  of  finite  element  as  originally 
introduced  by  Turner  et  at.  (1956)  attempts  to  overcome  this 
difficulty  by  assuming  the  real  continuum  to  be  divided  into  elements 
interconnected  only  at  a  finite  number  of  nodal  points  at  which  work 
equivalent  fictitious  forces  are  introduced.  If  such  an  idealization 
is  possible  the  problem  reduces  to  that  of  a  conventional  structural 
type  well  suited  to  numerical  treatment. 

Although  finite  element  formulations  can  be  based  on  either 
stress  fields  or  displacement  fields,  most  often  a  displacement  based 
finite  element  formulation  is  applied  in  practice  since  it  can  easily 
be  programmed  for  digital  computers.  The  procedure  of  the  finite 
element  method  as  applied  to  a  frame  can  be  described  by  the  following 
steps  (Chen  and  Atsuta  1977) . 

(a)  Each  member  is  separated  by  imaginary  lines  into  a  number  of 


finite  elements. 


(b)  The  elements  are  assumed  to  be  interconnected  at  a  discrete  number 
of  nodal  points.  The  displacements  of  these  nodal  points  become 
the  unknown  parameters  of  the  problem,  as  in  simple  structural 
analysis . 

(c)  A  set  of  functions  is  chosen  to  uniquely  define  the  state  of 
displacement  within  each  finite  element  in  terms  of  its  nodal 
displacements.  These  functions  are  so  chosen  that  they  ensure 
continuity  throughout  the  member. 

(d)  The  displacement  functions  uniquely  define  the  state  of  strain 
within  an  element  in  terms  of  nodal  displacements.  These  strains 
together  with  any  initial  strain  define  the  state  of  stress 
throughout  the  element  and  also  at  its  boundaries. 

(e)  A  system  of  forces  concentrated  at  the  nodes  of  the  element, 
equilibrating  the  boundary  stresses  and  any  distributed  loads, 
is  determined.  This  results  in  a  stiffness  relationship  of 
the  form 

[k]  {q}  =  {Q}  (3.1.1) 

in  which  { q }  is  the  vector  of  nodal  displacements,  {Q}  is  the 
vector  of  nodal  forces,  and  [k]  is  the  element  stiffness  matrix 
whose  typical  influence  coefficient  k  is  the  force  Q_^  due  to  a 
unit  displacement  q_.  when  all  other  displacements  (q's)  are  zero. 

(f)  Once  this  stage  has  been  reached  the  solution  procedure  can 
follow  the  standard  direct  stiffness  structural  analysis. 

Variational  principles  may  be  regarded  as  one  of  the  most 
important  bases  for  the  finite  element  method  (Washizu  1969).  They 
have  contributed  to  the  development  of  structural  analysis  by  leading 
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to  rigorous  finite  element  formulations.  Numerous  types  of  finite 
element  models  may  be  derived  based  on  variational  principles.  Since 
continuum  problems  cannot  usually  be  solved  ’exactly’,  the  variational 
method  provides  an  approximate  formulation  of  the  problem  which  yields 
a  solution  compatible  with  the  assumed  degree  of  approximation.  The 
variational  technique  used  herein,  for  the  formulation  of  equations 
such  as  Eq.  3.1.1,  is  the  principle  of  virtual  work  (Fung  1965). 

3 . 2  Basic  Equations  and  Assumptions 

The  following  assumptions  are  used  to  formulate  the  basic 
equations  of  beams  and  beam-columns  throughout  this  work,  unless  other¬ 
wise  stated. 

1  -  The  member  is  straight,  prismatic  and  symmetric  about  the  plane 

of  the  frame  (the  y-z  plane) .  The  member  z  axis  (the  reference 
axis)  coincides  with  the  centroidal  axis  of  the  cross-section. 

2  -  Loads  are  applied  in  the  plane  of  the  frame  only  at  the  ends  of 

an  element. 

3  -  Only  in-plane  deformations  occur  and  sections  that  were  originally 

normal  to  the  z  axis  will  remain  undistorted  and  normal  to  the 
beam  axis  after  deformation. 

4  -  The  slope  at  any  point  along  the  reference  axis  is  given  by 

Av 

v’  =  -T-2-  =  sin6  (3.2.1) 

o  Az 

Assumption  4  is  basic  to  the  nonlinear  formulation  developed 
herein.  Since  a  Lagrangian  coordinate  system  is  used,  this  expression 
is  'exact'  when  the  elemental  length  of  beam  (Az)  does  not  change  in 
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length.  Since  the  axis  of  a  member  may  be  expected  to  undergo  only 
very  small  deformations  up  to  the  point  of  collapse  (of  the  order  of 
1  or  2  percent),  Eq.  3.2.1  permits  accurate  solutions  for  large 
displacement  problems  of  frames,  as  will  subsequently  be  demonstrated. 

The  above  assumptions  permit  the  displacements  u  and  v  of 
an  arbitrary  point  A  on  a  beam  cross  section  to  be  expressed  in  terms 
of  the  displacements  of  the  reference  axis  of  the  beam.  Thus,  referring 
to  Fig.  3-1,  the  displacements  can  be  written  as 

u  =  u  -  y  sin0  (3.2.2a) 

o 

and 


v  =  v  -  y  (1  -  cos0) 

Assumption  4  permits  Eqs.  3.2.2  to  be  written  as 


and 


u  =  u  -  y  v’ 
o  o 


(3.2.2b) 


(3.2.3a) 


v  =  v  -  y  (1  -  cos@)  (3.2.3b) 

The  axial  strain  at  the  arbitrary  point  A  may  now  be  obtained 
from  the  large  displacement  strain-displacement  equations  (Saada 
1974)  as 


£  =  U'  +  i  [(U')2  +  (V’)2] 

Z  Z 


(3.2.4) 
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FIGURE  3-1  Element  Deformation 
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Substitution  of  Eqs.  3.2.3  into  Eq.  3.2.4  permits  the  strain  to  be 
determined  throughout  the  beam  in  terms  of  displacements  of  the 
reference  axes. 

A  variety  of  beam  formulations  are  now  available  depending 
on  the  degree  to  which  Eqs.  3.2.3  and  3.2.4  are  approximated  during 
the  formulation.  The  assumptions  that  yield  the  classical  linear 
stability  equations  are  followed  in  Sect.  3.3.  However,  a  large 
displacement  formulation,  which  retains  all  terms  implied  by  Eqs. 

3.2.1,  3.2.3  and  3.2.4,  is  developed  in  Sect.  3.5. 

3 . 3  Formulation  of  Linear  Stability  Equations 

In  this  section  a  linear  stability  formulation  is  derived 
and  a  solution  technique  to  solve  for  the  elastic  critical  load  of 
plane  frame  is  described. 

Detailed  derivation  of  the  equations  of  this  section  is 
carried  out  in  Appendix  A. 

3.3.1  Secant  Stiffness  for  Elastic  Element 

Within  the  limitation  of  "small"  displacements  and  elastic 
stresses  the  equilibrium  equation  for  an  element,  as  shown  in  Fig. 

3-2a,  can  be  written  as  (Ghali  and  Neville  1972) 

[k]  {q}  =  {Q}  (3.3.1) 

in  which  {q}  is  the  vector  of  nodal  displacements,  (Q }  is  the  vector 
of  nodal  forces,  and  [k]  is  the  element  stiffness  matrix  given  by 


Eq.  A. 2. 6c. 
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(a)  Local  Nodal  Displacements  and  Forces 
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FIGURE  3-2  Nodal  Displacements  and  Forces 
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When  the  element  is  subjected  to  both  axial  force  and  bending 
the  stiffness  of  the  element,  which  in  this  case  may  be  called  the 
secant  stiffness  matrix,  becomes  a  function  of  the  axial  load.  The 
equation  of  equilibrium,  Eq.  3.3.1,  in  this  case  becomes  (Chajes  1974) 

[k  ]  {q}  =  {Q}  (3.3.2) 

s 

in  which 

[k  ]  =  [k]  +  [k  ]  (3.3.3) 

The  elastic  stiffness  [k]  is  that  commonly  used  for  frame  analysis 
(Eq.  A. 2. 6c).  For  finite  element  analysis  [k]  is  usually  based  on  a 
cubic  displacement  assumption  with  shear  deformations  ignored.  The 
matrix  [k  ]  is  called  the  geometric  stiffness  matrix  (Eq.  A.2.6d) 

o 

and  is  a  function  of  the  element  axial  force. 

3.3.2  Structure  Secant  Stiffness  Matrix 

The  vector  of  element  local  nodal  displacements  { q }  can  be 
expressed  in  terms  of  the  element  global  nodal  displacements  (r  }  , 

E  Cj 

shown  in  Fig.  3-2b,  as 

{q}  =  [T]  (rE>G  (3.3.4) 


in  which  [T]  is  a  transformation  matrix  given  in  Table  3-1. 


If 


'Vg1 


is  the  vector  of  nodal  forces  associated  with  nodal 


displacements  {r  }  then,  by  the  reasoning  of  Sect.  2.2 

E  Cj 


<q>  {Q} 


<rE>G 


(3.3.5) 


c  s  0 

s  -c  0 


0  0-1 


[T]  = 


c  s  0 

s  -c  0 

0  0-1 


C  =  COS  (p 

s  =  sin  (f) 


TABLE  3-1  Transformation  Matrix 
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Substituting  Eq.  3.3.4  into  Eq.  3.3.5  yields 

"Vg  [t]T  {q}  “  "Vo  {re}g  (3-3-6> 

or 

<rE>G  {[TlT  *  {REJg}  =  0  (3-3.7) 

Since  this  must  be  true  for  all  <r  >  ,  Eq.  3.3.7  implies  that 

E  G 

{Re)g  =  [T] T  {Q}  (3.3.8) 

Premultiplying  Eq.  3.3.2  by  [T]T  yields 

[T]T  [k  ]  {q}  =  [T]T  {Q}  (3.3.9) 

Substituting  Eqs .  3.3.4  and  3.3.8  into  Eq.  3.3.9  yields 

[T]x  [k  ]  [T]  {r}„  =  {R_}„  (3.3.10a) 

S  EG  EG 

or 

[VG  {rE}G  ■  {RE}G  (3.3.10b) 

The  formation  of  the  complete  stiffness  matrix  for  the  entire  structure 
is  obtained  by  direct  addition  for  all  the  interface  nodes.  The 
equilibrium  equations  for  the  entire  structure  can  then  be  written  as 

[K  ]  {r}  =  {R}  (3.3.11) 

s 


in  which  {r}  is  the  assembled  vector  of  global  nodal  displacements, 

{ R}  is  the  assembled  vector  of  global  nodal  forces,  and  [K  ]  is 

s 

the  structure  secant  stiffness  matrix  and  is  a  function  of  the 
element  axial  forces.  The  matrix  [Kg]  consists  of  two  parts  as 


65 


[K  ]  =  [K]  +  [K  ]  (3.3.12) 

s  G 

in  which  [K]  is  the  structure  stiffness  matrix,  and  [K  ]  is  the  geometric 

G 

stiffness  matrix  assembled  for  the  overall  structure. 

Equation  3.3.11  is  an  inhomogeneous  equilibrium  equation  that 
can  be  used  to  obtain  a  unique  solution  for  the  frame  behavior  under  a 
set  of  loads  {R}.  In  addition,  the  equation  can  be  used  to  obtain  the 
frame  critical  load  by  the  reasoning  discussed  in  Sect.  2.2.2. 

3.3.3  Frame  Critical  Load 

Solving,  symbolically,  for  the  displacement  vector  {r}  in 
Eq.  3.3.11  yields 

{r}  =  [K  ]-1  {R}  (3.3.13) 

s 

The  displacement  vector  {r}  increases  without  limit,  for  finite 
values  of  {r},  only  when  the  inverse  of  the  secant  stiffness  matrix 
becomes  infinitely  large.  Since  the  inverse  of  the  matrix  is  obtained 
by  dividing  the  adjoint  matrix  by  the  determinant,  the  elements  of  the 
inverse  will  become  indefinitely  large  when  the  determinant  is  equal 
to  zero.  Thus 

I  I  K  11=0  (3.3.14) 

ii  s  ii 

gives  the  stability  criterion  as  discussed  in  Sect.  2.2.2. 

Equation  3.3.14  can  be  rewritten  as 


||  [K]  +  [Kg]  | 


0 


(3.3.15) 
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In  dealing  with  linear  stability  problems  it  is  assumed  that  the 
buckling  deformations  are  independent  of  all  deformations  prior  to 
instability.  This  leads  to  the  possibility  of  expressing  the  axial 
force  for  each  element,  at  the  buckling  condition,  as  a  multiple  of 
the  axial  force  occurring  for  some  reference  loading.  Since  the 
critical  load  is  unknown  a  factor  X,  an  arbitrary  multiple  of  the 
reference  load  vector  {r},  is  introduced  to  represent  the  relative 
magnitude  of  the  applied  loads.  Thus 

{R}  =  X  {R}  (3.3.16) 

The  factor  X  is  called  the  instability  factor  or  eigenvalue.  Since 
the  geometric  stiffness  matrix,  derived  from  Eq.  A.2.6d,  is  propor¬ 
tional  to  the  internal  axial  forces  it  follows  that 

[Kg]  =  X  [Kg]  (3.3.17) 


in  which  [K  ]  is  the  geometric  stiffness  matrix  for  the  reference 

Cj 

value  of  the  applied  loading.  Thus  Eq.  3.3.15  can  be  written  as 

||  [K]  +  X  [Kg]  ||=0  (3.3.18) 

which  is  the  requirement  for  nontrivial  solution  of  the  homogeneous 
form  of  Eq.  3.3.11  and  gives  the  eigenvalue,  X  ,  that  represents 
the  critical  load. 

Equation  3.3.18  can  be  solved,  using  the  determinant  search 
technique  discussed  in  Sect.  2.2.2,  for  the  critical  load,  P^. 

However,  since  Eq.  3.3.18  is  a  classical  eigenvalue  problem  an 
inverse  power  iteration  method  becomes  more  practical. 
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3.3.4  Inverse  Power  Iteration  Method 

The  inverse  power  iteration  method  is  described  by  Householder 
(1964)  and  Wilkinson  (1965) .  It  is  one  of  the  most  efficient  methods 
available  for  solving  for  the  lowest  eigenvalue  in  the  case  of  large 
systems.  It  has  the  advantage  of  providing  the  buckling  eigenvalue 
A  and  the  associated  buckling  mode  simultaneously.  Stability  eigen- 
systems  of  over  20,000  degrees  of  freedom  have  been  handled  efficiently 
with  the  use  of  this  method  (Skogh  et  at.  1972) . 

To  solve  for  the  critical  load  of  a  frame  the  homogeneous 
form  of  Eq.  3.3.11  is  used.  That  is 

[K  ]  {r}  =  {0}  (3.3.19) 

s 

or 

[  [K]  +  [Kg]  ]  (r)  =  {0}  (3.3.20) 

Substituting  Eq.  3.3.17  into  Eq.  3.3.20  yields 

[  [K]  +  A[Kg]  ]  {r}  =  {0}  (3.3.21) 

which  can  be  written  in  the  form 

[K]  {r}  =  -X[KG]  {r}  (3.3.22) 

To  solve  for  the  lowest  value  of  X,  Eq.  3.3.22  can  be  rearranged  as 

j  [K]  {r}  =  -  [KGHr>  (3.2.23) 

At  this  stage  the  inverse  power  iteration  method  is  performed, 

using  'Cholesky  Decomposition'  (Elwi  and  Murray  1977),  to  solve  for 

the  eigenvalue  A  and  the  associated  buckling  mode  (r).  The  method 

cr 
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can  be  described  by  the  following  procedure. 

1  -  The  structure  elastic  stiffness  matrix  [K]  is  assembled. 

2  -  For  a  particular  set  of  loads,  a  simple  linear  elastic  analysis 

is  performed  and  the  axial  force  in  each  element  is  determined. 

3  -  Using  the  axial  forces  obtained  in  Step  2,  the  geometric  stiffness 

matrix  [K  ]  can  be  assembled. 

Cj 

4  -  An  approximate  normalized  buckling  mode  {r}  is  assumed. 

5  -  By  substituting  the  vector  {r}  into  the  right  hand  side  of  Eq. 

3.3.23,  the  equation  can  be  rewritten  as 


[K]  if}  =  {Rg} 


(3.3.24a) 


in  which 


[R  ]  =  -  [K„]  {r} 

g  b 


(3.3.24b) 


6-  Equation  3.3.24a  can  be  put  in  the  form 


[L]  [D]  [L]T  {f} 


(3.3.25) 


or 


[L]  [G]  {f} 


(3.3.26a) 


in  which 


[G]  =  [D]  [L] 


T 


(3.3.26b) 


The  solution  is  then  obtained  as 


if}  -  [G]-1  {?} 


(3.3.27a) 
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in  which 

{r}  =  [L]"1  {R  }  (3.3.27b) 

c3 

By  forward  substitution  one  can  solve  Eq.  3.3.27b,  for  the  vector 
{r}.  Then  substituting  {r}  into  Eq.  3.3.27a  and  back  substitution 
yields  {y} . 

7  -  The  largest  value  in  the  new  vector  is  determined  and  the  vector 

is  divided  by  this  value  to  obtain  a  new  normalized  buckling  mode 
{r}. 

8  “  The  new  buckling  mode  {r},  obtained  from  Step  7,  is  compared  with 

the  previous  one. 

9  -  Steps  5  to  8  are  repeated  using  the  new  vector  {r},  until  the 

difference  between  the  buckling  modes  from  two  successive  iterates 
is  arbitrarily  small. 

10  -  The  largest  value  in  vector  {r},  obtained  in  Step  7  after  convergence, 
is  equal  to  and  the  critical  load  on  the  frame  is  obtained  by 
multiplying  the  particular  set  of  loads  assumed  in  Step  2  by  the 
inverse  of  this  value. 

3.3.5  Computer  Program 

A  computer  program,  called  Eigenvalue  Stability  Analysis  of 
_Frames  (ESTANF) ,  has  been  developed  to  perform  an  eigenvalue  analysis 
based  on  the  formulation  presented  in  this  section.  The  results  for 
some  example  applications  are  presented  in  the  next  section. 


. 
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3 . 4  Numerical  Results  of  Elastic  Buckling  Problems 

The  formulation  given  in  Sect.  3.3  enables  the  solution  of  a 
wide  variety  of  problems.  At  this  stage  it  is  appropriate  to  examine 
the  numerical  accuracy  that  may  be  attained  in  the  solution  and  also 
to  check  the  efficiency  of  the  computer  program.  In  general  the  ade¬ 
quacy  and  validity  of  numerical  formulation  may  be  measured  by  comparing 
its  performance  on  problems  for  which  accurate  solutions  have  been 
derived  by  classical  methods. 

3.4.1.  Sample  Solutions 

In  order  to  demonstrate  the  efficiency  of  the  finite  element 
method  in  the  solution  of  linear  stability  problems  several  examples 
are  solved  in  this  section  and  compared  with  classical  solutions.  For 
each  problem  the  cpu  computer  time  is  given  as  a  measure  of  the  cost 
of  running  the  program. 

E-l  Critical  Buckling  Load  of  Columns 

Three  column  types,  shown  in  Fig.  3-3,  are  analyzed.  Each 
type  is  solved  twice,  by  dividing  the  member  into  two  elements  and 
then  four  elements.  As  can  be  observed  from  the  results,  in  Table  3-2, 
the  difference  between  the  classical  solutions  given  by  Timoshenko  and 
Gere  (1961)  and  the  present  analysis  is  less  than  1.0%.  For  all  types 
of  columns  analyzed  the  cpu  time  ranges  between  .15  and  0.2  seconds. 

E-2  One  Story  -  One  Bay  Example  Frame 

Each  of  the  frames  shown  in  Fig.  3-4  is  solved  using  different 
values  of  G  where  G  is  given  by 


' 


mZ 

(a) 

(b) 

(c) 

L  =  120  in 

I  =  100  in4 

E  =  30,000  ksi 

FIGURE  3-3  Types  of  Columns 
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r  pi 


E  =  30,000  ksi  Frame  (a) 

L  =  h  =  1 20  in 

Ic=  100  in4 

lb  =  50,100  or  200  in4 

ip  pi 


Frame  (b) 


FIGURE  3-4  One  Storey  -  One  Bay  Frame 
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G 


I  /L 
c  c 

I  /L 
g  g 


(3.3.1) 


in  which  I  and  I  are  the  moment  of  inertia  of  the  column  and  beam, 
c  g 

and  and  are  the  length  of  column  and  beam  respectively.  For  each 
value  of  G  the  frame  is  solved  twice,  first  assuming  the  member  as  one 
element  and  second  by  dividing  each  member  into  three  elements  with 
the  central  element  being  80%  of  the  length  of  the  member.  The  reason 
for  the  unusual  subdivision  with  three  elements  will  be  given  later  in 
the  thesis,  when  inelastic  behavior  is  considered. 

The  results  for  frame  (a) ,  with  hinged  bases,  is  given  in 
Table  3-3  while  the  numerical  results  for  frame  (b) ,  with  fixed  bases, 
is  shown  in  Table  3-4.  The  differences  in  critical  loads  between  the 
present  analysis  and  the  classical  solutions  tabulated  by  Lu  (1962)  are 
less  than  1%  even  when  treating  each  member  as  one  element.  It  is 
noted  that  the  present  analysis  gives  an  upper  bound  to  the  classical 
solution.  The  cpu  time  ranges  between  0.25  to  0.35  seconds. 

E-3  Two  Story  -  Two  Bay  Example  Frame 

The  frame  shown  in  Fig.  3-5  is  solved  treating  each  member  as 
one  element.  An  energy  formulation  derived  by  Johnson  (1960)  gives  a 
critical  load  of  780.2  kips  while  the  present  formulation  provides 
a  critical  load  of  778.4  kips  with  a  difference  of  0.23%.  Unfortunately 
there  was  not  an  exact  solution  available  for  this  problem.  The  computer 
cpu  time  is  about  0.5  sec. 

E-4  Multistory  Multibay  Example  Frame 

A  classical  method  of  analysis  based  on  the  energy  formulation, 
for  the  large  frame  shown  in  Fig.  3-6,  is  presented  by  Johnson  (1960). 


- 


75 


<3 

co 

a 

3 

3 

ll 

54 

00 

CO 

3 

CNI 

r— 1 

rH 

4-1 

• 

• 

• 

co 

4-4 

o 

o 

o 

4-1 

•H 

c 

Q 

W 

4-1 

m 

rH 

m 

o 

54 

• 

• 

• 

3 

CTn 

o 

43 

CO 

• 

cu 

co 

CO 

ON 

•H 

o 

MT 

CO 

CM 

CO 

r— 1 

<3 

z 

c 

< 

3 

u 

r— 1 

3 

3 

3 

3 

3 

CO 

ii 

3 

cn 

r-» 

3 

3 

CO 

rH 

54 

44 

• 

• 

• 

P4 

co 

44 

o 

o 

o 

4J 

•H 

c 

Q 

3 

e 

0) 

w 

4-1 

CO 

m 

40 

o 

u 

• 

• 

• 

3 

o 

o 

40 

• 

C4 

<r 

00 

O0 

o 

CO 

CM 

2 

r— 1 

3 

3 

a 

o 

•H 

•H 

CO 

mt 

rH 

CO 

4-1 

3 

• 

• 

• 

CO 

3 

3 

00 

ON 

40 

3 

rH 

PL, 

CO 

r^. 

cn 

r— 1 

o 

Ml- 

co 

CM 

CJ 

CO 

m 

O 

O 

o 

• 

• 

• 

o 

rH 

CM 

3 


3 

I 

5-1 

p4 


4-1 

o 

co 


3 

CO 

a) 


3 

a 

•H 

a) 


c**) 

i 

CO 

W 

CQ 

< 

E-4 


76 


0) 

cn 

a 

d 

OJ 

II 

u 

CN 

o 

OJ 

m 

in 

4-1 

• 

• 

• 

CO 

4-1 

o 

o 

O 

4-1 

•H 

d 

Q 

QJ 

6-2 

w 

4-4 

<r 

i"- 

r-* 

o 

5m 

• 

• 

• 

a 

m 

<r 

cn 

co 

• 

Pm 

in 

in 

•H 

o 

r-'- 

m 

CN 

CO 

Sn 

rH 

d 

2 

rH 

rH 

rH 

d 

< 

a) 

■u 

r— 1 

a 

d 

d 

O) 

OJ 

CO 

n 

d 

m 

o 

CTi 

a; 

QJ 

cn 

a> 

CTi 

Jm 

4—t 

• 

• 

• 

Pm 

CO 

4M 

o 

o 

O 

hi 

•H 

d 

Q 

OJ 

E 

a> 

w 

4-1 

rH 

CJl 

o 

5m 

• 

• 

• 

CJ 

cn 

• 

Pm 

r>. 

<r 

m 

o 

m 

CN 

2 

rH 

rH 

rH 

rH 

d  d 
o  o 

cn 

o 

o 

•H  *H 

• 

• 

• 

CO  4-1 

*4 

in 

m 

m 

CO  3 

J 

m 

cn 

m 

CTJ  rH 

Pm 

r-- 

in 

CN 

t— 1  o 
o  cn 

rH 

rH 

rH 

m 

o 

o 

O 

• 

• 

• 

o 

rH 

CM 

77 


i1 


P  1 2  P  P 

21  i  21 


l 


I  =  75  in4 
L  =  96  in 
E  =  30,000  ksi 


FIGURE  3-5  Two  Storey  -  Two  Bay  Frame 
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FIGURE  3-6  Multistory  -  Multibay  Frame 
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It  gives  a  critical  load  of  1189.6  kips  compared  to  1187.4  kips 
obtained  from  the  present  analysis  furnishing  a  difference  of  0.2%. 

No  exact  solution  was  available.  The  cpu  computer  time  is,  approxi¬ 
mately,  1  second. 

3.4.2  Comments  on  Results 

The  examples  presented  in  this  section  have  been  selected  to 
test  the  ability  of  the  computer  program  to  analyze  a  variety  of 
rectilinear  framed  structures.  The  analysis,  as  it  can  be  seen  from 
the  previous  examples,  exhibits  good  results. 

3 . 5  Formulation  of  Nonlinear  Equations 

In  this  section  a  technique  for  nonlinear  stability  formula¬ 
tion  of  elastic  multistory  frames  is  presented.  The  formulation  of 
nonlinear  equilibrium  equations  is  derived,  based  on  the  principle  of 
virtual  work.  Nonlinearities  enter  the  formulation  as  a  result  of  non 
linear  strain-displacement  relationships,  which  contain  strain  product 
of  the  same  order  of  magnitude  as  the  engineering  strains. 

The  Newton-Raphson  method  is  used  to  solve  for  the  overall 
load-displacement  characteristics  of  multistory  building  frames. 

3.5.1  Description  of  Basic  Formulation 

The  basic  equations  from  which  a  nonlinear  formulation  may 
be  derived  have  been  given  in  Sect.  3.2.  In  Appendix  A  approximations 
were  made  which  produced  a  linear  stability  formulation.  An  effort  is 
made  in  this  section  to  retain  all  nonlinear  terms.  Retaining  the 
(1  -  cos0)  term  in  Eq.  3.2.2b,  the  expression  for  v*  required  in  Eq. 
3.2.4  is  obtain  by  differentiating  Eq .  3.2.2b.  This  yields 


. 
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v?  =  v*  -  y  sin  9 


d0 

dz 


(3.5.1) 


To  evaluate  the  term  d9/dz,  differentiate  the  expression  for 
appearing  in  Eq.  3.2.1  to  obtain 


it  d0 

V  =  cos  0  — 
o  dz 


(3.5.2) 


from  which 


d© 

dz 


_T  I 


COS0 


(3.5.3a) 


The  physical  interpretation  of  contained  in  Fig.  3-lb  allows  cos  0 
to  be  approximated  as 


cos  9  =  /  1  -  (v') 

o 


f  \  2 


(3.5.3b) 


Substituting  Eqs.  3.5.3  into  Eq.  3.5.1,  and  using  Eq.  3.2.1  yields 


v 


i 

v0 


y  v’  vM 
o  o 


/  i  -  (v;> 


(3.5.4) 


In  contrast  to  the  linear  formulation  (Sect.  A.l),  both  non¬ 
linear  terms  in  Eq.  3.2.4  are  now  retained.  Substituting  Eqs.  3.5.4 
and  3.2.3a  into  Eq .  3.2.4  yields  the  strain-displacement  equation  for 
the  nonlinear  formulation  as 


■  »;♦* 


K>2  +  ^ 


-  y  v 


(v  ’  )  2 

1  +  u*  +  -7====== 

O  /  1-  (vl)2 


o'  — 1 


+  \  y2(v”)2 
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(3.5.5) 
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The  variation  in  strain  6e  can  be  obtained  from  Eq.  3.5.5,  noting  that 

z 

Uo  *  Vo*  an<^  Vo  may  vary  iridePendently.  The  result  is 


6e 
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(3.5.6) 


As  in  Appendix  A,  the  principle  of  virtual  work  may  be  written 


as 


6W 


G  6£  dV  -  <Q>  {6q} 
z  z 

* 

v 


o 


(3.5.7) 


Substituting  Eq.  3.5.6  into  Eq.  3.5.7  results  in  the  equilibrium 
equation. 
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then  Eq. 


in  which 


and 


Let  stress  resultants  be  defined  as 


n  = 


a  •  dA 
z 

A 


(3.5.9a) 


m  = 


O  •  y  •  dA 
z 

4 

A 


(3.5.9b) 


m*  = 


a  •  y2  •  dA 
z 

a 

A 


(3.5.9c) 


3.5.8  can  be  rewritten,  using  the  definitions  in  Eqs.  3.5.9,  as 
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(3.5.10a) 
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Equations  3.5.10  are  approximate  only  to  the  extent  intro¬ 
duced  by  assumption  4  of  Sect.  3.2. 

If  the  displacements,  u  and  v  ,  are  functions  of  a  discrete 

o  o 

set  of  N  displacement  coordinates,  Eq.  3.5.10a  may  be  written  as 


0 


in  which 


3vT 

o 

a2  3q. 


+  a 


3v" 

o 


3  3q. 

iJ 


dz 


(3.5.11) 


(3.5.12) 


and  for  which  i  has  a  range  of  N.  (In  this  and  the  following  equations 
N  =  4) .  Equations  3.5.11  are  nonlinear.  If  Eqs .  3.5.11  are  not  satis¬ 
fied,  corrections  to  the  q^  may  be  obtained  by  Newton-Raphson  iteration 
according  to  (Murray  et  at.  1979) 

3^. 

M>.  =  ^Aq.  =  -ip.  (3.5.13) 

1  9q.  3  1 


in  which  the  summation  convention  is  used  for  repeated  indices  and  j 
also  has  a  range  of  N. 

Substituting  Eq.  3.5.12  into  Eq.  3.5.13  results  in 
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(3.5.14) 
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Equation  3.5.14  is  the  basic  Newton-Raphson  equation  and  may  be 

adapted  to  geometric  and  material  nonlinearity. 

Once  the  expressions  for  u  and  v  are  defined  in  terms  of 

o  o 

the  q’s,  all  terms  in  Eq.  3.5.14  may  be  evaluated  and  numerical  inte¬ 
gration  with  respect  to  z  may  be  carried  out  to  form  an  incremental 
tangent  stiffness  matrix  [k^]  and  unbalanced  load  vector  {AQ}.  Equation 
3.5.14  may  then  be  written  symbolically  as 


[kT]  (Aq}  =  {AQ} 


(3.5.15) 


in  which  the  influence  coefficient  (k^)„  is  evaluated  as 


(Vij  = 
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dz 


(3.5.16) 


;  and  the  unbalanced  load  AQ^,  is  evaluated  as 
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r  3u’ 
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dv") 
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a,  .  +  a 

[  1  3q.  2 

3q.  +  a3 
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dz 


(3.5.17) 


All  equations  derived  above  are  applicable  to  inelastic  as 
well  as  elastic  behavior.  The  equations  are  specialized  for  elastic 
response  in  Appendix  B  resulting  in  a  tangent  stiffness  matrix  [k^] , 
given  by  Eq.  B.4.4  for  a  flexural  element,  or  Eq.  B.4.11  for  a  truss 
element.  The  unbalanced  load  vector  { AQ }  is  given  by  Eqs.  B.3.4  and 
B.3.11  for  flexural  and  truss  elements,  respectively. 


3.5.2  Finite  Element  Model 

For  computational  purposes  a  finite  element  model  is  developed 
to  obtain  a  set  of  algebraic  equations  using  a  beam  element  function. 

In  the  finite  element  approach  the  displacements  at  any  point  are 
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represented  by  interpolating  functions  denoted  as  {<{>},  and  the  general¬ 
ized  displacements  {u}  and  {v}  at  the  nodes.  Thus 

u  =  «J»  {U}  (3.5.18a) 

o 

and 

vq  =  «I»  (V>  (3.5.18b) 

in  which  {u}  and  {v}  are  defined  in  Fig.  3-7  for  a  flexural  element  and, 
in  Fig.  3-8  for  a  truss  element. 

The  formulation  developed  in  this  study  uses  two  types  of 
interpolating  functions  for  two  types  of  element,  namely  a  flexural 
element  and  a  truss  element.  The  functions  selected  for  the  flexural 
element  are  cubic  polynomials  that  can  be  expressed  as  linear  combina¬ 
tion  of  the  four  base  functions  illustrated  in  Fig.  A-2  and  given  by 
Eq .  A. 2. 3b. 

For  a  truss  element,  linear  functions  are  used  that  may  be 
expressed  as  combination  of  the  two  base  functions  illustrated  in 
Fig.  A-l  and  given  in  Eq.  A. 2. 3a. 

Equations  3.5.18  are  used  in  Appendix  B  to  evaluate  the  finite 
element  equations,  Eq.  3.5.15. 

3.5.3  Assembly  of  Finite  Element  Equations 

Equation  3.5.15  is  the  incremental  equilibrium  equation,  for 
an  element.  In  this  equation  the  element  stiffness  matrix  has  been 
evaluated  with  respect  to  nodal  displacements  {q},  referenced  to  a 
local  coordinate  system  shown  in  Fig.  3-7  for  a  flexural  element  and 
Fig.  3-8  for  a  truss  element.  The  local  nodal  displacements  can  be 


written  as 
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{U}T  =  <uI,(u,Iz),  uJ,(u^z)> 
{V}T  =  <vI,  0j> 
{q1}1  =  <uI,(u,Iz),  v1,  9l> 


FIGURE  3-7  Local  Nodal  Displacements  for  Flexural  Element 
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{UJT  =  <u’,  uJ> 
(V}T=  <v',  vJ> 
{q'lT  =  <u‘,  v‘> 


FIGURE  3*8  Local  Nodal  Displacements  for  Truss  Element 
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<q>  =  <  <U>,  <V>  > 


(3.5.19a) 


or 


<q>  = 


<  1 

<  u  >  ^ 


J  ,9u.  J 
U  ’  (9l} 


v1  ,  01  ,  vJ  ,  0J  >  (3.5.19b) 


for  a  flexural  element;  and 


.  .  s  1  J  1 

<q>  =  <  u  ,  u  ,  v  ,  v  > 


(3.5.19c) 


for  a  truss  element;  and  I  and  J  are  nodal  numbers. 

In  frame  analysis  it  is  convenient  to  select  a  different  set 

of  reference  axis  for  the  global  system  of  nodal  displacements  {r  }  as 

E  G 

shown  in  Fig.  3-9.  The  global  nodal  displacements  can  be  written  as 


<r  >  =  <  r"^  .  r^  > 

EG  E  *  E  G 


(3.5.20a) 


in  which 


.3:.  .  TTI  „I  flI  /9Ux  I  ,9V.  I 

<rE>G  =  <  U  ’  V  >  6  >  3Z}  >  (3Y}  > 


(3.5.20b) 


The  element  displacements  with  respect  to  the  local  coordinate 
system  can  be  related  to  those  in  the  global  coordinate  system  by  the 
transformation 

{q}  =  [T]  {rE>G  (3.5.21) 

in  which  the  transformation  matrix  [T]  for  a  flexural  element  is  derived 
in  Appendix  C  and  given  in  Table  3-5.  Also,  the  transformation  matrix 
[T]  for  a  truss  element  is  given  in  Table  3-6. 
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=  <u\  v'.e'.tu.zj'.fv.v)1),. 


FIGURE  3-9  Global  Nodal  Displacements 
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Using  the  transformation  procedure  described  in  Sect.  3.3.2 
results  in 

[kT]G{arE}G  =  {iRE}G  (3-5'22) 

The  element  stiffness  matrices  can  now  be  assembled  to  form 

[Kt]  {Ar}  =  {AR}  (3.5.23) 

in  which  [K^]  is  the  structural  tangent  stiffness  matrix  assembled  for 
the  entire  structure,  {Ar}  is  the  assembled  vector  of  incremental  nodal 
displacements,  and  {Ar}  is  the  assembled  vector  of  incremental  nodal 
forces,  called  the  unbalanced  forces. 

Once  Eq.  3.5.23  is  assembled  the  Newton-Raphson  method  can 
be  used  to  solve  for  the  load-deformation  characteristics  of  the 
structure . 

3.5.4  Newton-Raphson  Solution  Procedure 

The  Newton-Raphson  method  has  proven  to  be  one  of  the  most 
useful  solution  techniques  available  for  nonlinear  analysis  (Rajasekaran 
and  Murray  1973) .  Many  investigators  have  adapted  the  method  quite 
successfully.  Walker  and  Hall  (1968)  used  it  to  study  large  deflec¬ 
tions  of  beams  while  Brebbia  and  Connor  (1969)  used  it  to  study  the 
geometrically  nonlinear  behavior  of  arbitrary  shells. 

The  Newton-Raphson  solution  technique  can  be  best  described 
with  reference  to  Fig.  3-10  and  the  following  algorithm. 

1  -  For  any  approximate  {r}^,  the  stiffness  matrix  [K^]  and  the  vector 

of  unbalanced  forces  {Ar}  are  evaluated. 

n 


‘ 
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FIGURE  3-10  Newton  -  Raphson  Procedure 
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2  -  The  increment  in  displacements  {Ar}  necessary  to  equilibrate  the 

unbalanced  forces,  evaluated  in  step  1,  are  evaluated  using 
Eq.  3.5.23  as 

{Ar}  =  [K  ] _ 1  {AR}  (3.5.24) 

n  T  n  n 

3  -  The  displacements  are  updated.  That  is 

{r}  =  {r}  +  {Ar}  (3.5.25) 

n+1  n  n 

4  -  Steps  1  to  3  are  repeated  until  the  vector  of  unbalanced  forces 

{AR}  is  arbitrarily  small.  At  this  stage  a  point  on  the  load- 
displacement  diagram  is  obtained. 

5  -  A  new  load  level  is  applied  and  steps  1  to  4  are  repeated  to  obtain 

the  load-deformation  characteristics  of  the  structure. 

3.5.5  Summary 

A  large  deformation  theory  for  plane  frame  problems  has  been 
formulated  in  this  section.  A  finite  element  method  has  been  developed 
and  a  Newton-Raphson  solution  procedure,  to  solve  for  the  complete 
load-deformation  curve,  is  explained. 

3.5.6  Computer  Program 

To  show  the  capability  of  the  theory,  formulated  above,  a 
computer  program,  called  Nonlinear  Elastic  Analysis  (NONELA) ,  is  developed 
to  solve  the  set  of  equations  presented  in  this  section  for  the  case  of 
linearly  elastic  material.  The  method  of  Newton-Raphson  is  used;  the 
stiffness  matrix  [K^]  being  evaluated  at  each  step  of  the  procedure. 

Some  applications  are  presented  in  the  next  section  to  demonstrate  the 
capability  of  the  theory  and  the  efficiency  of  the  computer  program. 
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3 . 6  Nonlinear  Elastic  Solutions 

The  examples  presented  in  Sect.  3.4  have  dealt  with  elastic 
buckling  or  bifurcation  problems.  When  large  deformation  theory, 
presented  in  Sect.  3.5,  is  considered  the  problem  becomes  nonlinear. 
For  this  type  of  problem  the  load-deformation  response  is  obtained  by 
applying  the  load  incrementally  and  iterating  to  find  the  equilibrium 
configuration  for  each  loading  condition  as  described  in  Sect.  3.5.4. 

3.6.1  Examples  of  Beam-Columns 
NL-1  'The  Elastica'  Problem 

This  example  considers  the  elastic  post-buckling  behavior  of 
a  simply  supported  beam-column.  A  small  transverse  load  was  applied 
to  act  as  an  initial  imperfection.  Results  are  shown  in  Fig.  3-11. 
Agreement  with  the  'Elastica'  solution,  discussed  in  Sect.  2.3.1, 
(Timoshenko  and  Gere  1961)  is  very  good.  The  number  of  elements 
used  was  4  elements.  The  average  number  of  iterations  to  get  from  one 
load  level  to  the  next  with  a  relative  error  of  1  x  10”4  was  4,  with 
an  average  cpu  time  per  iteration  of  0.1  sec. 

NL-2  Cantilever  Under  Pure  Moment 

In  this  example  a  cantilever  beam  acted  upon  by  a  moment  at 
its  free  end  is  analyzed.  A  very  good  agreement  with  the  solution  of 
Epstein  and  Murray  (1976) ,  was  obtained  with  four  elements  for  defor¬ 
mations  up  to  the  semicircumference.  Then  the  number  of  elements  was 
increased  to  six  and  deformations  pursued  up  to  the  complete  circum¬ 
ference.  Fig.  3-12  shows  some  of  the  exact  deformed  configuration 
with  the  corresponding  positions  of  the  nodes  as  obtained  from  the 
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FIGURE  3-1 1  Solution  for  the  ‘Elastica’ 
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FIGURE  3-1 2  Solution  for  Cantilever  Under  Pure  Moment 
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computer  program.  The  number  of  iterations  averaged  4  iterations  per 
each  load  increment  with  an  average  cpu  time  of  0.1  sec . /iteration. 

NL-3  Cantilever  Beam  With  One  Concentrated  Load 

The  results  summarized  in  Fig.  3-13  are  for  a  cantilever  beam 
loaded,  with  a  concentrated  load,  at  its  tip.  The  solid  line  represents 
the  exact  solution,  given  by  Bisshop  and  Drucker  (1945).  Two  solutions 
for  the  present  formulation  are  introduced.  One  using  2  elements  and 
the  other  using  4  elements.  It  is  seen  that  the  solution  obtained 
using  4  elements  is  in  excellent  agreement  with  the  exact  solution. 

NL-4  Cantilever  Beam  With  Two  Concentrated  Loads 

A  comparative  study  by  Ebner  and  Ucciferro  (1972)  gives 
numerical  results  obtained  by  several  different  finite  element  techni¬ 
ques  for  the  cantilever  beam  shown  in  Fig.  3-14.  The  results  reported 
are  reproduced  in  Table  3-7  along  with  results  obtained  by  the  present 
formulation.  It  is  seen  that  the  present  solution  with  two  elements 
is  as  accurate  as  the  solution  by  various  other  versions  of  the  finite 
element  approach. 

3.6.2  Examples  of  Frames 
NL-5  One  Story  Frame 

The  results  obtained  for  a  simple  one  story  bent  are  shown  in 
Fig.  3-15.  The  results  are  compared  with  an  approximate  solution  pre¬ 
sented  by  Oran  and  Kassimali  (1976).  It  is  important  to  note  that  the 
formulation  developed  as  part  of  the  present  investigation  did  function 
properly,  with  no  significant  convergence  difficulties,  up  to  even 
beyond  the  critical  load  of  the  limiting  case  of  the  perfect  bent.  It 
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FIGURE  3-1 3  Load  Deflection  Curves  for  Cantilever  Beam 
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FIGURE  3-14  Cantilever  Beam  with  Two  Lateral  Loads 
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may  be  of  interest  to  note,  in  this  connection,  that  an  alternate 
iterative  technique,  suggested  by  James  et  at,  (1974) ,  failed  to  con¬ 
verge  at  substantially  lower  load  levels.  A  lack  of  convergence  in 
the  iterative  process  was  viewed  by  those  authors  as  an  indication  of 
instability.  In  the  light  of  data  presented  herein,  such  a  conclusion 
would  not  necessarily  be  justified. 

The  solution  obtained  averaged  about  4  iterations  per  each 
load  increment  with  average  cpu  time  per  iteration  of  0.1  sec. 

NL-6  Multistory  Frames 

The  two  frames  shown  in  Fig.  3-5  and  Fig.  3-6  are  analyzed 
here.  The  results  for  the  two  story-two  bay  frame  are  shown  in  Fig. 
3-16.  In  this  example  the  average  number  of  iterations  to  get  from 
one  load  level  to  the  next  was  5  with  an  average  cpu  time  per  iteration 
of  0.2  second.  The  results  obtained  for  the  multistory  multibay  frame 
are  illustrated  in  Fig.  3-17.  In  this  case  the  number  of  iterations 
averaged  6  iterations  per  load  increment  with  an  average  cpu  time  of 
1.2  seconds  per  iteration. 

Unfortunately  there  is  no  exact,  or  approximate,  solution 
available  for  these  frames  to  compare  with.  However,  according  to  the 
results  obtained  from  the  previous  examples,  the  author  strongly  feels 
that  the  present  formulation  provides  good  results  for  these  frames. 

3.6.3  Comments  on  Results 

Results  obtained  in  this  section  compared  favorably  to  the 
exact  and  approximate  solutions  available.  The  examples  demonstrate 
that  the  computer  program  developed  is  capable  of  analyzing  a  wide 
variety  of  nonlinear  frame  problems  with  very  good  results. 
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FIGURE  3-1 6  Load  Deformation  Curves  for  Two  Story-Two  Bay  Frame. 
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FIGURE  3-1 7  Load-Deformation  Curve  for  Multistory-Multibay  Frame 


CHAPTER  IV 


BEHAVIOR  OF  ELASTIC  -  PLASTIC  FRAMES 

The  determination  of  the  complete  load-deformation  response 
for  plane  frames,  loaded  to  the  inelastic  range,  is  possible  by  a 
computerized  step-by-step  elastic-plastic  analysis,  described  in 
Sect.  2.4.5.  In  this  chapter,  features  of  elastic-plastic  analysis  of 
frames  are  presented.  Two  methods  have  been  adapted,  one  is  an  incre¬ 
mental  solution  while  the  other  is  an  iterative  technique.  For  each 
technique  first  and  second  order  analyses  are  presented.  A  comparison 
between  the  methods  is  furnished  and  results  are  compared  to  published 
reports . 

4 . 1  Basic  Assumptions  and  Limitations 

In  addition  to  the  first  three  assumptions  presented  in  Sect. 
3.2,  the  following  assumptions  are  also  considered. 

1  -  The  material  is  assumed  to  be  linearly  elastic-perf ectly  plastic 

(i.e.  strain  hardening  is  not  considered). 

2  -  Residual  stresses  are  ignored. 

3  -  Plastic  deformations  are  concentrated  in  generalized  hinges. 

4  -  Plastic  moment  capacity,  M  ,  in  the  presence  of  axial  load,  P, 

at  a  particular  section  is  reduced  in  accordance  with  the  inter¬ 
action  yield  curve  shown  in  Fig.  4-1. 
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M  =  Mp 

M  =1.18  (1-P/Py)  Mp 


For  0<P/Py<.15 
For  .15<P/Py<1.0 


FIGURE  4-1 


Interaction  Curve 
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4 . 2  The  Incremental  Method  of  Analysis 

Equation  3.3.11  is  the  equilibrium  equation,  for  a  structure, 
assembled  in  a  global  coordinate  system.  In  that  equation  the  geometric 
stiffness  matrix  [K  ]  depends  on  the  axial  forces  in  the  individual 
members  which  in  turn  depend  on  the  level  of  loading. 

Since  the  geometric  stiffness  matrix  [K  ]  varies  between  two 

Cj 

levels  of  loading,  Eq.  3.3.11  can  be  put  in  an  incremental  form  as 

[K]  {Ar}  +  [Kg]  {Ar}  +  [AKg]  {r}  =  {AR}  (4.2.1) 

Equation  4.2.1  can  be  written,  symbolically,  as 

[K^  {Ar}  =  {AR}  -  [AKg]  {r}  (4.2.2) 

which  is  an  incremental  equilibrium  equation  that  can  be  used  in  the 
step-by-step  elastic-plastic  analysis  to  obtain  the  complete  load- 
deformation  characteristics  of  a  frame. 

4.2.1  Incremental  First  Order  Analysis 

Consider  a  frame  subjected  to  two  sets  of  loads,  one  of  which 
is  constant,  {Rq}»  while  the  other,  {AR^},  varies  proportionately.  In 
conducting  a  first  order  analysis  the  geometric  stiffness  is  not  con¬ 
sidered  and  Eq.  4.2.2  becomes 

[Kt]  {Ar}  =  {AR}  (4.2.3) 


in  which  [K^,]  is  the  tangent  stiffness  matrix  not  including  the  geometric 
nonlinearity. 
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The  step-by-step,  elastic-plastic  first  order  analysis,  to 
obtain  the  load-deformation  response  of  a  frame,  can  be  best  described 
using  the  following  procedure. 

1  -  An  elastic  analysis  is  performed  based  on  the  original  elastic 

stiffness  matrix  of  the  structure.  The  vectors  of  nodal  displace¬ 
ments  {r  }  and  {Ar  }  as  well  as  the  vectors  of  member  forces  {Q  } 
c  v  c 

and  {AQv},  due  to  the  two  sets  of  loads  considered  above,  are  obtained. 

2  -  The  moment  capacity  of  each  critical  section  is  calculated  as 


(M  )  .  =  (M  )  .  -  (M  )  . 

pv  1  pc  1  Cl 


(4.2.4) 


3  - 


in  which  (M  ) ^  is  the  moment  capacity  that  would  resist  the  moment 

at  the  ith  section  due  to  the  variable  set  of  loads  {AR  } ,  (M  ) .  is 
—  v  pc  1 

the  plastic  moment  capacity  of  the  ith  section  and,  (M  ) ^  is  the 

moment  at  the  ith  section  due  to  the  constant  set  of  loads  {R  }. 

c 

The  load  factors  X.  are  calculated,  for  each  section, as 

l 


X 


(M  )  . 

pv  1 

i  (AM  )  . 

v  i 


(4.2.5) 


in  which  (AM  ) .  is  the  moment  at  the  ith  section,  due  to  the  variable 
vi  — 

set  of  loads,  {AR^}.  The  section  having  the  smallest  load  factor  will 
be  the  first  section  to  be  plastified.  This  load  factor  then  is 
called  X^  which  will  cause  the  first  hinge  to  be  formed. 

4  -  The  applied  loads,  nodal  displacements,  and  member  forces  are  up¬ 
dated,  using  the  factor  X^,  as 

{R.}  =  (R  }  +  X.  {AR  }  (4.2.6a) 

1  civ 

{rl}  =  <rc}  +  X1  {Arv} 


(4.2.6b) 
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and 


{Q, }  =  (Q  }  +  A  {AQ  }  (4.2.6c) 

1  c  1  v 

in  which  {R^},  {r^},  and  {Q^}  are  vectors  of  applied  loads,  nodal 
displacements,  and  member  forces  respectively  at  the  formation  of 
the  first  plastic  hinge. 

5  -  Now  the  first  plastic  hinge  has  been  formed.  A  real  hinge  is 

placed  at  this  section  and  the  incremental  stiffness  matrix  of 
the  structure  is  modified  by  changing  the  particular  member  stiff¬ 
ness  matrix  where  the  plastic  hinge  has  appeared. 

6  -  An  increment  in  the  variable  set  of  loads  {AR^}  is  applied  to 

the  structure  and  a  new  elastic  analysis  is  performed  on  the 
basis  of  the  modified  stiffness  matrix.  The  member  forces,  nodal 
displacements  and  plastic  deformations  due  to  the  load  increment 
can  then  be  evaluated. 

7  -  The  moment  capacity  remaining  in  each  section  may  then  be 

calculated  as 


(M  )  .  - 
pc  1 


<Vi 


(4.2.7) 


in  which  (M^K  is  the  moment  at  the  ith  section  developed  at  the 
formation  of  the  first  plastic  hinge  and  can  be  obtained  from  the 
updated  vector  of  member  forces  { } ,  given  in  step  4. 

8  -  The  search  for  the  next  section  to  be  plastified  is  conducted, 
as  described  in  step  3,  and  the  applied  loads,  nodal  displace¬ 
ments  and  member  forces  are  updated  accordingly  as 


(r2>  =  (r1)  +  a2  (arv> 


(4.2.8a) 


Ill 


{r0}  =  (r  }  +  X  {Ar  }  (4.2.8b) 

2  1  2  v 

and 

{Q2>  =  +  X2  {AQ^}  (4.2.8c) 

in  which  X^  is  the  smallest  load  factor  obtained. 

The  above  procedure  is  continued  until  a  failure  mechanism 
is  formed.  The  failure  criterion,  as  far  as  the  incremental  analysis 
is  concerned,  is  the  singularity  of  the  global  tangent  stiffness  matrix 
of  the  structure.  When  a  mechanism  forms,  a  zero  element  appears  on  the 
main  diagonal  of  the  triangularized  structure  stiffness  matrix  and  the 
solution  is  completed. 

4.2.2  Incremental  Second  Order  Analysis 

In  the  case  of  a  second  order  analysis,  geometric  nonlinearity 
is  considered  and  Eq.  4.2.2  is  used.  The  same  procedure  for  first  order 
analysis,  described  in  Sect.  4.2.1,  is  used  for  second  order  analysis 
with  few  changes. 

First,  a  new  geometric  stiffness  matrix  [K  ]  has  to  be  assembled 

Cj 

at  the  beginning  of  each  load  increment.  This  matrix  is  a  function  of 

the  axial  forces  in  the  members  up  to  this  stage  of  loading. 

The  second  change,  that  has  to  be  introduced,  is  to  modify  the 

vector  of  incremental  loads  {AR^}  to  account  for  the  quantity  [AK^Hr} 

in  Eq.  4.2.2.  The  incremental  geometric  stiffness  matrix  [AK  ]  is  a 

function  of  the  incremental  axial  forces  in  the  members  due  to  {AR  } 

v 

while  {r}  is  the  vector  of  total  nodal  displacements  obtained  from  a 
previous  step. 
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To  account  for  the  reduction  in  plastic  moment  capacity  due  to 
the  presence  of  axial  forces  in  the  members,  the  interaction  curve  shown 
in  Fig.  4-1  is  used.  Two  alternative  methods  are  suggested.  The  first 
method  is  to  estimate,  approximately,  the  final  axial  force  in  each 
member  and  reduce  the  moment  capacity  according  to  this  estimate.  The 
other  alternative  is  to  calculate  the  moment  capacity  based  on  the 
axial  force  developed  at  the  beginning  of  the  load  increment,  neglecting 
any  change  in  this  value  due  to  any  change  in  axial  force  once  a  hinge 
has  formed  at  this  particular  location.  Both  alternatives  suggested 
above  are  approximate  solutions  and  result  in  some  error.  The  second 
suggestion  is  adapted  in  the  analysis  presented  in  this  section.  A 
more  precise  method  is  presented  in  the  following  section. 

Once  the  above  modifications  are  introduced,  the  procedure 
described  in  Sect.  4.2.1  is  followed  until  a  mechanism  is  formed. 

When  a  mechanism  forms  the  determinant  of  the  decomposed  structure  stiff¬ 
ness  matrix  will  have  a  negative  sign.  At  this  stage  a  negative  incre¬ 
ment  in  load  is  applied  to  the  structure  to  obtain  the  descending  line 
of  load-deformation  characteristic,  shown  in  Fig.  1-1  (Curve  F) . 

4 . 3  The  Iterative  Method  of  Analysis 

In  this  section  the  equilibrium  equation,  Eq.  3.3.11,  is 
satisfied  at  each  equilibrium  position,  (i.e.  at  the  formation  of 
each  hinge) . 

4.3.1  First  Order  Analysis 

For  first  order  analysis,  the  geometric  nonlinearity  that  appears 
in  Eq.  3.3.12  is  not  considered.  In  such  a  case,  Eq.  3.3.11  can  be  written 


as 
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[K]  {r}  =  {R} 


(4.3.1) 


in  which  [K]  is  the  elastic  stiffness  matrix  assembled  for  the  entire 


structure . 


The  procedure  presented  in  Sect.  4.2.1  can  be  used  here,  with 
Eq.  4.3.1,  up  to  and  including  step  5.  However,  the  new  stiffness  matrix 
in  this  case,  obtained  from  step  5,  is  a  tangent  stiffness  matrix. 

The  loads  that  produced  the  first  plastic  hinge  are  then  applied 
to  the  modified  structure  together  with  an  applied  moment  equal  to  the 
plastic  moment  capacity,  M  ,  at  the  location  of  the  hinge.  The  first 
set  of  loads  (R  }  increases  until  the  formation  of  the  second  hinge 
while  the  applied  moment  which  can  be  assembled  into  the  load  vector 
(Rc)  is  kept  constant.  It  is  important  to  note  that  because  the  stiff¬ 
ness  matrix  is  a  tangent  stiffness, the  set  of  displacements  {r^}  and 
{r^}  are  fictitious  displacements.  However,  the  difference  is  the  real 
displacement  vector  {r}  defined  in  Fig.  4-2. 

The  maximum  moment  capacity  left  in  any  section  after  the  first 
hinge  was  formed  can  be  obtained  from  Eq.  4.2.4  and  the  search  for  the 
second  section  to  be  plastified  is  based  on  the  smallest  factor  X  defined 


by  Eq.  4.2.5. 

Once  X  is  determined,  the  applied  loads,  nodal  displacements 
and  member  forces  are  updated.  Also,  the  constant  set  of  forces  {R^} 
is  modified  to  include  the  plastic  moment  capacity  as  an  applied  moment 
at  the  second  hinge  preparing  for  the  next  equilibrium  position. 

The  above  procedure  is  continued  until  a  mechanism  forms. 


* 
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FIGURE  4-2  Load-Deformation  Response 
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4.3.2  Iterative  Second  Order  Analysis 

The  first  order  analysis  as  described  in  Sect.  4.3.1  neglects 
the  effects  of  axial  forces  in  the  members  as  they  influence  the  geo¬ 
metric  nonlinearity  of  the  structure. 

To  account  for  the  geometric  nonlinearity,  the  geometric 
stiffness  matrix  is  included  in  the  analysis  (Eq.  3.3.12).  In  this 
case,  a  new  geometric  stiffness  matrix  is  formed  at  each  load  level 
based  on  the  previous  load  level.  This  geometric  stiffness  matrix  is 
to  be  corrected  through  an  iterative  process  up  to  the  new  level  of 
loading  to  be  consistent  with  the  updated  applied  loads  {R^}. 

The  reduction  in  moment  capacity  due  to  the  presence  of 
axial  force  P  at  any  section  is  considered  based  on  the  interaction 
curve  shown  in  Fig.  4-1.  The  ratio  P/P^  at  any  section  is  checked  to 
ensure  that  the  moment  capacity  of  the  section  does  not  exceed  the  plastic 
moment  capacity,  M  c< 

A  new  moment  capacity  M  should  be  calculated  at  each  itera¬ 
tion  to  account  for  the  change  in  the  axial  forces.  Consequently,  the 

constant  set  of  forces  {R  },  described  in  Sect.  4.3.1,  must  be  modified 

c 

in  accordance  with  the  new  moment  capacity,  M  ,  at  each  iteration. 

The  same  procedure,  described  in  Sect.  4.3.1,  is  then  used 
for  an  iterative  second  order  analysis. 

4 . 4  Comparison  Between  Incremental  and  Iterative  Approaches 

In  the  incremental  method  the  geometric  stiffness  matrix  is 
evaluated  at  the  beginning  of  the  load  increment  as  a  function  of  the 
total  axial  forces  in  the  members  up  to  this  step.  In  the  iterative 
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approach,  however,  the  geometric  stiffness  matrix  is  re-evaluated  with 
each  iteration  as  a  function  of  the  total  member  axial  forces.  This 
results  in  a  more  accurate  geometric  stiffness  matrix  and »consequently, a 
more  accurate  solution. 

As  discussed  in  Sect.  4 . 2 . 2 ,  the  incremental  procedure  does 
not  allow  any  change  in  the  plastic  moment  capacity,  M^c , of  a  section 
due  to  any  change  in  the  axial  force,  once  a  hinge  has  formed  at  this 
particular  section.  Two  alternative  solutions  were  suggested  in  Sect. 
4.2.2  to  overcome  this  problem.  However,  both  alternatives  are  approxi¬ 
mate  and  result  in  some  error  when  evaluating  the  load-deformation 
response  of  the  structure.  In  the  iterative  procedure , the  plastic 
moment  capacity  of  each  section  is  re-evaluated  with  each  iteration  to 
account  for  any  change  in  the  axial  forces.  This  results  in  a  more 
accurate  analysis. 

The  effect  of  the  approximations  made  in  the  incremental  method 
will  be  demonstrated  by  some  example  applications  in  Sect.  4.6. 


4 . 5  Computer  Programs 

Two  computer  programs  have  been  developed  to  solve  for  the  load- 
deformation  characteristics  of  frames  by  plastic  analysis.  The  first 
program  performs  an  incremental  plastic  analysis  of  frames,  based  on 
the  procedure  described  in  Sect.  4.2,  and  is  called  INPLAF.  The  second 
program  performs  plastic  analysis  of  frames  using  the  iterative  technique, 
described  in  Sect.  4.3,  and  is  called  PLAFIT.  Each  program  has  the 
option  to  perform  first  or  second  order  analysis. 


« 
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4 . 6  Elastic-Plastic  Solutions 

Three  example  frames  are  checked  in  this  section  using  the 
incremental  and  iterative  procedures  presented  in  Sects.  4.2  and  4.3, 
respectively.  The  load-deformation  characteristics  for  each  frame 
are  obtained  using  the  different  analyses  and  the  results  are  compared 
with  published  results.  For  each  analysis  the  cpu  computer  time  is 
given  as  a  measure  of  the  cost  of  running  the  program. 

4.6.1  Sample  Solutions 
EP-1  Four  Story  Frame 

The  frame  shown  in  Fig.  4-3  was  analyzed  by  Korn  and  Galambos 
(1968) .  In  the  analysis  the  loads  were  increased  proportionately. 
Figure  4-4  illustrates  the  load-deformation  curves  obtained  by  the 
different  analyses.  Curves  (1),  (2),  (3),  and  (4)  represent  iterative 
first  order,  iterative  second  order,  incremental  first  order,  and 
incremental  second  order  analysis,  respectively.  The  results  obtained 
by  Korn  and  Galambos  (1968)  are  identical  to  curves  (1)  and  (2). 

The  cpu  computer  time  for  first  order  analysis  is  about  1.7 
seconds  using  the  incremental  approach  and  5  seconds  for  the  iterative 
approach.  For  second  order  analysis  the  incremental  technique  used 
about  1.3  sec.  compared  to  20  sec.  for  the  iterative  technique. 

As  can  be  observed  from  Fig.  4-4  the  incremental  second  order 
analysis  (curve  4)  is  considerably  different  from  the  iterative  second 
order  analysis  (curve  2) .  This  difference  basically  occurs  during  the 
formation  of  the  first  hinge  as  the  geometric  stiffness  matrix  in  the 
incremental  analysis  is  formulated  at  the  beginning  of  the  load  incre¬ 
ment  which  in  this  case  results  in  a  null  geometric  stiffness  matrix. 


■ 
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All  Members  are  Ficticious 

A  =11.75  in2 
I  =144  in4 
Mp=  160  ft-  K 
E  =  30,000  Ksi 
<Ty  =  36  Ksi 


FIGURE  4-3  Frame  EP-1 
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FIGURE  4-4  Load-Deformation  Relationship  for  Frame  EP-1 
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At  this  stage  it  was  thought  that  better  results  could  be 
obtained  if  an  iterative  procedure  were  applied  only  up  to  the  develop¬ 
ment  of  the  first  hinge  followed  by  an  incremental  procedure  up  to  the 
collapse  mechanism.  This  procedure,  which  may  be  called  iterative- 
incremental  second  order  analysis,  results  in  a  load-deformation 
relationship  given  by  curve  (5)  in  Fig.  4-4  and  compares  favorably  to 
curve  (2).  The  cpu  computer  time  used  to  obtain  curve  (5)  is  about  5  sec. 
which  is  considerably  less  than  that  used  to  obtain  curve  (2) . 

EP-2  Three  Story  Frame 

McNamee  and  Lu  (1972)  performed  a  second  order  analysis  on  the 
frame  shown  in  Fig.  4-5  increasing  the  loads  proportionally.  The  load- 
deflection  curve  they  obtained  is  curve  (2)  in  Fig.  4-6.  Curves  (1) 
and  (3)  represent  the  iterative  and  incremental  second  order  analysis, 
respectively,  presented  in  this  chapter.  While  the  cpu  computer  time 
for  the  iterative  technique  is  about  4  seconds,  it  is  1.2  sec.  for  the 
incremental  technique.  It  is  important  to  note  that  an  iterative- 
incremental  analysis,  suggested  in  the  previous  example,  results  in  a 
load-deflection  curve  identical  to  curve  (1)  and  uses  only  about  2  sec. 
cpu  computer  time. 

EP-3  One  Story  Three  Bay  Frame 

In  the  frame  shown  in  Fig.  4-7  the  vertical  loads,  P,  are  kept 
constant  while  increasing  the  horizontal  force,  H.  The  results  of  first 
and  second  order  analyses  presented  by  Galambos  and  Lay  (1965)  are  given 
as  curves  (3)  and  (6) ,  respectively,  in  Fig.  4-8.  The  present  study 
gives  an  iterative  first  order  analysis  (curve  1)  and  an  incremental 
first  order  analysis  (curve  2)  which  are  almost  identical.  The  cpu  time 
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FIGURE  4-5  Frame  EP-2 
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FIGURE  4-6  Load-Deflection  Curve  for  Frame  EP-2 
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FIGURE  4-8  Load-Deformation  Characteristics  for  Frame  EP-3 
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is  about  0.6  sec.  for  the  iterative  technique  and  0.45  sec.  for  the 
incremental  one. 

The  second  order  analyses  result  in  curves  (4)  and  (5)  represent¬ 
ing  the  iterative  and  incremental  technique,  respectively.  The  cpu 
computer  time  is  about  2.4  sec.  and  .90  sec.  for  the  iterative  and 
incremental  technique,  respectively.  When  an  iterative- incremental 
procedure  is  performed  curve  (7)  is  obtained  with  a  cpu  time  of  1.4  sec. 

4.6.2  Comments  on  Results 

The  examples  presented  above  indicate  that  the  iterative 
method  of  analysis,  presented  in  Sect.  4.3,  exhibits  good  results. 

However,  an  iterative  method  up  to  the  development  of  the  first  hinge 
followed  by  an  incremental  procedure  until  collapse  gives  approximately 
the  same  results  and  does  not  cost  as  much  as  the  iterative  method. 


CHAPTER  V 


NONLINEAR  INELASTIC  FORMULATION  AND  SOLUTIONS 


Nonlinear  formulations  relating  kinematic  and  mechanical 
variables  may  fall  into  one  of  the  three  categories:  geometric  non¬ 
linearity,  material  nonlinearity,  or  combined  geometric  and  material 
nonlinearities.  The  formulation  presented  in  Chapter  3  falls  into  the 
first  category  while  this  chapter  is  devoted  to  problems  falling  into 
the  third  one. 

In  this  chapter  an  inelastic  stability  analysis  of  multistory 
frames  is  presented  based  on  a  stiffness  formulation  which  accounts 
for  geometric  as  well  as  material  nonlinearity.  The  formulation  is  an 
extension  of  that  presented  in  Sect.  3.5  to  include  material  nonlinearity. 
The  formulation  results  in  finite  element  equations  and  the  Newton- 
Raphson  method  is  then  used  to  solve  for  the  overall  load-deformation 
characteristics  of  the  structure. 

In  the  analysis  presented  herein,  the  effect  of  axial  load  on 
the  stiffness  and  strength  of  the  individual  members  is  considered. 

Partial  plastif ication  of  sections  is  taken  into  account.  The  influence 
of  residual  stresses  as  well  as  strain  hardening  of  the  material  is  in¬ 
cluded  in  the  analysis.  The  formulation  also  permits  consideration  of 
extended  regions  of  plastic  deformations  rather  than  discrete  hinges 
in  beams  and  beam-columns. 

Also  in  this  chapter,  the  formulation  is  checked  by  comparing 
computed  results  with  some  of  the  available  experimental  test  results. 
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5 . 1  Basic  Formulation  for  Inelastic  Behavior 

Equation  3.5.14  is  the  basic  incremental  equilibrium  equation, 
derived  in  Sect.  3.5.1,  and  is  applicable  to  elastic  and  inelastic 
response.  Equation  3.5.14  can  be  written  in  matrix  form,  for  inelastic 
response,  as 

[k*]  {Aq}  =  {AQ1}  (5.1.1) 

in  which  [k^]  is  the  inelastic  element  tangent  stiffness  matrix  and  is 
a  function  of  geometric  and  material  nonlinearities  and  {AQ1}  is  the 
vector  of  incremental  forces,  evaluated  for  an  inelastic  element. 

The  difference  between  the  equilibrium  equation  for  an  elastic 
element,  Eq.  3.5.15,  and  that  for  an  inelastic  element,  Eq.  5.1.1,  basic¬ 
ally  arises  in  the  evaluation  of  the  tangent  stiffness  matrix  and  the 
incremental  load  vector. 

5.1.1  Main  Assumptions 

In  deriving  the  tangent  stiffness  for  a  partially  plastic 
section,  it  is  assumed  that  no  strain  reversal  occurs.  If  one  considers 
a  section  which  has  already  yielded  and  an  infinitesimal  increment  in 
bending  moment  is  applied,  the  resistance  to  this  moment  is  a  measure 
of  incremental  bending  stiffness.  The  yielded  zones  increase  somewhat, 
while  some  strain  reversal  may  take  place  in  the  already  yielded  zones. 
These  effects  are  neglected  in  this  study  and  the  additional  moment  is 
assumed  to  be  resisted  by  all  parts  of  the  cross  section  with  the  stress 
increment  determined  from  the  tangent  modulus  of  the  material  at  the 
current  value  of  strain.  The  section  properties  relating  increments 
in  stress  resultants  to  increments  in  displacements  may  then  be  computed 
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by  a  transformed  area  concept  as  described  in  Sect.  D.l.  This  concept 
is  applicable  to  thin  wall  open  beam  sections  since  the  stress  is 
essentially  uniaxial  and  the  effect  of  shear  stress  on  the  principal 
stresses  may  be  neglected. 

5.1.2  Evaluation  of  Incremental  Forces  {AQ^} 

To  evaluate  the  vector  of  incremental  forces  {AQ^},  in  Eq.  5.1.1, 
it  is  necessary  to  evaluate  the  stress  resultants  n,  m,  and  m*,  defined 
by  Eqs.  3.5.9,  at  any  section  prior  to  the  load  increment.  Equations 
B.3.4  and  B.3.11  give  the  vectors  of  incremental  forces,  for  a  flexural 
element  and  truss  element,  respectively,  assuming  an  elastic  response. 

The  same  equations  are  applicable  to  inelastic  response  provided  that 
the  stress  resultants  are  evaluated  by  direct  integration  of  stresses 
as  detailed  in  Appendix  D,  Eqs.  D.3.1  through  D.3.3. 


5.1.3  Evaluation  of  Element  Tangent  Stiffness  Matrix 

The  element  tangent  stiffness  matrix  for  an  elastic  element 
is  given  by  Eq.  B.4.4  for  a  flexural  element,  and  by  Eq.  B.4.11  for  a 
truss  element.  These  equations  are  applied  to  inelastic  response  with 
two  modifications. 

First  the  stress  resultants  n,  m  and  m*  are  numerically 
evaluated  by  direct  integration  of  stresses  (Eqs.  D.3.1b,  D.3.2b, 
and  D. 3 . 3b) . 

The  second  modification  that  should  be  introduced  for  inelastic 


analysis  is  in  the  evaluation  of  the  incremental  stress  resultants 
3n 


3m  ,  3m* 

,  and 


These  increments  are  evaluated  in  detail  in 


3qj  *  3qj  *  3qj 

Sect.  B . 2  for  an  elastic  element.  For  an  inelastic  element  detailed 


evaluation  of  the  incremental  stress  resultants  is  given  in  Appendix  E 
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and  the  incremental  stress  resultant  vectors  are  given  in  Tables  E-l 
and  E-2  for  a  flexural  element  and  a  truss  element,  respectively. 

Once  the  stress  resultants  are  evaluated  and  the  incremental 
stress  resultant  vectors  are  formed  the  same  formulation,  presented 
in  Sect.  B.4  to  derive  the  elastic  tangent  stiffness  matrix  can  be 
used  to  derive  the  inelastic  tangent  stiffness  matrix  for  an  inelastic 
element.  It  should  be  noted  that  the  finite  element  shape  functions 
selected  for  an  elastic  analysis  (Eqs.  A. 2. 3)  are  also  applicable 
to  inelastic  response  evaluation. 

5 . 2  Transformation  and  Assembly 

The  incremental  equilibrium  equation  (Eq.  5.1.1)  is  formulated 
in  the  y-z  coordinate  directions  which  are  the  reference  axes  of  the 
member.  The  local  incremental  displacements  {Aq}  may  be  related  to  the 
incremental  displacements  in  the  global  system  {Ar}  for  a  typical 
element  by  the  transformation  matrices  given  in  Tables  3-5  and  3-6  for 
a  flexural  element  and  a  truss  element,  respectively.  Assembly  of  the 
incremental  equilibrium  equation  then  proceeds  by  the  direct  stiffness 
method,  presented  in  Sect.  3.3.2,  to  yield 

[K*]  {ir}  =  {AR1}  (5.2.1) 

in  which  [K^]  is  the  inelastic  tangent  stiffness  matrix  for  the  entire 
structure  and  is  obtained  by  direct  addition,  {Ar}  is  the  assembled 
vector  of  incremental  displacements,  and  {AR'*'}  is  the  assembled  vector 
of  incremental  forces  (called  the  unbalanced  load  vector) . 

Equation  5.2.1  is  the  equilibrium  equation  for  the  structure 
and  becomes  the  basis  for  the  Newton-Raphson  solution. 
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5.3  Newton-Raphson  With  Gradient  Test 

The  Newton-Raphson  technique,  described  in  Sect.  3.5.4,  is 
used  in  this  chapter  to  trace  the  behavior  of  inelastic  structures. 

The  procedure  is  identical  to  that  of  Sect.  3.5.4  except  for  the 
detailed  evaluation  of  tangent  stiffness  matrix  and  the  vector  of 
unbalanced  forces.  The  computational  procedure  and  the  associated 
incremental  equilibrium  equation  (Eq.  5.2.1)  are  applicable  up  to  the 
peak  of  the  load-deflection  curve. 

To  be  able  to  trace  the  load-deflection  characteristics 
beyond  the  peak  of  the  curve,  a  special  treatment  must  be  adapted 
(Rajasekaran  and  Murray  1973).  The  descending  branch  of  the  load- 
deflection  curve  in  the  inelastic  region  is  characterized  by  a  negative 
definite  stiffness  matrix.  This  means  that  the  structure  can  only 
withstand  a  decrease  in  the  load. 

To  obtain  the  descending  branch  of  the  load-deflection  curve, 
the  determinant  of  the  tangent  stiffness  matrix  is  checked  for  sign 
change  to  determine  if  the  stiffness  matrix  is  still  positive  definite. 
If  this  sign  is  negative  the  load  is  decremented.  The  test  for  sign 
is  called  a  'gradient  test'  and  will  determine  exactly  where  to  start 
decreasing  the  applied  loads  in  order  to  maintain  equilibrium. 

5 . 4  Residual  Stresses  and  Strain  Hardening  of  the  Material 

Residual  stresses  are  caused  by  a  variety  of  factors  and  in 
some  cases  they  may  be  as  large  as  the  yield  stresses  (Huber  and  Beedle 
1954;  Tall  1964).  They  are  a  result  of  thermal  effects  and  plastic 
deformation  during  the  manufacturing  process.  Tall  (1964)  ,  Huber  and 
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Beedle  (1954) ,  and  Beedle  and  Tall  (1960)  provide  a  large  number  of 
illustrations  showing  the  measured  residual  stresses  caused  by  hot 
rolling,  welding,  or  flange  cutting  for  steel  cross  sections  of  various 
shapes  and  sizes. 

In  this  investigation  residual  stresses  are  considered. 

These  may  have  a  considerable  effect  on  the  behavior  of  the  structure, 
especially  during  the  transition  from  the  fully  elastic  to  the  partly 
yielded  condition.  Full  treatment  of  residual  stresses  is  discussed  in 
Sect.  D.l. 

In  most  of  the  elastic-plastic  formulations  developed  for 
frame  analysis  the  material  is  assumed  to  be  ideally  elastic-plastic  and 
the  influence  of  strain  hardening  on  the  ultimate  load  is  neglected.  In 
the  formulation  presented  in  this  chapter  an  attempt  is  made  to  include 
the  strain  hardening  of  the  material  in  the  analysis.  In  this  case  the 
idealized  stress-strain  relationship  shown  in  Fig.  D-2  is  assumed.  The 
treatment  is  detailed  in  Appendix  D  (Sect.  D.l). 

5 . 5  Plastif ication  of  the  End  Zones 

Because  of  the  distribution  of  the  Gauss  sampling  points, 
discussed  in  Sect.  B.3.1,  plastif ication  at  the  ends  of  members  does 
not  immediately  show  its  effect  on  the  element  tangent  stiffness  matrix 
when  treating  each  member  as  one  element.  For  this  reason  better 
results  may  be  obtained  by  subdividing  the  member  into  three  elements 
with  the  central  element  being  80%  of  the  length  of  the  member.  The 
result  is  that  a  number  of  Gauss  points  are  located  on  each  end  within 
the  outer  10%  of  member  length.  Therefore,  the  influence  of  inelastic 
strains  is  felt  shortly  after  they  are  initiated  at  the  member  ends. 
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An  element  tangent  stiffness  matrix  is  formed  for  each  element 
and  a  static  condensation  scheme,  developed  by  Elwi  and  Murray  (1977) 
is  used  to  eliminate  the  internal  degrees  of  freedom  and  to  then 
assemble  only  the  member  end  degrees  of  freedom  into  the  structure 
stiffness  matrix. 

5 . 6  Computer  Program 

A  computer  program,  called  Inelastic  Stability  Analysis  of 
Frames  (INSTAF) ,  has  been  developed  to  solve  for  the  load-deformation 
characteristics  of  an  inelastic  structure.  The  program  can  analyze 
braced  and  unbraced  frames  and  has  the  provision  of  including  residual 
stresses  and  strain  hardening  of  the  material  in  the  analysis. 

5. 7  Inelastic  Solutions 

To  demonstrate  the  accuracy  and  to  test  the  efficiency  of  the 
inelastic  formulation  presented  in  this  chapter,  some  of  the  available 
experimental  test  results  are  checked  using  the  computer  program.  The 
results  are  presented  in  this  section. 

5.7.1  Sample  Examples 
IN-1  Beam-Column  Test 

A  W4  x  13  beam-column  was  tested  by  Van  Kuren  and  Galambos 

(1964),  in  which  the  axial  load,  P,  was  kept  constant  while  the  applied 

moment,  M  ,  was  increased  until  failure  took  place.  Figure  5-1  shows 
o 

the  comparison  between  the  test  results  (curve  1)  and  the  theoretical 
moment-rotation  curve  obtained  using  the  present  formulation  (curve  2) . 
Residual  stress  data  were  not  given  by  Van  Kuren  and  Galambos  and  therefore 
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(1)  Test  Results 

(2)  Theoretical  Study 

(3)  Theoretical  Study  with  Residual 
Stresses  and  Lateral  Load 


FIGURE  5-1  Results  of  Example  IN-1 
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were  not  included  in  the  theoretical  study.  This  may  explain  the 
difference  in  ultimate  capacities  between  the  two  curves.  The  fact  that 
the  slopes  of  the  moment-rotation  curves  are  different  indicates  that 
there  was  geometric  imperfection  in  the  member  before  testing.  Using 
an  analysis  of  the  beam-column  where  a  small  lateral  load  (0.001  kips) 
was  applied  at  midspan,  the  geometric  imperfection  was  simulated.  An 
assumed  residual  stress  distribution  was  used,  with  a  maximum  compres¬ 
sive  value  of  0.3  The  analysis  gave  curve  3,  which  compares  favor¬ 

ably  to  the  test  results. 

Twelve  load  increments  were  used  to  draw  the  complete  moment- 
rotation  curve,  obtained  in  the  present  analysis,  with  an  average  of 
5  iterations  for  each  load  increment  and  0.15  second  per  iteration. 

IN-2  One  Story  Frame  with  Fixed  Base 

The  one  story-one  bay  frame  shown  in  Fig.  5-2  was  tested  by 
Arnold  et  at.  (1968).  The  mechanical  and  member  properties  are  given 
in  Tables  5-1  and  5-2  respectively.  Two  sets  of  load  were  applied  to 
the  frame.  The  vertical  set  of  load  was  kept  constant  while  increasing 
the  lateral  load,  H.  Figure  5-4  illustrates  the  results  obtained  using 
different  formulations.  The  load-deformation  curves  (1),  (2),  (3)  and 
(4)  represent  the  inelastic  analysis  presented  in  this  chapter,  the 
experimental  test  results,  the  second  order  elastic-plastic  analysis 
and  the  first  order  elastic-plastic  analysis,  respectively.  In  the 
inelastic  analysis  the  effects  of  residual  stresses  and  strain  hardening 
were  considered.  These  effects  are  not  included  in  the  elastic-plastic 
analyses  which  were  performed  using  the  iterative  method,  described  in 
Sect.  4.3.  The  results  show  good  agreement  between  the  behavior  of 
the  frame  tested  and  the  inelastic  analysis. 
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FIGURE  5-2  Example  Frame  IN-2 
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FIGURE  5-3  Example  Frame  IN-3 
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Section 
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TABLE  5-1  Mechanical  Properties 


Section 

b, 

in 

inches 

d, 

in 
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f, 

in 
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w, 
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b 

t 
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Z, 

in 
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38.57 

1070.6 

— 

28.0 

5WF18.5 

5.025 

5.12 

0.420 

0.265 

11.96 

19.3 

104.5 

2.16 

1.28 

56.17 

640.34 

558.4 

1 1.4 

TABLE  5-2  Member  Properties 
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FIGURE  5-4  Load-Deformation  Curves  of  Frame  IN-2 
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IN- 3  One  Story  Frame  With  Hinged  Base 

Schilling  et  al.  (1956)  tested  the  frame  shown  in  Fig.  5-3. 
During  the  test  the  applied  loads  on  the  frame  were  increased  propor¬ 
tionally.  Good  agreement  between  the  test  results  (curve  2)  and  the 
inelastic  analysis  (curve  1)  is  observed  from  Fig.  5-5.  The  figure 
also  shows  the  first  order  elastic-plastic  analysis  (curve  4)  and  the 
second  order  elastic-plastic  analysis  (curve  3) ,  obtained  using  the 
iterative  solution  presented  in  Sect.  4.3. 

In  examples  IN-2  and  IN-3  about  twelve  load  increments  were 
needed  in  the  inelastic  analysis  to  obtain  the  load-deformation  character 
istics  of  the  frames.  The  average  number  of  iterations  was  5  iterations 
for  each  load  increment  with  an  average  of  0.25  second  per  iteration. 

5.7.2  Comments  on  Results 

Results  obtained  from  the  inelastic  analysis  compare  favorably 
to  the  experimental  test  results  available.  The  results  show  the 
capability  of  the  computer  program  to  analyze  inelastic  frame  problems 
efficiently. 
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FIGURE  5-5  Load-Deformation  Curve  of  Frame  IN-3 


CHAPTER  VI 


BEHAVIORAL  STUDIES 

In  Chapter  5  an  inelastic  stability  analysis  based  on  a  stiffness 
formulation  was  presented.  The  validity  of  the  theory  was  then 
evaluated  by  comparing  the  results  with  some  sample  problems.  The 
primary  objective  of  this  chapter  is  to  compare  the  results  obtained 
using  the  formulation  given  in  Chapter  5  with  those  of  the  simple 
plastic  theory  given  in  Chapter  4  (i.e.  elastic-plastic  analysis). 

To  this  end,  a  cantilever  beam-column  and  single  story  frames  have 
been  analyzed. 

In  addition  to  the  basic  study  outlined  above,  the  effects 
of  residual  stresses  and  strain  hardening  on  the  behavior  of  single 
story  and  multistory  frames  are  examined.  Although  these  effects 
have  been  investigated  by  many  researchers,  their  influence  on 
frame  behavior  and  strength  has  not  been  considered  in  detail. 

Mild  steel  with  a  yield  stress  of  36  ksi  is  used  in  the 

following  analyses,  unless  otherwise  specified.  A  maximum  compressive 

residual  stress  of  an  =  0.3  a  is  considered,  and  the  residual  strain 

R  y 

distribution  shown  in  Fig.  D-lb  is  assumed.  A  strain  hardening  modulus 
of  900  ksi  is  used,  and  the  strain-hardening  strain  is  assumed  equal 
to  twelve  times  the  yield  strain.  Although  the  strain  hardening 
modulus  is  rather  high  it  is  thought  to  represent  commercially 
available  structural  steel  grades. 
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6 . 1  Behavioral  Study  of  a  Beam-Column 


The  cantilever  beam-column  shown  in  Fig.  6-1  is  used  in 
this  study.  The  axial  force,  P,  is  kept  constant,  while  the  lateral 
load,  Q,  is  increased  monotonically . 

6.1.1  Inelastic  Analysis  Vs.  Rajasekaran' s  Method 

This  example  problem  was  previously  analyzed  by  Rajasekaran 

(1971)  for  a  slenderness  ratio  (— )  equal  to  14.4  and  P/P  of  0.6. 

r  y 

Rajasekaran  used  an  iterative  incremental  technique  based  on  an 
equilibrium  balance  to  solve  for  inelastic  beam-column  problems. 
Figure  6-2  shows  the  load-deformation  curves  that  were  obtained 
using  different  solution  approaches.  It  should  be  noted  that  in 
the  inelastic  analysis,  the  residual  stresses  are  not  considered. 

The  results  obtained  by  Rajasekaran  are  identical  to  curves  3  and  4 
which  were  obtained  using  the  present  analysis. 

6.1.2  Inelastic  Vs.  Elastic-Plastic  Analysis 

To  compare  the  inelastic  method  to  the  elastic-plastic 
solution,  the  cantilever  shown  in  Fig.  6-1  was  analyzed  using  each 
method.  In  each  approach  the  axial  force  was  kept  constant  and  the 
corresponding  ultimate  value  of  Q  is  calculated  using  the  inelastic 
as  well  as  the  elastic-plastic  analysis.  In  the  inelastic  approach, 
however,  residual  stresses  and  strain  hardening  are  not  considered, 
as  they  are  not  accounted  for  in  the  elastic-plastic  formulation.  Th 
ratio  p  between  the  ultimate  values  of  Q  obtained  by  the  two  methods 
has  then  been  computed  as 
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FIGURE  6-1  Cantilever  Beam  -  Column 


L/r  =  14.4 
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FIGURE  6-2  Inelastic  Vs.  RajaseKaran  (1971) 
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P  = 


(Vin 

(Vel-pl 


x  100  (in  percent) 


(6.1.1) 


in  which  (Q  ) .  and  (Q  )  .  ,  are  the  ultimate  values  of  Q  obtained 

u  in  u  el-pl 

using  the  inelastic  and  the  elastic-plastic  analysis,  respectively. 
Figure  6-3  shows  the  relationship  between  p  and  P/P^  for  different 
values  of  L/r.  The  results  show  that  the  ultimate  lateral  load  as 
predicted  by  the  inelastic  analysis  is  less  than  that  obtained  from 
the  elastic-plastic  approach.  The  plot  also  shows  that  p  is  a 
function  of  L/r  and  P/P  . 

y 

The  conclusions  that  can  be  drawn  from  Fig.  6-3  may  be  mis¬ 
leading.  The  figure  actually  compares  the  lateral  strength  of  the 
cantilever,  as  predicted  by  the  inelastic  solution,  to  that  obtained 
from  the  elastic-plastic  approach.  The  overall  strength  of  the 
cantilever  as  predicted  by  both  methods  is  analyzed  more  appropriately 
through  interaction  curves  as  shown  in  Fig.  6-3a.  In  these  curves 

the  nondimensionalized  first  order  moment,  QL  /M  ,  at  the  maximum 

c  p 

strength  of  the  beam-column,  is  plotted  against  the  nondimensionalized 
axial  load,  P/P  .  Figure  6-3a  indicates  that  the  difference  between 
the  overall  strength  given  by  the  inelastic  method  and  that  obtained 
from  the  elastic-plastic  analysis  is  very  small.  The  reason  for  the 
difference  in  results  between  the  two  solutions  will  be  examined 


in  Chapter  7  (Sect.  7.2) 


6.1.3  Effects  of  Residual  Stresses  and  Strain  Hardening 

A  cantilever  beam-column  with  different  L/r-values  is  used 
to  study  the  effects  of  residual  stresses  and  strain  hardening. 

The  data  are  illustrated  in  Figs.  6-4  through  6-7.  These  figures 
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"IGURE  6-3 


Inelastic  Vs.  Elastic-Plastic  (Cantilever) 
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Figure  6-3a.  Column  Interaction  Curves 
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FIGURE  6-4  Load-Deformation  Characteristics;  Cantilever 
(L/r  =  14.4;  P/Py  =  0.15) 
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FIGURE  6-6  Load-Deformation  Characteristics;  Cantilever 
(L/r  =  80;  P/Pv=  0.15) 
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FIGURE  6-7  Load-Deformation  Characteristics;  Cantilever 
(L/r  =  80;  P/Pv  =  0.3) 
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show  that  when  only  the  lateral  strength  is  considered,  residual 
stresses  will  have  significant  effect  on  the  maximum  lateral  capacity, 
especially  for  members  with  large  axial  loads.  The  results  of  the 
ultimate  strength  analysis,  including  the  effect  of  residual  stresses 
are  plotted  as  dashed  lines  in  Fig.  6-3a.  The  results  show  that 
when  the  overall  strength  is  considered,  the  effect  of  residual 
stresses  is  relatively  insignificant. 

Except  for  the  case  of  low  values  of  L/r  and  P/P^  (Fig.  6-4) , 
the  data  demonstrate  that  strain  hardening  of  the  material  has  no 
effect  on  the  maximum  value  of  Q.  Actually,  it  was  found  that  the 
strain  hardening  tends  to  affect  only  the  slope  of  the  descending 
line  of  the  load-deformat ion  characteristics.  However,  for  the  case 
of  L/r  =  14.4  and  P/P^  =  0.15,  the  effect  is  pronounced.  This  can  be 
attributed  to  the  low  slenderness  ratio  of  the  member. 

6 . 2  Behavioral  Study  of  Single  Story  Frames 
6.2.1  Inelastic  Vs.  Elastic-Plastic  Analysis 

A  study  similar  to  that  of  Sect.  6.1.2  has  been  performed  for 
the  frame  shown  in  Fig.  6-8,  to  compare  the  inelastic  and  the  elastic- 
plastic  analyses  when  applied  to  frames.  The  results  are  then 
illustrated  on  plots  similar  to  Fig.  6-3. 

As  in  the  case  of  the  cantilever,  the  elastic-plastic  method 

tends  to  overestimate  the  lateral  strength  of  the  frame  (Fig.  6-9). 

The  results  of  the  ultimate  strength  analysis  are  given  in  the  form 

of  interaction  curves,  where  the  nondimensionalized  first  order 

moment,  QL  /2M  ,  at  the  maximum  strength,  is  plotted  again  the  non- 
c  p 

dimensionalized  axial  load,  P/Py.  Three  different  values  of  are 
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FIGURE  6-8  Single  Story  Frame 
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FIGURE  6-9  Inelastic  Vs  Elastic  -  Plastic 
(Single  Story  -  Hinged  Bases) 


considered  (0,  1,  and  3),  where  GT  is  evaluated  at  the  top  joint 
from  Eq.  3.3.1.  The  results  are  given  in  Fig.  6-10  (L/r  =  40)  and 
Fig.  6-11  (L/r  =  60).  The  effective  length  factors  corresponding  to 
G^  values  of  0,  1,  and  3,  with  the  bottom  end  of  the  column  pinned 
(G  =  °°) ,  are  2.0,  2.3,  and  2.92,  respectively.  This  gives  a  range 

D 

of  effective  slenderness  ratios  from  KL/r  =80  (K  =  2.0,  L/r  =  40) 
to  KL/r  =  175.2  (K  =  2.92,  L/r  =  60). 

Generally,  the  results  show  that  the  strength  decreases  as 
the  slenderness  ratio  increases.  Also,  the  results  of  Figs.  6-10 
and  6-11  indicate  that  the  difference  in  overall  strength,  as 
predicted  by  the  two  methods,  is  insignificant. 

6.2.2  Effects  of  Residual  Stresses  and  Strain  Hardening 

The  effects  of  residual  stresses  and  strain  hardening  are 
evaluated  in  a  study  of  the  frame  shown  in  Fig.  6-8.  The  gravity 
loads  are  kept  constant,  and  the  relationship  between  Q  and  A/L  is 
developed  using  different  values  of  L/r  and  P/P^.  The  results  are 
given  in  Figs.  6-12  through  6-15. 

It  is  clear  that  the  residual  stresses  have  a  significant 
effect  on  the  lateral  strength  of  the  frame,  particularly  for  frames 
with  high  gravity  loads  (high  P/P^  values).  In  Fig.  6-12,  where 
L/r  =  40  and  P/P  =  0.15,  the  reduction  in  ultimate  lateral  load  due 

y 

to  residual  stresses  is  about  2%.  When  the  value  of  P/P  is 

y 

increased  to  0.6  (L/r  =  40)  the  ultimate  lateral  load  of  the  frame 
(Q^)  is  reduced  by  about  30%  (Fig.  6-13).  The  reduction  in  ultimate 
lateral  load  is  about  4%  for  the  frame  with  L/r  =  60  and  P/P^  =  0.15 
(Fig.  6-14).  When  P/P  is  increased  to  0.4,  for  the  same  frame 
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Figure  6-10.  Maximum  Strength  of  Beam-Columns 
in  Portal  Frames  (L/r  =  40) 
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Figure  6-1 1.  Maximum  Strength  of  Beam-Columns  in  Portal  Frames  (L/r  =  60) 
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FIGURE  6-12  Load-Deformation  Characteristics;  Single  Story  Frame 

(L/r  =  40;  P/Py=  .15) 
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FIGURE  6-13  Load-Deformatin  Characteristics;  Single  Story  Frame 

(L/r  =  40;  P/Py  =  0.6) 
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FIGURE  6-14  Load-Deformation  Characteristics;  Single  Story  Frame 

(L/r=  60;  P/Pv=  0.15) 
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FIGURE  6-15  Load-Deformation  Characteristics;  Single  Story  Frame 

(L/r  =  60;  P/Pv  =  0.4) 
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(L/r  =  60) ,  the  lateral  load  capacity  is  reduced  by  about  35% 
(Fig.  6-15). 


The  interpretation  of  the  results  shown  in  Figs.  6-12  through 
6-15  may  be  misleading.  In  these  figures  only  the  lateral  load 
capacity  of  the  frame  is  considered.  The  results  of  the  ultimate 
strength  analysis,  including  the  effect  of  residual  stresses,  are 
plotted  as  dashed  lines  in  Figs.  6-10  and  6-11.  The  results  show 
that  the  effect  of  residual  stresses  is  actually  small  in  comparison 
with  the  overall  capacity.  The  figures  also  show  that  the  residual 
stress  effect  decreases  as  the  level  of  axial  load  decreases,  and 
as  the  slenderness  ratio  increases. 

Figures  6-12  through  6-15  show  that  strain  hardening  of  the 
material  have  no  effect  on  the  ultimate  lateral  load  of  the  frame. 

It  may  have  an  effect  only  on  the  slope  of  the  descending  portion 
of  the  load-deformation  characteristic.  The  change  in  slope  is 


particularly  noticeable  for  frames  with  low  values  of  P/Py. 


6 . 3  Behavioral  Study  of  Multistory  Frames 
6.3.1  Frames  and  Loading 

The  formulation  developed  in  Chapter  5  is  intended  primarily 
for  the  analysis  of  multistory  frames.  In  the  following,  the  three 
frames  shown  in  Figs.  6-16  through  6-18  are  used  for  this  purpose. 
Material  properties,  column  lengths,  and  girder  lengths  are  shown 
for  each  frame.  The  frames  are  designated  by  a  numeral  indicating 
the  number  of  stories,  followed  by  a  numeral  indicating  the  number 
of  bays.  The  loading  on  each  frame  is  the  service  load,  and  it  is 
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FIGURE  6-16  Configuration  and  Working  Loads;  Frame  8-1 
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FIGURE  6-17  Configuration  and  Working  Loads;  Frame  6-2 
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FIGURE  6-18  Geometry;  Frame  24-3. 


■ 


165 


assumed  to  be  increased  proportionally  until  collapse.  Frames  8-1 
and  6-2,  analyzed  by  Korn  and  Galambos  (1968),  are  subjected  to 
concentrated  loads,  as  shown  in  Figs.  6-16  and  6-17.  The  loading  for 
both  frames  is  given  in  ton  units  (1  ton  =  2.24  kips).  Frame  24-3, 
analyzed  by  Parikh  (1966),  is  subjected  to  distributed  loads  given 
in  Table  6-1.  Beam  and  column  sections  for  Frame  24-3  are  given  in 
Tables  6-2  and  6-3,  respectively. 

6.3.2  Presentation  of  Results 

The  load-deformation  characteristics  for  the  three  frames 
are  shown  in  Figs.  6-19,  6-20,  and  6-21.  Each  frame  was  analyzed 
using  different  methods  of  analysis,  as  indicated  on  the  figures. 

It  can  be  observed  that  residual  stresses  have  a  limited  effect  on 
the  ultimate  capacity  of  frames.  However,  they  have  a  slight  effect 
on  the  load-deformation  curve  during  the  transition  from  the  elastic 
into  the  fully  plastic  condition.  It  is  believed  that  the  effect 
of  residual  stresses  for  large  frames  will  become  more  significant 
with  increasing  gravity  loads. 

The  effect  of  strain  hardening  was  found  to  be  more  signi¬ 
ficant  for  large  frames  than  for  single  story  frames.  When  considered, 
it  caused  an  increase  in  the  ultimate  load  ranging  from  5  to  10%.  This 
is  due  to  the  large  number  of  sections  that  have  to  be  plastified 
in  order  for  the  frame  to  reach  its  collapse  load.  When  plasticity 
develops  in  the  last  few  critical  sections,  the  strain  at  a  sufficient 
number  of  the  sections  that  reached  the  full  plastic  moment  earlier 
may  have  attained  the  strain  hardening  value.  This  does  not  occur 
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Unit 

Working 

Load 

Roof 

k/f  t 

3.00 

Floor  AB 

k/ft 

4.36 

Floor  BC 

k/f  t 

4.73 

Floor  CD 

k/ft 

4.06 

Wind 

K 

5.76 

Wall 

K 

24.5 

Column 

K 

7.5 

TABLE  6-1  Loading;  Frame  24-3 
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[^^Beams 

Floor'''-- 

AB 

BC 

CD 

1 

W  14  x  26 

M  12  x  11.8 

W  16  x  45 

2 

W  16  x  36 

W  12  x  16.5 

W  8  x  55 

3 

do 

do 

do 

4 

do 

do 

do 

5 

do 

W  16  x  31 

do 

6 

W  16  x  45 

W  16  x  40 

do 

7 

do 

do 

do 

8 

W  18  x  50 

W  18  x  50 

do 

9 

do 

do 

do 

10 

W  21  x  55 

W  21  x  55 

W  21  x  55 

11 

do 

do 

do 

12 

W  21  x  62 

W  21  x  62 

W  21  x  62 

13 

do 

do 

do 

14 

W  21  x  68 

W  21  x  68 

W  21  x  68 

15 

do 

do 

do 

16 

W  24  x  68 

W  24  x  68 

W  24  x  68 

17 

do 

do 

do 

18 

W  24  x  76 

W  24  x  76 

W  24  x  76 

19 

do 

do 

do 

20 

do 

do 

do 

21 

W  24  x  84 

W  24  x  84 

W  24  x  84 

22 

do 

do 

do 

23 

do 

do 

do 

24 

W  27  x  84 

W  27  x  84 

W  27  x  84 

TABLE  6-2  Beam  Sections;  Frame  24-3 
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'^•^Column 

Story^-"-^ 

A 

B 

C 

D 

1 

- 

3 

W 

12 

X 

40 

W 

12 

X 

40 

W 

12 

X 

58 

W 

12 

X 

58 

3 

- 

5 

W 

12 

X 

58 

W 

12 

X 

58 

W 

12 

X 

79 

W 

12 

X 

79 

5 

- 

7 

W 

14 

X 

78 

w 

14 

X 

78 

w 

14 

X 

111 

W 

14 

X 

111 

7 

- 

9 

W 

14 

X 

111 

w 

14 

X 

111 

w 

14 

X 

136 

W 

14 

X 

136 

9 

- 

11 

W 

14 

X 

127 

w 

14 

X 

127 

w 

14 

X 

158 

W 

14 

X 

158 

11 

- 

13 

w 

14 

X 

136 

w 

14 

X 

142 

w 

14 

X 

193 

W 

14 

X 

158 

13 

- 

15 

w 

14 

X 

142 

w 

14 

X 

167 

w 

14 

X 

211 

w 

14 

X 

184 

15 

- 

17 

w 

14 

X 

167 

w 

14 

X 

193 

w 

14 

X 

246 

w 

14 

X 

202 

17 

- 

19 

w 

14 

X 

211 

w 

14 

X 

237 

w 

14 

X 

314 

w 

14 

X 

246 

19 

- 

21 

w 

14 

X 

246 

w 

14 

X 

264 

w 

14 

X 

342 

w 

14 

X 

287 

21 

- 

23 

w 

14 

X 

287 

w 

14 

X 

314 

w 

14 

X 

370 

w 

14 

X 

314 

23 

- 

25 

w 

14 

X 

314 

w 

14 

X 

342 

w 

14 

X 

398 

w 

14 

X 

320 

TABLE  6-3  Column  Sections;  Frame  24-3 
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FIGURE  6-19  Load-Deformation  Characteristics;  Frame  8-1. 
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FIGURE  6-20  Load-Deformation  Characteristics;  Frame  6-2. 
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FIGURE  6-21  Load-Deformation  Characteristics;  Frame  24-3 
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in  the  single  story — single  bay  frame,  as  the  critical  sections 
develop  plasticity  within  a  short  time  of  one  another. 

It  is  important  to  note  that  the  results  of  the  first  and 
second  order  elastic-plastic  analyses  for  Frames  8-1  and  6-2, 
obtained  using  the  iterative  procedure  described  in  Chapter  4,  are 
identical  to  those  obtained  by  Korn  and  Galambos  (1968) .  The 
results  obtained  by  Parikh  (1966)  for  Frame  24-3  are  very  similar 
to  curve  (2)  of  Fig.  6-21  (up  to  the  peak  load).  It  is  noted  that 
Parikh  took  into  consideration  residual  stresses  in  the  columns. 

In  Figs.  6-19  through  6-21,  a  comparison  between  the 
inelastic  response  (curve  1)  and  the  elastic-plastic  response 
(curve  5)  shows  a  slight  difference  only  near  the  ultimate  load. 
This  difference,  ranging  between  2  to  4%,  may  be  attributed  to  the 
gradual  penetration  of  yielding.  If  the  strain  hardening  effect 
is  included  in  the  inelastic  analysis,  there  is  no  apparent 
difference  between  ultimate  loads.  Thus,  when  the  overall  behavior 
of  the  frame  is  considered,  local  losses  in  frame  stiffness  due  to 
inelastic  action  are  often  compensated  for  by  the  behavior  of  other 
members  and  by  the  effects  of  strain  hardening. 


CHAPTER  VII 


SUMMARY  AND  CONCLUSIONS 

An  investigation  of  inelastic  behavior  of  multistory  steel 
frames  has  been  presented.  The  investigation  includes  a  review  of 
stability  theory,  with  a  discussion  of  the  different  approaches  to 
formulating  stability  and  strength  analyses.  An  elastic  buckling 
formulation  with  a  solution  technique  to  determine  the  elastic 
critical  load  of  frames  is  also  given. 

A  nonlinear  method  of  frame  analysis,  based  on  a  simplified 
large  deformation  theory  that  is  applicable  to  both  inelastic  and 
elastic  responses  has  been  developed.  A  finite  element  model  has 
been  used  for  this  purpose,  and  a  solution  method  to  solve  for  the 
response  of  frames  using  the  finite  element  equations  is  discussed. 

In  the  inelastic  analysis  the  effect  of  axial  loads  on  the  stiffness 
of  the  structure  is  considered.  Gradual  penetration  of  yielding 
through  the  cross  section,  the  spread  of  inelastic  zones  along  the 
member  length,  the  presence  of  residual  stresses,  and  strain 
hardening  of  the  material  are  also  accounted  for. 

The  investigation  also  includes  the  features  of  the  elastic- 
plastic  response  of  frames  for  first  and  second  order  analyses.  The 
incremental  and  the  iterative  solution  techniques  are  discussed  in 
this  context. 

During  the  course  of  the  investigation  computer  programs 
have  been  developed  for  all  types  of  analyses.  A  variety  of  problems 
are  solved  to  determine  the  accuracy  and  efficiency  of  the  proposed 
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techniques  and  to  demonstrate  their  performance  in  comparison  with 
other  methods. 

The  behavior  of  a  number  of  frames  under  the  action  of 
combined  vertical  and  lateral  loads  has  been  presented.  Comparisons 
are  made  between  the  inelastic  and  the  elastic-plastic  methods  of 
analysis.  The  effects  of  residual  stresses  and  strain  hardening  are 
also  studied. 

The  following  sections  outline  the  major  conclusions. 

7 . 1  Nonlinear  Strain-Displacement  Relationship 

In  an  energy  formulation,  deformations  affect  the  equilibrium 
equations  through  the  nonlinear  terms  in  the  strain-displacement 
equations  for  virtual  strains.  Additional  nonlinear  terms  must  also 
be  included  in  the  formulation  if  the  large  geometric  effects  on 
stresses  are  to  be  considered.  Thus,  in  the  formulation  presented 
in  Sect.  3.5.1,  the  effect  of  the  deformations  on  the  equilibrium 
equations  arises  from  the  nonlinear  terms  of  the  expression  for  6ez 
in  the  virtual  work  equation  (Eq.  3.5.7).  The  nonlinear  terms  as 
they  affect  stresses  arise  from  the  nonlinear  strain-displacement 
equations.  Thus,  a  stability  formulation  arises  from  the  first  effect, 
while  the  second  effect  simply  gives  a  more  accurate  estimate  of 
stresses  in  the  deformed  configuration  of  the  structure. 

7.2  Inelastic  Vs.  Elastic-Plastic  Analysis 

The  elastic-plastic  method  of  analysis  is  found  to  over¬ 
estimate  the  ultimate  capacity  of  frames  under  the  action  of  lateral 
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loads,  together  with  large  axial  loads  and  column  slenderness  ratios. 
When  the  overall  capacity  is  considered,  the  difference  between  the 
inelastic  method  and  the  elastic-plastic  solution  is  fairly  small. 

This  difference  arises  primarily  because  of  two  factors. 

Firstly,  the  assumed  interaction  curve  in  the  elastic-plastic 
analysis  (Fig.  4-1)  overestimates  the  capacity  of  the  section.  This 
curve  is  applicable  only  for  strength  analysis,  since  it  applies  to 
the  case  of  L/r  =  0.  However,  the  actual  ultimate  strength  inter¬ 
action  curves  are  functions  of  the  slenderness  ratio,  L/r  (Fig.  2-10). 
The  true  interaction  curve  values  are  the  ones  calculated  as  stress 
resultants  in  the  inelastic  method  of  analysis  as  presented  in  this 
investigation. 

Secondly,  the  P-A  moments  which  are  included  in  the  elastic- 
plastic  analysis,  through  the  geometric  stiffness  matrix,  are  under¬ 
estimated.  For  a  cantilever  member  with  an  axial  compressive  force, 

P,  and  a  sway  displacement  at  the  top.  A,  the  P-A  method  considers 
the  maximum  secondary  moment  at  the  bottom  of  the  cantilever.  This 
implies  a  triangular  distribution  of  secondary  moments  along  the 
member.  The  true  distribution,  however,  is  a  parabola  with  zero 
moment  at  the  top  and  PA  at  the  bottom.  This  difference  becomes 
more  significant  for  members  with  large  axial  loads  and  slenderness 
ratios . 

7 . 3  Gradual  Penetration  of  Yielding 

As  the  results  of  the  examples  presented  in  Chapter  5  show, 
the  load-deformation  curves  obtained  using  the  inelastic  and  the 


176 


elastic-plastic  methods  are  almost  identical  for  frames  with  low 
L/r  values  and  P/P^  ratios.  This  indicates  that  gradual  penetration 
of  yielding  has  a  negligible  effect  on  the  behavior  of  the  frames. 
However,  for  the  large  frames  presented  in  Sect.  6.3,  there  is  a 
difference  between  the  two  analyses  as  the  ultimate  load  is  approached 
This  may  be  attributed  to  the  gradual  penetration  of  yield.  The 
effect  is  apparent  in  large  frames  because  of  the  large  number  of 
sections  that  develop  fully  plastic  behavior  before  the  ultimate 
frame  load  is  reached.  In  single  story  frames,  however,  only  a 
few  sections  had  to  be  plastified  for  the  frames  to  form  mechanisms. 
Thus,  partial  plastif ication  of  sections  does  not  show  its  effect 


immediately  on  the  element  tangent  stiffness  matrix. 


7 . 4  Effects  of  Residual  Stresses  and  Strain  Hardening 

The  effect  of  residual  stresses  on  the  overall  strength  of  a 
frame  is  found  to  be  very  small.  It  decreases  as  the  axial  load 
decreases,  and  as  the  slenderness  ratio  increases. 

Strain  hardening  of  the  material  is  more  significant  for 
large  frames  than  for  single  story  frames.  Although  it  only  affects 
the  slope  of  the  descending  portion  of  the  load-deformation  curve 
for  a  single  story  frame,  it  causes  an  increase  in  the  ultimate 
strength  of  large  frames  by  up  to  10%. 


7 . 5  Incremental  Vs.  Iterative  Method  in  the  Elastic-Plastic  Analysis 


Although  more  expensive,  the  iterative  technique  is  found 
to  give  more  accurate  results  than  the  incremental  method.  However, 
good  results  can  be  obtained  if  an  iterative  procedure  is  used 
until  the  first  hinge  is  developed,  followed  by  an  incremental 
approach  until  the  collapse  mechanism  has  been  reached.  This 
procedure  has  been  called  the  iterative- incremental  technique,  and 
it  gives  approximately  the  same  results  as  the  iterative  method, 
but  the  cost  of  producing  the  solution  is  significantly  less  than 


that  of  the  latter. 
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APPENDIX  A 


ELEMENT  STIFFNESS  MATRIX  FOR 
BEAM-COLUMNS 


In  this  Appendix  the  "small-displacement"  element  beam-column 
stiffness  matrix  is  derived. 

A. 1  Virtual  Work  Equations  for  Linear  Stability  Problems 

To  formulate  the  virtual  work  equations  for  a  linear  stability- 
formulation  it  is  necessary  to  simplify  the  strain-displacement  equations 
presented  in  Sect .  3.2.  Assuming  that  9,  of  Fig.  3-1,  remains  "small" 
such  that 


1  -  cos0  i  0  (A. 1.1) 

Eqs.  3.2.3  become 

u  =  u  -  y  v'  (A. 1.2a) 

o  o 

v  =  v  ( A . 1 . 2b) 

o 

In  addition,  the  (u’)2  term  in  Eq.  3.2.4  is  assumed  "small"  with  respect 

to  the  other  two  terms  so  that  Eq.  3.2.4  becomes 

e  =  u'  +  i  (v')2  (A. 1.3) 

z  2 

Substituting  Eqs.  A. 1.2  into  Eq.  A. 1.3  yields 

e  =  u’  -  y  v"  +  \  (v’)!  (A. 1.4) 

z  o  J  O  2  o 
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Equation  A. 1.4  is  the  strain-displacement  equation  from  which  classical 
linear  stability  formulations  arise. 

The  principle  of  virtual  work  may  be  written  for  the  length 
of  beam  illustrated  in  Fig.  3-1  as 


6W 


* 

G  6e  dV  -  <Q>  { 6 q }  =  0 

z  z 


(A. 1.5) 


From  Eq.  A. 1.4 


6e  =  6u*  -  y  6vM  +  v’  6v'  (A. 1.6) 

Z  O  O  o  o 

Substituting  Eq .  A. 1.6  into  Eq .  A. 1.5;  integrating  over  the  area  of 
cross-section  A,  and  defining  the  stress  resultants  n  and  m  as 


n 


r 

G  dA 
z 

A  . 


and 


(A. 1.7a) 


m 


G 

z 


A 


ydA 


yields 


(A. 1.7b) 


6W 


(n  6u'  -m  5v"  +  n  vf  6vf)  dz  -  <Q>  { 6 q } 
o  o  o  o 


0  (A. 1.8) 


For  linear  elastic  response  it  is  assumed  that 

G  =  E  e  (A. 1.9) 

z  z 

Substituting  Eq.  A. 1.4  into  Eq.  A. 1.9,  and  the  result  into  Eqs.  A. 1.7, 
the  stress  resultants  may  be  expressed  in  terms  of  displacements  as 


. 
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n  =  EA  u'  (A. 1.10a) 

o 

and 

in  =  -  El  v"  (A.  1.10b) 

o 

in  which  I  is  the  moment  of  inertia  of  the  cross-section,  and  the  first 
moment  of  area  vanishes  because  the  axis  is  the  centroidal  axis.  A  term 
involving  (v^)7  has  been  discarded  from  Eq.  A. 1.10a  for  simplicity  of 
calculation. 

To  form  the  classical  equation  of  linear  stability,  Eqs.  A. 1.10 
are  substituted  into  the  first  two  terms  of  Eq.  A.1.8,  resulting  in  the 
equilibrium  requirement  that 


£ 


(EA  u’  6uT 
o  o 


+  El  vM  6vM)  dz 
o  o 


+ 


£ 


n  vf  6vf  dz  -  <Q>  { 6q }  =  0 

o  o 


(A. 1.11) 


Equation  A. 1.11  is  now  used  to  determine  finite  element  matrices. 


A.  2  The  Finite  Element  Model 


The  displacements  u^ 
displacement  coordinates  <q>. 


and  v^  may  be  approximated  in  terms  of  the 
of  Fig.  3-2a,  as 


u 

o 


(A. 2.1) 


and 


* 
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{ 

q. 


v  =  <d)> 

o 


(A. 2. 2) 


Standard  linear  shape  functions  (Fig.  A-l)  are  used  for  and 

standard  cubic  shape  functions  (Fig.  A-2)  are  used  for  <cj)>.  In 
terms  of  the  nondimens ional  coordinate  4,  shown  in  Figs.  A-l  and  A-2, 
the  shape  functions  are 


«f»  =  <  -  2  (C 


1)  ,-(£  +  !)> 


(A. 2. 3a) 


and 


«t»  =  <  \  (C  +  2)  (5  -  l)2  ,  |  (5  +  1)  (5  -  l)2, 

i  (2  -  ?)  (5  +  l)2  ,  |  (c  -  1)  (5  +  l)2  >  (A. 2. 3b) 


The  derivatives  required  for  Eq.  A. 1.11  are  now  evaluated  by 


differentiation  of  Eqs.  A. 2.1  and  A. 2. 2  to  yield 


u l  =  <  (<f>) 1  >  <  qx,  q4  > 


o 


(A. 2. 4a) 


v 


t  _ 


o 


<  4> '  >  <  q^y  q^,  q^ ,  q^> 


vM  =  <  tp"  >  <  q^y  q^ ,  q^,  q^  > 


(A. 2. 4b) 
(A. 2. 4c) 


in  which 


<  (4>)  ’>  = 


i  <  -1  ’  1  > 


(A. 2. 5a) 


<  4)'  >  = 


\<\  (35*  -  3)  ,  |  (3?2  -  2?  -  1)  , 


|  (3?2  -  3)  ,  |  (352  +25-1)  > 


(A. 2. 5b) 


<  >  =  j?  <  6?,  2(35  -  1)  ,  -  65,  2(35  +  1)  >  (A. 2. 5c) 


194 
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FIGURE  A-1  Basic  Functions  for  Linear  Displacements 
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A 

1 

'  f 


FIGURE  A-2  Basic  Functions  for  Cubic  Polynomial 
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Substituting  Eqs.  A. 2.4,  and  their  derivatives,  into  Eq.  A. 1.11; 
using  Eqs.  A. 2. 5;  carrying  out  the  integration;  and  recognizing  the 
arbitrary  nature  of  <6q>  ;  Eqs.  A. 1.11  may  be  written  as 


[k  ]  {q}  =  {Q} 

b 


(A. 2 . 6a) 


in  which  the  "secant  stiffness  matrix",  [k  ],  can  be  subdivided  into 

s 

the  "elastic  stiffness  matrix",  [k] ,  and  the  "geometric  stiffness  matrix", 
[k2] .  That  is 


‘V 


[k]  +  [k  ] 

o 


(A. 2. 6b) 


in  which 


[k]  = 


EA/£ 

0 

0 

-EA/£ 

0 

0 


12  El 
£3 

6  El 

£2 

0 

■12  El 

£3 

6  El 
£2 


4EI 

£ 

0 

■6  El 

2  El 

£ 


(SYM.) 


EA 

£ 

0 

0 


12  El 

£3 

-6  El 

£2 


4  El 

£ 


(A. 2. 6c) 


and 
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0 

0 


6 

5 

A 

10 

0 


(SYM.) 

2l2 

15 


0  0 


0 

0 


A 

5 

A 

10 


A 

10 

-A 

30 


0 

0 


6 

5 

A 

10 


2  £2 
15 


(A. 2 . 6d) 


The  matrix  [k]  is  the  standard  elastic  stiffness  matrix  commonly  used 

for  frame  analysis  (Ghali  and  Neville  1972) .  The  matrix  [k  ]  is  called 

§ 

the  geometric  stiffness  matrix  (Chajes  1974;  Porter  and  Powell  1971) 
and  is  a  function  of  the  element  axial  force. 
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APPENDIX  B 


DETAILS  OF  NONLINEAR  ELASTIC  FORMULATION 


This  appendix  presents  the  analytical  development  of  the  basic 
nonlinear  equations  used  in  this  study.  Also,  a  finite  element  model 
is  developed  and  the  finite  element  equilibrium  equations  are  formulated. 


B . 1  Basic  Equations 


In  Sect.  3.5.1  the  incremental  equilibrium  equation  for  an 
element  (Eq.  3.5.14)  was  written  symbolically  as 


[kT]  (Aq)  =  {AQ} 


(B.1.1) 


in  which  the  influence  coefficients  (kT)  _  of  the  tangent  stiffness 
matrix  were  evaluated  by  Eq.  3.5.16  which  may  be  rewritten  as 


(Vij 


3uf  3v’ 

(e.  .  -r-2-  +  e  -r-2-  +  e 

a  lj  8qi  2j  9qi 


9v" 

o 

3j  dq± 


)  •  dz 


(B.1.2) 


in  which  e. e0.,  and  e„ .  are  defined  as 
lj  2j  3j 


eu 


(B.1.3a) 


e2j 


^2 

dqi 


(B.  1 . 3b) 


and 


e„  . 

3j 


9a3 


(B . 1 . 3c) 
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The  quantities  a^,  a^ ,  and  a^  are  given  by  Eqs.  3.5.10.  To  evaluate 

the  e  terms  Eqs.  3.5.10  are  differentiated  with  respect  to  q.,  which 

J 

results  in 


(1  +  u')  |S- 
o  9q_ 


9uT 

+  n  -  v" 


9q . 


o  9q. 

J 


~  9v" 

9m  o 

”  m  r\ 

3qj 


(B . 1 . 4a) 


-2j 


rv  9V  ' 

=  V-  ^  + 

3qj 


v'  vM 
o  o 


o  9q . 
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2  + 


(V) 

o 


i  \  2 


/l-(v')2  ^  1  -  (vi)2J 

O 


9m 

3qj 


-  m 


v 


o 


2  + 


(v’) 

o 


1  -  (v') 

o 


9v" 

o 

9q, 


V 


If 


/l-(v’)1 2 

o 


»  \  4 


3(vf ) 
o 


o 


5(vo)2  1 

9v ' 
o 

v'  vM 
.  o  o 

(V-)2  1 

1  1  ° 

9m* 

1_  (vo)2J 

3qjJ 

i-(v')2 

o 

1-  (v’)2 
k  o  J 

3qj 

+  m* 


2  v'  v" 
o  o 


l-(v;)2 


1  + 


(v;) 


M2 


1 


i-(V)2 


3v" 


ms  2 


9q. 


(v") 
o  +  o 


1-  (v*)2 
o 


+ 


4 (v’ ) 4  5(v ’ ) 2  >  9v' 

_ 2 _  +  _ ° _  _ ° 

(1-  (v,)2)2  +  1-  (v*)2  9q  . 


(B . 1 . 4b) 


and 


•3j 


r  (v;)2 

°  /l-(v')2 


V" 

9q .  o 


M2 


1  + 


(v’) 
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l-(v’)^ 
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9m* 

3qj 


-  m 


9uf  v' 

3qj  /  1  -  (v^)  2 


(v')z 

2  +  1-  (v')z 
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9v' 

o 

9q 4 
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+  m* 


1  + 


{vo) 


9v"  2v'  v' 

o  .  o  o 


*  2 


l-(v^)2 


1  + 


(V) 

o 


9v' 


1  -  (vo>  2 


3qj 


(B . 1 . 4c) 


To  evaluate  the  incremental  force  AQ^,  of  Eq.  B.1.1,  Eq.  3.5.17  is 
used.  That  is 


AQ. 

l 


r  9u'  9v'  9v"'> 

o  ,  o  ,  o 

3.,  — —  +  a„  — —  +  a 


% 


1  9q 


2  9q.  3  9q. 

l  ni' 


dz 


(B.1.5) 


All  equations  derived  above  apply  to  elastic  and  inelastic 


response . 


B . 2  Nonlinear  Elastic  Response 


To  specialize  the  problem  to  elastic  response  the  material 

is  assumed  to  be  linearly  elastic  and  the  stress-strain  relationship 

obeys  Hooke's  law  (Eq.  A. 1.9). 

To  evaluate  the  stiffness  coefficients  (k  ) . .  and  the 

T  13 

incremental  forces  AQ_^,  it  is  necessary  to  evaluate  first  the  stress 

„  ,  ,  .....  9n  9m  ,  9m* 

resultants  n,  m  and  m*  and  also  their  derivatives  -z —  ,  -z —  ,  and  -r 

’  oq.  9q.  9q . 

3  3  3 

that  appear  in  Eqs .  3.5.10  and  B.1.4. 

Substituting  Eq.  A. 1.9  into  the  stress  resultants  defined  in 
Eqs.  3.5.9  yields 


n  = 


E  •  £  •  dA 

z 


(B . 2 . la) 
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m  = 


E  •  e  •  y  •  dA 
z 


(B .  2 .  lb) 


and 


m 


* 


E  •  e  •  y  dA 
z 


(B . 2 . lc) 


Let  the  cross  section  properties  for  a  symmetric  section  be 
defined  as 


dA  =  A 


(B . 2 . 2a) 


y  •  dA  =  0 


(B  .  2 . 2b) 


y2  dA  =  I 


(B . 2 . 2c) 


y  dA  =  0 


(B .  2 . 2d) 


y 4  dA  =  I 


(B.2.2e) 


Substituting  Eq.  3.5.5  into  Eqs .  B.2.1  and  integrating,  using 
the  definitions  given  by  Eqs.  B.2.2  leads  to 


n 


=  EA 


(u')2  (v')2- 

.,o,o 

uo  +  — 2 —  — 2 — 


+  j  El 


(v’)2  (v")2> 

(v")2  +  —2 _ 2_ 

'Vo;  l-Cv')2 


(B  .  2 . 3a) 
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m  =  -  El 


(v')2  V"  >! 

u'  V M  +  V"  +  - 2 - — 

°  0  °  /T  -  (v’)2 


(B  . 2 . 3b) 


and 


m*  =  El 


M2 


(u’) 

u'  +  — 2 
o  2 


+ 


(v')2^ 
o 


+  2  EI4 


(v')2  (v")2-. 


o 


(v7+tt-(v.)2 


(B.2.3c) 

Differentiating  Eqs.  B.2.3  with  respect  to  the  generalized 


displacement  yields 


and 


9n 

3qj 


=  EA 


9u' 


9v'  -i 


(1  +  u’) 


+  v' 


o  9q.  o  9q. 
1  J 


+  El 


r  K  (v">2 

o  o 

(v')2  1 
1  1  ° 

9v f 

°  +  v" 

r  (vM2 

i  ^  ° 

9v" 

o 

1-  (v')2 
o 

i-<v;)2J 

9q .  o 

J 

CM 

V 

-  o 

> 

1 

I— 1 

H 

QJ 

►a 

t_i. 

1 _ 

(B . 2 . 4a) 


9m 

9q  . 
J 


=  -El 


9u '  9v" 

^  +  (1  +  u’)  y-2-  + 
o  9q.  o  9q. 


.it 


v '  vM 
o  o 


/l-  (v')2 
o 


2  + 


(vp 


1-  (v')2 
o 


9v ' 
o 

8qj 


+ 


(Vq) 


/  1-  (v 1  )  2 
o 


9v" 

o 

9q . 
J 


(B . 2 . 4b) 


9m* 

3qj 


n 


EI  L(1  +  3q 


3u'  3v'-, 

O  ,  ,  o 

-  +  v'  -r - 

o  9q . 

J  -1 
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+  El, 


v' 


1  + 


(V) 

o 


1  *  (vo>2 


3v' 

^ —  +  v" 
9q .  o 


1  + 


(V;) 


1-  (v’)2 
O 


9v" 

o 

9q . 

J 


(B.2.4c) 


Equations  B.2.3  and  B.2.4  are  expressions  for  the  stress 
resultants  and  their  derivatives,  respectively,  derived  for  an  elastic 
element  subjected  to  axial  force  and  bending  moments. 

For  an  elastic  element  subjected  only  to  axial  force  (truss 
element)  the  bending  effect  is  eliminated  and  the  displacements  at 
any  point  along  the  element,  Eqs .  3.2.3,  become 


u  =  u 

o 


and 


v  =  v 

o 


Consequently 


(B . 2 . 5a) 


(B  .  2 . 5b) 


and 


u 


i  _ 


=  u 


o 


T  _ 


=  V 


o 


(B  .  2 . 6a) 


(B . 2 . 6b) 


Substituting  Eqs.  B.2.6  into  the  strain  expression  given  by  Eq.  3.2.4 
yields 


£ 

z 


K>2  + 


(B.2.7) 


Substituting  Eq.  B.2.7  into  Eqs.  B.2.1  and  integrating  using  Eqs.  B.2.2 
yields  the  stress  resultants 


■’  ;  '  !.  «*|  |  v.:J  :vTnxt*<jU  J 
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n  =  EA 


u 


(u’) 

o 


o 


+ 


(v ' ) 2— r 

o 


m  =  0 


and 


m*  =  El 


O')2  O')2- 

1,0,0 
11  ’  +  — - —  +  — - — 
o  2  2 


Differentiating  Eqs.  B.2.8  with  respect  to  q  yields 


3n 

3qj 


=  EA 


3u '  3v'-r 

(1  +  u’)  -T-2-  +  v’ 

o  dq  .  o  dq. 

J  1 


3m 

3qj 


=  0 


and 


3m* 

3q . 


=  El 


(1  +  up 


3u'  3v'- 

_ o  ,  o 

3q.  Vo  3q. 


J”1 


(B. 2.8a) 


(B .  2 . 8b) 


(B . 2 . 8c) 


(B  .  2 . 9a) 


(B  .  2 . 9b) 


(B. 2 . 9c) 


B . 3  Evaluation  of  Unbalanced  Forces  Vector  {Ar}  For  Finite  Element  Model 

The  incremental  force  AQ^  as  given  by  Eq.  B.1.5  can  be  evalu¬ 
ated  knowing  the  quantities  u1 ,  vf  and  v",  their  derivatives  with 

o  o  o 

respect  to  q^,  and  the  stress  resultants  n,  m  and  m*. 

The  quantities  u^ ,  v^  and  v”  can  be  evaluated,  from  Eqs. 

3.5.18  as 


u 


t  _ 


o 


<  <f>'  >  {u} 


(B . 3 . la) 
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v»  =  <  4> *  >  {v} 

o 


(B . 3 . lb) 


and 


v"  =  <  cj)n  >  {V} 

o 


(B  .  3  .  lc) 


Differentiating  Eqs.  B.3.1  with  respect  to  yields 


9u' 

o 

3qi 


4>T 

qi 


(B. 3 . 2a) 


9vf 

o 

9q  . 

l 


(B.3.2b) 


and 


9v" 

o 

3qi 


(B . 3 . 2c) 


in  which  (b  or  cb"  is  the  term  of  the  differentiated  shape  functions 
a . 

l  x 

associated  with  the  generalized  displacement  q^  in  { U }  or  {v}  as 
defined  in  Fig.  3-7. 


B.3.1  Vector  of  Incremental  Forces  {AQ}  of  Flexural  Element 
Substituting  Eqs.  B.3.2  into  Eq.  B.1.5  yields 


AQi 


=  Q,  - 


<J>* 

qi 


dz 


(B . 3 . 3a) 


for  i  =  1  to  4 ;  and 


AQ. 

l 


Qi 


i 


U2  ♦q. 

1 


+  a3  ^  ] 


dz 


(B. 3 . 3b) 


for  i  =  5  to  8;  and  the  quantities  a^,  a^  and  a^  are  defined  in 
Eqs.  3.5.10b,  3.5.10c,  and  3.5.10d,  respectively. 

Arranging  Eqs.  B.3.3  in  matrix  form  leads  to 


{AQ} 


{Q}  - 


& 


i-L 


a2  (cf)’ }  +  a^  {$"} 


dz 


(B.3.4) 


which  can  be  written,  symbolically,  as 
{AQ}  =  {Q}  -  {QR} 


(B.3.5) 


in  which  { AQ }  is  the  vector  of  incremental  forces  (8  x  1)  associated 
with  the  vector  of  nodal  displacements  defined  by  Eq.  3.5.19b,  {Q} 
is  the  vector  of  applied  loads  and  {Q  }  is  called  the  vector  of 
resisting  forces  and  is  defined  as 


(qr) 


n 


a^  {(J)’} 
a2  {(I)' }  +  a^ 


dz 


(B.3.6) 


To  evaluate  the  vector  {Q  }  Gaussian  integration  method  is 
used.  The  basic  assumption  in  Gaussian  numerical  integration  is  that 
(Bathe  and  Wilson  1976) 


b 

» 

{F(r ) }  dr  =  a.  (F(r.)}  +  a0  (F(r0)}+  ...+  a  {F(r  )} 

11  z  z  n  n 

1  (B.3.7) 


in  which  (F(r)}  is  the  matrix  to  be  integrated,  a,,  an ,  ...  and  a  are 

l  Z  n 

weighting  factors,  r^,  r2,  ...  and  r^  are  sampling  points  and  {F(ri>} 
is  the  matrix  (F(r)}  evaluated  at  the  ith  sampling  point. 


■ 
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According  to  Eq.  B.3.7, 


n 


I 

i=l 


a 


Eq.  B.3.6 

+  {(J)"} 


can  be  written  as 


i 


(B.3.8) 


in  which  n  is  the  number  of  sampling  points  (Gauss  points)  assumed 
along  the  element  and  a_^  is  the  weight  associated  with  the  iUi  Gauss 
point.  Many  references  have  tabulated  a  variety  of  sampling  points 
and  the  corresponding  weights  (Loxan  et  al.  1942;  Bathe  and  Wilson 
1976) . 


In  this  study  four  sampling  points  are  assumed  along  the 

element.  For  each  Gauss  point  assumed  the  matrices  { cf) T }  and  {({)"} 

are  evaluated.  Also,  the  quantities  u'  v’ ,  and  vM  are  calculated 

o  o  o 

using  Eqs.  B.3.1.  At  this  point  the  stress  resultants  can  be  evalu¬ 
ated  from  Eqs.  B.2.3.  Consequently  the  quantities  a^,  a2»  and  a^  can 
be  calculated  from  Eqs.  3.5.10b,  3.5.10c,  and  3.5.10d.  Substituting 
the  results  into  Eq.  B.3.8  and  using  the  appropriate  weight  yields 
the  contribution  of  one  Gauss  point  to  the  vector  {Q  }.  Summing  over 
all  the  Gauss  points  gives  the  vector  of  resisting  force  (Q  }  which 
can  be  subtracted  from  the  vector  of  applied  loads  { Q }  to  obtain  the 
vector  of  incremental  forces  {AQ}  as  defined  in  Eq.  B.3.5. 


B.3.2  Vector  of  Incremental  Forces  {AQ}  of  Truss  Element 

For  an  element  subjected  only  to  axial  force  the  linear  inter¬ 
polation  functions  of  Eqs.  A. 2. 3a  are  used.  Thus  Eqs.  B.3.2  become 


3u* 

o 


(0 


(B . 3 . 9a) 
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3v' 

o 

9q, 


=  <5„  )f 


qi 


(B  .  3 . 9b) 


and 


3v" 

o 

8qi 


=  0 


(B . 3 . 9c) 


Substituting  Eqs.  B.3.9  into  Eq.  B.1.5,  and  using  Eqs.  3.5.10b, 
3.5.10c,  and  3.5.10d,  yields 


AQi 


=  Q,  - 


£ 


n(l  +  u ' )  (d)  )  '  dz 

o'  vyq/ 


(B . 3 . 10a) 


for  i  =  1,2;  and 


AQi 


Qi  - 


n  v’  (cf)  )  ’  dz 
£  °  qi 


(B. 3.10b) 


for  i  =  3,4. 


Writing  Eqs.  B.3.10  in  matrix  form  leads  to 


{AQ}  =  {Q}  - 


£ 


n(l  +  u^)  {  ((}))'  } 
n  { ($) ' } 


dz 


(B.3.11) 


which  can  be  written  symbolically  as 


Ri 


{AQ}  =  {Q}  -  {Q  } 


(B.3.12) 


in  which  {AQ}  is  the  vector  of  incremental  forces  associated  with  the 
vector  of  nodal  displacements  defined  by  Eq.  3.5.19c  and  it  is  a  4  x  1 


vector . 
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Again  the  vector  {Q  }  can  be  evaluated  by  numerical  integration 
using  Gauss  method  explained  above. 

B.3.3  Assembly  of  Unbalanced  Forces 

To  solve  for  the  incremental  displacements  {Ar}  in  Eq.  3.5.24 
the  vector  of  unbalanced  forces  {Ar}  must  be  evaluated.  The  incremental 
forces  {AQ}  are  first  obtained  in  local  coordinate  system  by  Eq.  B.3.5. 
These  forces  can  be  assembled  into  the  global  coordinate  system  using 
Eq.  3.3.8  and  the  appropriate  transformation.  Assembling  Eq.  B.3.5 
or  Eq.  B.3.12  into  the  global  coordinate  system  yields 

{AR}  =  {R}  -  {RR}  (B.3.13) 

in  which  {Ar}  is  the  vector  of  unbalanced  forces,  {r}  is  the  assembled 
vector  of  applied  loads  and  {R  }  is  the  assembled  vector  of  resisting 
forces . 


B . 4  Evaluation  of  Stiffness  Matrix  [Kt] 


To  solve  for  the  incremental  displacements  {Ar},  in  Eq.  3.5.24, 
the  structure  tangent  stiffness  matrix  [K^]  must  be  evaluated.  In 
this  section  two  different  element  stiffness  matrices  are  developed, 
one  for  a  flexural  element  and  the  other  for  a  truss  element.  The 
element  stiffness  matrices  are  then  assembled  to  form  the  structure 
stiffness  matrix  [K^] . 


B.4.1  Element  Stiffness  Matrix  for  a  Flexural  Element 


The  influence  coefficient  (k^)^  must  be  evaluated  to  form 
the  element  stiffness  matrix  [kT] .  Substituting  Eqs.  B.3.2  into 


Eq.  B.1.2  yields: 


. 
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(a)  for  1  <  i  <  4 


(V±j 


u 


4>  ’  •  dz 


Z 


qi 


(B. 4. la) 


and 


(b)  for  5  <  i  <  8 


(Vij 


z 


[e?  •  <|>f  +  e  •  4>"  ]  dz 

2J  <l±  3J  q± 


(B .  4  .  lb) 


in  which  e1 . ,  e~  .  ,  and  e„  .  are  defined  in  Eqs.  B.1.4.  Thus,  the  ith 
lj  2j’  3j  H  J  — 

column  in  the  element  stiffness  matrix  [k^]  can  be  written  as 


Vi 


(4)'} 


li 


£ 


}  e2 j  +  ^"^3j 


dz 


(B.4.2) 


The  full  matrix  [k^]  can  be  written  as 


[kTl 


£L 


{(J)1}  <e1> 


(4)' }  <e2>d"  (4)"}  <e3> 


dz 


(B.4.3) 


or  it  can  be  written,  symbolically,  as 


[KjJ  = 


[\] 


[  k^  ]  +  [^3  ] 


(B.4.4) 


in  which  [1^]  is  an  8  x  8  matrix,  while  [k^]  ,  [k^]  and  [k^]  are  4x8 


matrices  defined  as 


' 
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[\] 


i 


{*'}  F(1  +  u')  <  I1  >  -  v"  <  I51  > 
o  dq  o  dq 


+  n  <  (j)|  >  -  in  <  $2  >j  •  dz 


(B . 4 . 5a) 


[k2] 


,  .  dn  . 

v'  <  -=—  > 
o  dq 


v ?  vM 
o  o 


/l-  (v’)2 


2  + 


o 


'l 


1”  (Vq)2, 


<  ^  > 
9q 


Vn  (Vn)2  r 

O  O 

1-  (v’)2 

o 


I  \  2 


1  + 


(V) 

o 


i-(v)2 


<  fr  >  -  n  "  *2  =• 


-  m 


I  \  z 


O 


l  /T-  (v;>2 


2  + 


(v') 

o 


i-  (V’)2 
O 


<  ^  > 


V 


It 

o 


(v;) 


i  2 


/I  -  (v')2  l1+  1  -  (v’)2 


2  + 


3^)2  ^ 


1  -  (v;)2j 


<  V2  > 


+  m* 


r 

2 

v’  v" 

r 

o  o 

1  _L_ 

1 

-  (v*)2 

l  1 

o 

i\2 


(v;) 


'')2, 
O  ' 


<  *”2  > 


+ 


(<) 

o 


1  \  2 


1  -  K  > 2 

O 


1  + 


(Vq) 


i  2 


1  -  <v;> 2 


i+ 


4(vo) 


- 

\ 

j  <  *2  > 

dz 


and 


(B . 4 . 5b) 


[k3] 


r 

(v')2  •, 

— 

1  +  u'  + - 2 - 

n 

1  0  y  i-  (v’)2J 

o 

<  I51  > 

3q 


t  a  2 


+  V 


IV 


1  + 


(v«) 

o 


1  -  (v;)  2 


<  > 
d  q 


-  m 


<  4>{  > 
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/  1-  (v')* 
o 


2  + 


(v') 

o 


1x2 


i-  (V»)2 
O 


<  > 


+  m* 


2  v'  v" 
o  o 


M2 


1-  (V)2 

o 


1  + 


(V) 

o 


i-(v;)2 


<  r2  > 


t  ^  2 


1  + 


1-  (V’)2 
o 


<  *”2  > 

1 — 

J 

dz 


(B. A . 5c) 


In  Eqs.  B.4.5,  <(j)|>,  and  are  1x8  row  vectors 


defined  as 


<  cj)|  >  =  <^<  4> f  >  •  <  0  >^> 


(B . 4 . 6a) 


<  t'2  > 


=  <(<  0  >  •  <  (f>* 


(B  .  4 . 6b) 


and 


<  ^2  > 


=  <^:  0  >  •  <  4>"  >^> 


(B .  4 . 6c) 


3  n.  3m  3m^ 

;  and  the  vectors  <-7T-> ,  >  and  <~r — >  are  also  1x8  row  vectors  which 

3q  dq  dq 

are  called  the  incremental  stress  resultant  vectors. 

In  order  to  be  able  to  use  the  Gaussian  integration  method 
to  evaluate  [k-^]  ,  [ ]  >  and  [k^]  the  incremental  stress  resultant  vectors 
must  be  evaluated  at  each  Gauss  point.  This  can  be  done  as  follows: 
Substituting  Eqs.  B.3.2  into  Eqs.  B.2.4  yields 

(a)  for  1  <  j  <  4 


3n 

3q, 


EA  (1  +  uf)  (j)’ 

°  qj 


(B.4.7a) 


■ 
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3m 

3qj 


-  EI  v"  <p' 

0  q. 

J 


(B.4. 7b) 


3m* 

3q . 

J 


EI  (1  +  u')  <j>' 

°  q. 

3 


(B.4.7c) 


and , 


(b)  for  5  <  j  <  8 


3n 

3qj 


EA  v'  (j)'  +  EI 

0  q3 


V* 

(v")2 

o 

o 

1  - 

(  V  *  )  2 

o 

1  + 


(v') 

o 


i  2 


(vo)2 


<K 


+  v" 
o 


1  + 


(v’) 

o 


2 


1-(V’)2 


(j)" 

q. 


(B . 4 . 7d) 


3m 

sir 


=  -  EI 


ll\  2 


V ?  (v") 
o  o 


/r  -jv'j? 

o 


2  + 


(vf) 

o 


t  \  2 


1_  (vo)2 


r 


qi 


2  i 

+ 

1  +  u*  +  - - 

( p " 

L  °  /l-(v')2J 

o 

q-iJ 

(B.4.7e) 


3m* 

3qj 


EI  v’  4>'  +  EI/ 

o  q .  4 

J  ' 


v1 

(v")2 

o 

o 

1  - 

(V’)^ 

o 

1  + 


(v;> 


M2 


i-(v;)2 


4>: 


?  \  2 


+  v" 
o 


1  + 


(v') 

o 


1  ■  (vo)2 


♦'q. 


(B . 4 . 7f ) 


Equations  B.4.7  can  be  written  in  a  column  vector  form  as 


listed  in  Table  B-l. 


. 


215 


/ 


3q 


■< 


EA  (1  +  u’)  {<J)'} 
o 


(  V*  (v") 2 

ea  v; +  ei  i  °  (vV 

o 


1  + 


(v')2 

o 


1  -  (v')2 
o 


+  EI  v" 
o 


1  + 


(vf ) 
o 


'l 


1  -  K)2 

o 


{<)>'} 


3q 


=  < 


-  EI  v"  1<J>'} 
o 


-  EI 


v?  v" 
o  o 


/  1-  (v')2 
o 


2  + 


(vo)2 


'I 


1  -  (v;)2 


> 


-  EI 


1  +  u'  + 

o 


(v’) 

o 


I  \  2 


/  1-  (v')2 
o 


TABLE  B-l  Incremental  Stress  Resultant  Vectors 


For  Flexural  Element 
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r9m*i 

1  9q 


El  (1  +  u^)  {(f)’} 


vf  (v")* 

El  v'  +  El  -°  (  °  2 
o  4  1  -  (v  ) 1 
o 


1  + 


(v') 

o 


TV2" 


1  “  (vo)2 


M 


+  El.  v 
4  o 


1  + 


(v’) 

o 


1  -  (v;)2 


{f } 


Note : 


is  8  x  1  column  matrix 
9q 


{(J)'}  and  {(J)"}  are  4x1  column  matrices 


TABLE  B-l  (Cont.)  Incremental  Stress  Resultant 

Vectors  for  Flexural  Element 
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Now  to  evaluate  [k^] ,  the  quantities  u^,  and  v”  are  calcu¬ 

lated  from  Eqs.  B.3.1  at  each  Gauss  point.  Thus  the  stress  resultants 
can  be  evaluated  from  Eqs.  B.2.3  and  also  the  incremental  stress 
resultant  vectors  can  be  evaluated  as  listed  in  Table  B-l.  At  this 
stage  the  Gaussian  integration  method  explained  above  is  used  to  eva¬ 
luate  Eqs.  B.4.5  leading  to  the  element  stiffness  matrix  given  by 
Eq .  B  A .  4 . 


B.4.2  Element  Stiffness  Matrix  for  a  Truss  Element 

For  a  truss  element  the  influence  coefficients  (km) . .  are 

T  ij 

obtained  by  substituting  Eqs.  B.3.9  into  Eqs.  B.1.2  to  yield 


(a)  for  1  <  i  <  2 


(Vu 


% 


(?  ) 


dz 


and 

(b)  for  3  <  i  <  4 


(B  .  4 . 8a) 


e  ((f)  )  '  dz 

£  2j  qi 


(B .4 . 8b) 


The  jth  column  in  the  element  stiffness  matrix  can  then  be  written 
as 


({($)’}  e  ) 

]{(?)')  e2.[  dZ 


(B.4.9) 


and  the  full  matrix  [kT]  can  be  written  as 


. 
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[kT]  = 


{(?)’}  <e1> 

{  ($)  f  }  <e2>_i 


dz 


(B.4.10) 


or  it  can  be  written,  symbolically,  as 


[kT]  = 


L[k2] 


(B.4.11) 


in  which  [k^,]  is  a  4  x  4  matrix  and  [k^]  and  [ ]  are  2x4  matrices 
defined  by  Eqs.  B.4.5a  and  B.4.5b.  In  such  a  case  {($)'}  replaces  {cf)'} 
and  it  is  a  2  x  1  vector.  Also  <(^)T>,  < (^2)  ? > »  anc^  <(<J)2)">  rePlace 
<(f)|>,  <4>2>»  and  > »  respectively.  In  this  case  <($^)’>,  <(^)'>,  and 
<(^)">  are  1x4  matrices  given  as 


and 


<  (4>1)  ’  >  =  (4>)’  >  •  <  0  > 

<  (?2)'  >  =  <^<  0  >  •  <  (5)’  > 

<  ($9)"  >  =  S<  0  >]<  0  > 


(B .4 . 12a) 
(B . 4 . 12b) 

(B .4 . 12c) 


in  which  <(j)>  is  the  shape  functions  for  a  truss  element  given  by 
Eq.  A. 2. 3a;  and  the  incremental  stress  resultant  vectors  are  1x4 


row  vectors. 


To  evaluate  the  incremental  stress  resultant  vectors 
Eqs.  B.3.9  are  substituted  into  Eqs.  B.2.9  to  obtain 


(a)  for  1  <  j  <  2 


3n 

3qj 


EA  (1  +  u')  (<f>  )  ' 

o  q  j 


(B.4.13a) 


. 
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(B.4.13b) 


and 


9ni*  — 

fr-  =  ei  (l  +  u»)  («(>-  )' 

dq  o  q 

•J  J 


(b)  for  3  <  j  <  4 


dn 

dq. 


=  EA  v'  ($  ) f 

°  q . 

J 


dm 

dq  . 
J 


=  0 


dm* 

3qj 


EI  vf  ((J)  )  ' 

°  qj 


(B . 4 . 13c) 


(B. 4 . 13d) 


(B .  4 . 13e) 


(B .4 . 13f ) 


Equations  B.4.13  can  be  written  in  a  column  vector  as  given 
in  Table  B-2. 

Substituting  the  incremental  stress  resultant  vectors  as 
given  in  Table  B-2  and  the  stress  resultants  as  given  in  Eqs.  B.2.8 
into  Eqs.  B.4.5a  and  B.4.5b,  the  matrices  [k^]  and  [k^] ,  in  Eq .  B.4.11, 
can  be  redefined  as 


[k1] 


{($)’} 


% 


(1  +  u')  <  >  +  n  <  ((j)_)f  > 

o  dq  1 


dz  (B.4.14a) 


and 


[k2] 


{($)’}  [v-  <  >  +  n  <  ($2)  ’  >  j  dz 


£ 


(B  .  4 . 14b) 
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{-r} 

8q 


I  EA  (1  +  u^)  {(<(.)'} 

EA  v'  {($)'} 
o 


{f1}  =  {  o  } 

dq 


r8m* 

l3q 


} 


EI  (1  +  u')  { (4>)T} 
_ o 

EI  {($)'} 


Note 


{-r— }  is  4  x  1  column  matrix 
3q 


{(<(’)'}  is  2  x  1  column  matrix 


Incremental  Stress  Resultant 
Vectors  for  Truss  Element 


TABLE  B-2 
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Integration  leads  to  the  element  stiffness  matrix  [k^]  in 
Eq.  B.4.11. 

B.4.3  Assembly  of  Structural  Stiffness  [K-p] 

The  assembly  of  element  stiffness  matrix  into  the  structure 
stiffness  matrix  is  explained  in  detail  in  Sect.  3.3.2. 
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APPENDIX  C 


LARGE  DISPLACEMENT  TRANSFORMATION 


Consider  the  element  shown  in  Fig.  C-l.  The  relation  between 
the  reference  axis  in  local  and  global  coordinate  systems  can  be 
written  as 


Z  =  z  cos  a  -  y  sin  a 


(C. la) 


and 


Y  =  z  sin  a  -  y  cos  a 


(C . lb) 


Also  the  relation  between  the  local  displacement  u  and  v  and 

•V  "V 

the  global  displacements  U  and  V,  at  any  point,  can  be  expressed  as 


u  =  U  cos  a  +  V  sin  a 


(C. 2a) 


and 


v  =  U  sin  a  -  V  cos  a  (C.2b) 

in  which  u  and  v  are  functions  of  z  and  y.  Differentiating  Eqs.  C.2 
with  respect  to  z  yields 


3u 

3z 


3U  3Z  ,  3V  3Z  .  ,  9U  3Y  ^  3V  3Y  . 

3Z  37  cos  01  +  3Z  37  Sln  “  +  W  37  cos  “  +  W  li  Sln  01 


(C. 3a) 


and 
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FIGURE  C-1  Element  Transformation 
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3v 

3z 


3U  3Z 
3Z  dz 


sin  a  - 


3V  3Z  .  3U  3Y 

3Z  3z  cosa  +  3Y  3z 


sin  a 


3V  _3Y 
3Y  3z 


cos  a 


(C.3b) 


Differentiating  Eqs.  C.l  results  in 


3Z 

-tt-  =  cosa 
3z 


(C.  4) 


3Y 

-7T—  =  sma 

3z 


(C.  5) 


Also  there  is  no  change  in  angle  between  the  axes  after  deformation. 
This  constraint  has  been  approximated  by  imposing  the  condition  that 


3U  =  _  3V 

3Y  3Z 


(C .  6) 


Substituting  Eqs.  C.4,  C.5,  and  C.6  into  Eqs.  C.3  results  in 


and 


3u 

3z 


3U  2 
3Z  cos  a  + 


3  V  .  2 

sm  a 


(C.7a) 


3v 

3z 


3v 

3z 


3U 

3Z 


cos  a  sin  a 


3V 

3Y 


cos  a 


sin  a  - 


3V 

3Z 


(C. 7b) 


Introducing  the  notation 


0 


3V 

3Z 


(C .  8) 
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then  the  transformation  from  local  to  global  coordinate  system  may  be 
written  in  a  matrix  form  as 


f 

r  ^ 

u 

c 

s 

u 

f 

2  Z 

V 

u 

>  = 

C  S 

< 

e 

V 

s 

-c 

3U 

vf 

3Z 

-1 

cs  -cs 

3V 

9Y 

^  / 

(C .  9) 


in  which  c  and  s  denote  cos  a  and  sin  a,  respectively.  The  trans¬ 
formation  matrix  for  the  element  is  tabulated  in  Table  3-5. 
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APPENDIX  D 


CALCULATION  OF  SECTION  PROPERTIES  AND 
STRESS  RESULTANTS  FOR  I-SECTION 


To  evaluate  the  tangent  stiffness  matrix  for  an  elastic 
element,  it  is  necessary  to  know  the  cross  section  properties.  In 
the  case  of  inelastic  analysis,  it  is  important  to  evaluate  the  trans¬ 
formed  section  in  order  to  evaluate  the  tangent  stiffness  matrix. 

Once  the  transformed  section  has  been  determined,  the  evaluation 
of  the  cross  section  properties  is  identical  to  that  of  any  arbitrary 
elastic  section. 

In  this  Appendix  the  method  of  determining  the  transformed 
section  is  discussed.  The  equations  for  numerically  evaluating 
section  properties  are  presented  and  the  equations  used  to  evaluate 
the  stress  resultants  for  an  inelastic  section  are  given. 


D. 1  Determination  of  Transformed  Section 

For  any  set  of  nodal  displacements  the  strain  at  any  point  on 
the  element,  due  to  axial  and  bending  displacements,  may  be  evaluated 
as  (Eq.  3.5.5) 


£ 

z 


u 


u 


+ 


o 


V 


+ 


t  2 

o 


r 


o 


-  y  v 


If 


1  +  U* 

o 


T  ’  2 

o 


/  1- 


V 


»  2 


+ 


2 


(D.1.1) 


-228- 


229 


Adding  the  residual  strain  at  this  particular  point  results  in  the 
total  strain  given  as 

£  _  =  £  +  £  (D.1.2) 

zt  z  R 

in  which  £  is  the  total  strain  at  any  point  in  the  z  direction,  £ 

Z  L  Z 

is  the  strain  given  by  Eq.  D.1.1,  and  £  is  the  residual  strain. 

If  the  residual  strain  distribution  varies  in  some  arbitrary 
manner,  as  shown  in  Fig.  D-la,  it  can  be  approximated  by  linear 
segments,  shown  in  Fig.  D-lb. 

For  equilibrium 

r 

Cr  •  dA  =  0  (D.1.3a) 

< 

A 

and 

* 

a  •  y  •  dA  =  0  (D.1.3b) 

R 

A 

Consider  a  plate  segment  of  the  cross  section  with  a  linear 
variation  of  residual  strains.  By  superposition,  the  total  strain 
distribution  will  also  be  linear  in  the  segment. 

The  stress-strain  curve  used  in  the  analysis  is  the  tri- 
linear  curve  shown  in  Fig.  D-2.  To  approximate  the  curve  for  mild 
steel  curve  2,  the  dashed  curve,  is  used.  It  is  apparent  from  Fig.  D-2 
that  the  strain  hardening  of  the  material  can  approximately  be  accounted 


for  in  the  analysis. 


^  \ 


(a)  Actual  Residual  Strain 


(b)  Assumed  Residual  Strain 


FIGURE  D-1  Residual  Strain  Distribution 


Stress, 
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FIGURE  D-2  Trilinear  Stress-Strain  Diagram 
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Consider  a  plate  segment  under  a  linear  variation  of  strain 

as  shown  in  Fig.  D-3.  When  the  strain  at  the  two  ends  A  and  B  are  of 

the  same  sign,  £  and  £  may  be  in  any  of  the  three  strain  ranges 

A  B 

given  below: 


£  < 

y  - 

o  < 


'B 


> 


'B 


< 


£ 

St 


£ 

y 


£  < 

y  - 

o  < 


(D. 1 . 4a) 
(D.l.Ab) 

(D. 1 . 4c) 


in  which  £  and  £  are  the  yield  strain  and  strain  hardening  strain 
y  st 

respectively,  defined  in  Fig.  D-2. 

Thus,  there  are  nine  combinations  of  strain  distribution  when 

£  and  £  are  of  the  same  sign.  Also,  there  will  be  nine  other  combin- 
A  B 

at ions  when  they  are  of  the  opposite  sign.  In  each  segment  of  the 
cross  section  the  tangent  moduli  are  determined  from  the  strain 
distribution  assuming  no  strain  reversal.  This  divides  the  plate 
segment  into  at  most  five  regions  as  shown  in  Fig.  D-3,  in  which 
moduli  are  different  (Rajasekaran  1971) . 

If  each  element  of  area  is  transformed  such  that  the  product 
of  the  current  tangent  modulus,  E^,  times  the  original  element  of 
area,  A,  is  equal  to  the  original  modulus,  E,  times  the  element  of 
transformed  area,  A*" ,  then 


E  • 

t 


A 


(D.1.5) 


. 
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Region 


3 


4 


5 


B 


(a)  Plate  Segment 


(b)  Strain  Distribution 


j 

t 


i 

I 


Segment  Length 

◄- - — ► 


(c)  Transformed  Segment 


FIGURE  D-3  Transformed  Section  of  a  Plate  Segment  Region 
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In  order  to  calculate  the  tangent  properties  it  is  sufficiently  accurate 
to  modify  the  thickness  between  the  corresponding  regions  by  the  modular 
ratio  according  to  the  form 

b C  =  E_  b  /E  (D. 1.6) 

r  t  r 

in  which  b^  is  the  transformed  thickness  of  a  particular  region  r,  and 
b^  is  the  original  thickness  of  this  region. 

Now  the  section  is  transformed  and  the  cross  section  properties 
of  the  transformed  area  can  be  evaluated  as  indicated  in  the  following 
section. 


D. 2  Evaluation  of  Cross-Section  Properties 

Once  the  transformed  thickness  of  each  region  is  evaluated, 
the  area  and  moment  of  inertia  of  each  region  are  calculated  about 
the  centroidal  axis  of  the  transformed  segment.  Moments  of  inertia 
of  each  plate  segment  are  then  transformed  to  the  global  axis  orienta¬ 
tion  of  the  cross  section.  Summing  the  properties  for  all  segments 
gives  the  centroid  and  moment  of  inertia  of  the  total  cross  section. 

The  section  properties,  A*",  1^,  I^,  1^  and  1^ ,  needed  for 
the  inelastic  analysis  are  defined,  using  Fig.  D-4,  as 


(D. 2 . 2a) 


(D. 2 . 2b) 


. 


. 


. 
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FIGURE  D-4  Cross  Section 


°6 


FIGURE  D-5  Stresses  in  a  Typical  Plate  Segment 
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y2  ‘dA*" 


(D. 2 . 2c) 


y3  ‘dA1" 


(D. 2 . 2d) 


and 


y4 ‘dA^ 


(D. 2 . 2e) 


These  properties  can  be  written  in  numerical  form  with  respect  to 
the  x  and  y  axes,  of  Fig.  D-4,  as 


n 

l  \ 

k=l 


(D. 2 . 3a) 


n 


l 


k=l 


k  k 


(D. 2 . 3b) 


n 


n 


‘2  ^  (IxA  +  ^  \  Yk 

k=l  k=l 


(D. 2 . 3c) 


X3  =  3  X  (IKK)k\+  X  \ 


k=l 


k=l 


(D. 2 . 3d) 


and 


*4  ■  X  (i4->k + 6  X  (i->k  Yk + X  ^ Y 

k=l  k=l  k=l 


(D. 2 . 3e) 


in  which  (I  ),  is  the  moment  of  inertia  of  segment  k  about  the  centroidal 
xx  k 

axis  x-x  of  this  segment  and  i s  defined  by  Eq.  B.2.2e  for  the  kth 

segment  about  x-x  axis  of  the  segment. 
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D.3  Evaluation  of  The  Stress  Resultants 

From  the  stress-strain  relationship  shown  in  Fig.  D-2,  and 
assuming  the  strain  distribution  is  known,  the  stresses  in  each  region 
of  the  plate  segment  can  be  evaluated.  The  stress  resultants  may  then 
be  calculated  numerically  in  the  x-y  reference  system  from  their 
definitions.  Referring  to  Fig.  D-5,  the  numerical  evaluation  of  the 
stress  resultants  can  be  written  in  the  following  forms  (Rajasekaran 
1971) . 


or 


also 


n  = 


a  dA 
z 


(D.3. la) 


n 


n  = 


J  l  —  l  t  (a.  +  a.  ) 

c=l  r=l  2  r  k  ir  jr 


(D.3. lb) 


m  = 


a  y  dA 
z 


(D.3. 2a) 


or 


and 


n  5  £  t. 

-  -  I  I  -hr 

k=l  r=l 


CT.  (y.  +  2y .  )  +  a.  (y.  4-  2y .  ) 

ir  jr  •'ir  jr  J  ir  Jt 


(D.3. 2b) 


m*  = 


a  y  dA 
z 


(D.3. 3a) 
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or 


ra* 


n 


I 

k=l 


5 

l 


r=l 


Z  t 
r 

12 


2 


+  a. 

jr 


(4  y2.  +2 

jr 


4- 


(D. 3 . 3b) 


in  which  r  is  the  plate  region,  n  is  the  number  of  plate  segments,  k  is 
the  plate  segment  index,  t,  is  the  plate  segment  thickness,  Z  is  the 

K.  IT 

plate  region  length  and  i  and  j  refer  to  the  two  ends  of  each  region. 
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APPENDIX  E 


INCREMENTAL  STRESS  RESULTANT  VECTORS 
FOR  INELASTIC  ELEMENT 


To  evaluate  the  element  tangent  stiffness  matrix,  for  an 

inelastic  element,  it  is  necessary  to  evaluate  the  incremental  stress 

resultants  >  and  7^—  ,  in  Eqs.  B.4.5.  In  this  appendix  the 

3qj  9qj  3q j  H 

incremental  stress  resultant  vectors  are  derived.  The  transformed 


section  properties  needed  to  evaluate  such  matrices  are  determined  in 
Appendix  D. 


E.l  Incremental  Stress  Resultant  Vectors 


Consider  variations  in  the  stress  resultants  n,  m,  and  m*  at 
any  stage  of  loading.  These  variations  can  be  written,  from  Eqs.  3.5.9, 
as 


6n 


E  6e  dA 
t 

^A 


(E.l. la) 


6m 


E  6e  y  dA 
t  J 

J  A 


and 


6m* 


E  6e  y2  dA 
t  J 

A 


(E.l. lb) 


(E.l. lc) 


in  which  E^ 


is  the  tangent  modulus  at  this  particular  stage  of  loading. 
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From  the  transformed  area  concept,  discussed  in  Sect.  D.l 


it  can  be  written  that 


E  •  dA  =  E  dA 


(E 


in  which  E  is  the  original  modulus  and  A  and  A  are  the  original 
element  of  area  and  the  element  of  transformed  area,  respectively, 
Eqs.  E.1.1  can  be  written  as 


8n 

3qj 


E  dA1 


(E 


8m 

8q, 


t 

y  dA 


8q 


(E, 


and 


8m* 

8q . 


E  y2  dAC 
3qj 


(E 


Let  the  transformed  section  properties  be  defined  by  Eqs, 
Differentiating  Eq. 3.5.5  with  respect  to  q^  and  substituting  into 
Eqs.  E.1.3  and  integrating  using  Eqs.  D.2.2  yields 


8n 

8q . 

J 


=  EA  •  b. .  -  El.  •  b_.  +  El  •  b„. 

lj  1  2j  2  3j 


(E, 


8m 

8q . 


=  El.  •  b. .  -  El.  •  b9.  +  EI„  •  b„. 

1  lj  2  2j  3  3j 


1.2) 


Thus 


1.3a) 


1.3b) 


1.3c) 


D.2.2. 


1.4a) 


(E. 1. 4b) 
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and 


3m*  t  t  t 

=  El  •  b  -  El  •  b_.  +  El,  •  b_. 
dq  2  lj  3  2j  4  3j 


(E . 1 . 4c) 


in  which 


and 


lj 


9u' 

(1  +  U?)  tt-9-  +  V? 


9v' 

o 


o  9q 


o  3q . 
J 


(E . 1 . 5a) 


9u ' 


b0  . 
2j 


=  v 


°  3q. 


+ 


v'  v" 
o  o 


»  \  2 


/ 1  -  (v 


—  I2  + 
*  )  2  l 


(  V  *  ) 

o 


1  -  (vo)2 


3v 


o 


3qj 


+ 


(v')2 

1  +  u'  +  - 2 - 

°  /T-”(v')2 


o 


9v" 
_ o 

8qj 


(E . 1 . 5b) 


3j 


v’  (v"r 

o  o 

1  -  (v')2 
o 


1  + 


(V) 

o 


M2  'l 


1-  (v  * ) 2 

O 


9v' 


o 


9q. 


+  v" 
o 


1  + 


(v1) 

o 


1  -  ( V  *  )  2 

o 


3vM 

o 

3qj 


(E.1.5c) 


For  each  gauss  point  on  the  element  the  values  uT ,  v',  and 

o  o 

v"  are  evaluated  from  Eqs.  B.3.1.  Then  substituting  Eqs.  B.3.2  into 
o 

Eqs.  E.1.4  yields 


(a)  For  1  <  j  <  4 


9n 

9q . 

J 


=  [  EAC  (1  +  up  -  ElJ  vp  (jjp 


(E.l. 6a) 
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and 


5m 

3qj 


■  [  «i  (1  +  -  ei2  vo  1  *qj 


(E.1.6b) 


5m* 

3qj 


=  [  El!:  (1  +  u')  -  EI~  v"  ]  <J)’  . 

2  o  3  o  4  j 


(E . 1 . 6c) 


(b)  For  5  <  j  <  8 


5n 

3qj 


EA11  v? 
o 


-  El. 


vf  v" 
o  o 


/  1-  (v')2 
o 


2  + 


(v») 

o 


1  -  (v;>2 


in  2 


+  El 


v’  (v") 
too 


I  \  2 


2  l-(v’)2 

o 


1  + 


(V) 

o 


1-  (v’)2 

O  '  — 1 


d>'  . 


-  El. 


(v')2 

1  +  u*  + 

/  1  -  (v ' ) 2 


+  El  v" 
2  o 


1  + 


(vf) 

o 


»  \  2 


4>a, 


(E.  1 . 6d) 


5m 

3qj 


El 5”  v * 
1  o 


-  El, 


v’  v" 
o  o 


/l-  ( V  * )  2 
o 


2  + 


1-  (v  1  )  2 
o 


v  (v")2  (  ^v;~  ^ 

t  °  o  ,  .  o 

+  EI3  1-  ( v 1 )  2  l1  1-  (vf)2 

n  v.  x  n  ✓ 


(V) 


4>n 


+ 


-  El, 


(v’)Z 

l  +  uf  + - - - 

°  /i-(v;)2 


tt 


+  EI„  v 
3  o 


1  + 


(vf ) 
o 


(E.1.6e) 
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and 


8m* 

3qj 


ei!;  v' 

2  o 


v’  v" 


t  'i  ± 


-  EI. 


/  1-  (v')2 
o 


2  + 


(vM 

o 


1  "  (vo)2J 


+  EI 


v’  (vM)2 
t  o  o 

4  1  -  (vf )2 
o 


1  + 


(V1) 

o 


K. 


+ 


-  EI. 


1  +  u' 
o 


(vo)Z 


/  1  -  (v')2 
o 


4-  EI.  v" 
4  o 


1  + 


(V;) 


t  N  * 


1  ”  (Vq)2 


ft 


(E . 1 . 6f ) 


Equations  E.1.6  are  the  incremental  stress  resultant  equations  for  a 
flexural  element,  and  can  be  written  in  a  vector  form  as  in  Table  E-l. 

For  a  truss  element  the  centroid  of  the  transformed  area 
of  any  cross  section  remains  at  the  centroid  of  the  elastic  section 
because  the  section  is  subjected  to  axial  force  only  resulting  in  uni¬ 
form  stress  distribution.  This  leads  to 


1^  =  0  (E.1.7a) 

I*  =  0  (E . 1 . 7b) 

Also  for  a  truss  element  a  linear  interpolation  function  is  used  which 
leads  to  (Eq .  B.3.1c) 
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TABLE  E-l  Incremental  Stress  Resultant  Vectors 
for  Flexural  Element 
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{— — }  is  8  x  1  column  matrix 
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{(j)'}  and  {(j)"}  are  4x1  column  matrices 
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From  Eqs.  E.1.7,  Eqs .  E.1.6  can  be  modified  for  a  truss  element  as 


(a) 

8n 

3qj 


For  1  <  j  <  2 
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(b)  For  3  <  j  <  4 
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The  column  vectors  arising  from  Eqs.  E.1.8  are  given  in  Table 
E-2.  All  section  properties  needed  to  form  these  matrices  are  given 
by  Eqs.  D.2.3.  Once  the  incremental  stress  resultant  vectors  are 
formulated,  for  an  inelastic  element,  they  are  used  in  Eq.  B.4.4  or 
Eq.  B.4.11  to  evaluate  the  inelastic  tangent  stiffness  matrix. 
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